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the value of our portfolio is 


V(t, T, S) = O) {F(t, T, S) -N (di(t, T, S) -K - N (d2(t, T, S))} 
- p(t, T) - premium 


_ p(t, T) — pt, S) 
= Fars) COTO NGET, -K-N (dal T.S) 


— p(t, T) - premium 


where F(t, T, S) is the forward swap rate between t = T and t = S at 

time ¢ and 

In| Fee! + $0, 7,5)- (T=) 
o(t,T,S)-/T-t 

dy(t,T, S) =d; —o(t,T,S)- VT -t 


d\(t,T,S) = 


19.1.6 The Normal Black Model 


Usually the underlying security is assumed to follow a lognormal 
process (or Geometric Brownian Motion). However, there are some 
traders who believe that the normal process describes the real market 
more closely than that of lognormal counterpart. This model is also 
known as the Bachelier’s model. 

Let us assume that the current future price, strike price, risk-free 
interest rate, volatility and time to maturity as denoted as f, K,r,o 
and T -t respectively. Let us also assume that the current future price 
follows the following normal process: 


df = udt + odW;, 


where n is a constant drift. For instruments like swaptions, f repres- 
ents the forward rate. 

Let us start by studying the behaviour of the delta-hedged portfo- 
lio which consists of long delta shares of future contract and short 
one derivative in question. Say, call it II. Let us also denote the value 
of derivative by g. Then, the value of the delta-hedged portfolio is 
given by: 


Hege r 
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So applying Ito’s lemma using the earlier mentioned SDE into the 
changes of the portfolio value, one gets 


a 
dl = dg - 1 


3 dg 1 402 ə ə 
= (E +u +0 \dt+ 0o Saw -Ë (udt + caw) 
ðt af 2 af2 of of 


We want the aforementioned quantity to be a Q-martingale under the 
discounted expectation with risk-free rate. This is the same as stating 
that the previous quantity equals the gain from the risk-free interest 
rate for the portfolio value. So, we have: 


ail = rIldt 


Since it costs nothing to enter into a futures contract, one has II = g. 
Thus, we obtain the following PDE 


dg 1 407 
ue + o 8 = rg 
ot 2 df? 
In a risk-neutral world the process is giving as 
af = odV; 


with the trivial solution, from integration over the interval [t, T]: 
ITD) =f O +o (Vr- V 


We see that f is a Gaussian process; N[f;, o*(T —1)], that is, with mean 
f(t) and variance o?(T — t). By the application of Feynman-Kaé, we 
obtain the following solution 

T- F rA 


(T) - e 220- d 
na eal fr 
=00 


—r(T-t) 
= [fr =f,tovT t-z] A 


T (vT 
x [ome 20°71) oT — tdz 


—0O 


g(t, fr) = e” TPE? [®(T)] = 


eT- P 


ot 
= —_— | (T). edz 
J 20 
—0o 
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Here 
(fr -—K)*  foraCall 


a= | (K—fr)* foraPut 


For the Call we have 
rT F 2 
Mc(t) = J fr-K) -€7dz 
~y 210 EA 
-r(T-1) j a 
=ë f+o[T-t-z-KY -e 7dz 
270 a. 
orl) F 2 
= (f; tovT—t-z-K)-e€ 7dz=A-B 
200 | 
zo 
Set f, = F and with zo = = we get 


gnm (2 r(T-1) 
A= (F-K) | edze" ~F — K) - N[zo] 
~v 20T 
zo 
gov T 


erT) P 2 i 
B= ovVvT-t-z-€ 2dz=e 
J 20 : J 20 
0 


Then, the fair values of call C (payer swaption, PS) and put P (receiver 


iS) 


tl 


swaption, RS) are expressed as 
VT -t 
C= eT | (F-K): N(d) + Xe €? 
~v 2T 


Ow T-t zi) 
— e 


P = e” | (K - F) - N(-d) + 
V2 
where 
d= F-K 
oV/T-t 
and 
t = Tenor of swap in years (time between swaption maturity and 


swap maturity). 
F = Forward rate of the underlying swap. 
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K = Strike rate of the swaption. 

r = Risk-free interest rate. 

T = Time to swaption expiration in years. 

o = Volatility of the forward-starting swap rate. 
m = Compounding’s per year in swap rate. 


To apply this on swaptions we need, as before, to multiply C and P 
with the annuity 


1 
l- (+F/m) 
F 


The Greeks can easily be calculated by simple differentiations 


ac 


Ac= — =e”. Nad 

C= ap =e (d) 
aP 

Ap = — = e70 . N(-d 

pep ~ (-d) 


3C 8P eT 1 en 
= — . e 

ƏF? ƏF? oyT-t 2r 

W 22” gral ee 

do 00 J 20 


ac 
Ora eae 
Oa 


JIT = 
x (e-no? n T etn) 


V2- m 2V2- m-t 


Op = oP apr) 
ot 


F= 
x (e-on py a etn) 


v2- m 2V2- m-t 


19.1.6.1 Convexity Adjustments 


A standard bond or interest rate swap has a convex price-yield relation- 
ship. To price options with the Black-76 model when the underlying 
asset is a derivative security with a payoff function linear in the bond 
or swap yield, the yield should be adjusted for lack of convexity value. 
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Examples of derivatives where the payoff is a linear function of the 
bond or swap yield are constant maturity swaps (CMS) and constant 
maturity treasury swaps (CMT). The closed-form formula published by 
Brotherton-Ratcliffe and Iben (1993) assumes that the forward yield is 
lognormal distributed. 


1 ay. 
adj. = 5 aE We (es - 1) 


where 


P = Bond or fixed side swap value. 
yr = Forward yield. 

T = Time to payment date in years. 
o = Volatility of the forward yield. 


Example 19.1.3 
Consider a derivative instrument with a single payment 5 years from now that is 
based on the notional principal times the yield of a standard 4-year swap with an- 
nual payments. The forward yield of the 4-year swap starting 5 years in the future 
and ending 9 years in the future is 7%. The volatility of the forward swap yield is 
18%. Calculate the convexity adjustment of the swap yield. 

The value of the fixed side of the swap with annual yield is equal to the value 
of a bond where the coupon is equal to the forward swap rate/yield yr. 

c c c lt+c 


= T zt 3°" 4 
I+yr (l+ypy (+y (+yr) 


The partial derivative of the swap with respect to the yield is 
oP c 2c 3c 4(1 +c) 
Oye +y (tyr (+y (+yr) 
= {c = yp = 0.07} = -3.3872 


and the second partial derivative with respect to the forward swap rate is 
a?P 2c 6c 12c 20(1 +c) 
T 7t gF 5+ 6 
dye (tyry (Ut+yry (QA+yrP (1+yr) 
The convexity adjustment can now be found as: 
115.2933 7 gigs 
dj. =-= - 0.07 =° —1) = 0.001 
adj. = 7 agq9 ‘007 (¢ ) =0.0019 
The convexity-adjusted rate is then equal to 7.19% (0.07+0.0019). 


= 15.2933 
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19.1.6.2 Vega of the Convexity Adjustment 


By taking the derivative of the convexity adjustment, we get the 
convexity adjustment’s sensitivity to a small change in volatility 


v=- -yeoTe® T 
19.1.6.3 Implied Volatility From the Convexity Value 
in a Bond 
If the convexity adjustment is known, it is possible to calculate the 


implied volatility by simply rearranging the convexity adjustment 
formula 


dj. 1 
a +1 


afer fae). 2 a 
2 aye OYF YF 


19.1.7 European Short-Term Bond Options 


o= |in 


European bond options can be priced in the Black-76 model by using 
the forward price of the bond at expiration as the underlying asset. 


c=e' (F-N(d})—K - N(d2)) 
p = e” (K - N(-d2) -F - N(-d\)) 


where F is the forward price of the bond at the expiration of the 
option, and 


PS In(F/K) —(07/2)T _ 
2 = oVT = 


This model does not take into consideration the pull to par effect of the 
bond. At maturity, the bond price must be equal to principal plus the 
coupon. For this reason, the uncertainty of a bond will first increase 
and then decrease. 

The Black-76 model assumes that the uncertainty (variance) of the 
underlying asset increases linearly with time to maturity. Pricing of 


dı oVT 
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European bond options using this approach should thus be limited to 
options with short time to maturity relative to the time to maturity of 
the bond. A rule of thumb used by some traders is that the time to 
maturity of the option should be no longer than one-fifth of the time 
to maturity on the underlying bond 


Example 19.1.4 
Consider a European put option with 6 months to expiry and strike price 122 ona 
bond with forward price at option expiration equal to 122.5. The volatility of the 
forward price is 4%, and the risk-free discount rate is 5%. Calculate the option’s 
value. 

F = 122.5, K = 122,T =0.5,r = 0.05, o = 0.04. 


_ In(122.5/122) — (0.042/2) - 0.5 
: 0.0405 


N(-d\) = 0.4369, N(-d2) = 0.4481 


= 0.1587, dy = dı —0.04,/0.5 = 0.1305 
p = e995 [122N(-d>) — 122.5N(—-d1)] = 1.1155 


19.1.8 The Schaefer and Schwartz Model 


Schaefer and Schwartz (1987) modified the Black-Scholes model for 
pricing bond options to take into consideration that the price volatility 
of a bond increases with duration 


c =S - N(d1) - Ke”T N(d) 
p = Ke N(-d2) -S - N(-d}) 


where o = (as?!) D, and 


In(S/K) — (r + 07/2)T 
_ In(S/K) — ( WT di-oVT 
o/T 
Here D is the duration of the bond after the option expires. A is estim- 


ated from the observed price volatility og of the bond. « is a constant 
that Schaefer and Schwartz suggest should be set to 0.5. 


2 


pa 
go-l . D* 


where D* is the duration of the bond today. 
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Example 19.1.5 


Assume that the duration of the bond is 8 years and that the observed price 
volatility of the bond is 12%. This gives: à = 0.15. 

In Table 19.1 we use this value and compares the option prices from the Schaefer 
and Schwartz formula with option prices from the Black-76 formula. 


Table 19.1 Option prices from Schaefer and Schwartz and Black-76 


Bond Base Adjusted Black-76 Modified 
Duration Volatility (%) Volatility (%) Value Black-76 value 
1 12.0 1.5 5.5364 0.6929 

2 12.0 3.0 5.5364 1.3857 

3 12.0 4.5 5.5364 2.0783 

4 12.0 6.0 5.5364 2.7707 

5 12.0 7.5 5.5364 3.4628 

6 12.0 9.0 5.5364 4.1545 

7 12.0 10.5 5.5364 4.8457 

8 12.0 12.0 5.5364 5.5364 
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Convertibles 


20.1 Convertible Bonds 


A convertible bond is a security issued by a company that may be 
converted from debt to equity (and vice versa) at various prices and 
stages in the life cycle of the contract (e.g. the time to maturity). There 
are many types of convertible bonds with various conversion proper- 
ties and complex structures. Common examples of convertible bonds 
are Convertible Preferred Stock bonds, Zero-Coupon convert- 
ibles, Mandatory convertibles, to name but a few. The traditional 
(simplest) convertible bond is one that is a fixed coupon paying bond 
when the stock price S is below some predetermined conversion price 
K Ge. S < K), and may be converted to a predetermined number 
of stocks (the conversion ratio) when above (i.e. S > K). It may im- 
mediately become apparent that the traditional convertible bond is 
effectively a bond with an imbedded call warrant (option) on the stock 
of the issuing company with strike price K. 

Since a convertible bond is a liability of the issuer, if the company 
goes into liquidation, the convertible bondholder has priority over 
most other parties except pure bondholders. 

Returning to our traditional convertible bond example, for the mo- 
ment we ignore credit rating issues, the convertible price Beon is close 
to its fixed income value (bond) B of an equivalent pure bond from 
the same issuer when deep out-of-the-money. Hence, 


Bcn ~B forS < K. 
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When the stock price is very high and exceeds the conversion price 
K, the convertible bond becomes stock-like when deep in-the-money. 
Hence, 


Beon ~S forS >K. 


Moreover, given that the convertible is effectively a bond plus a war- 
rant on the stock, the price of a convertible must be comparable to the 
sum of the two individual components. From arbitrage considerations 
for options, we know that the price of a call warrant W must satisfy 


W > max[S- K, 0] 


prior to expiry. Ignoring issues related to credit risk, when the warrant 
element is out-of-the-money the convertible is worth the fixed income 
value B. Therefore, the convertible has to satisfy the relation 


Bcon > max[S- K, 0] + B 


given the hybrid nature of the convertible bond. Here the price of the 
convertible bond has been approximated as 


Beon = B+W 


In Fig. 20.1 we show the price track of a traditional convertible bond. 
The convertible price (solid red line) is always above the stock price 
track as already discussed. The red dotted line demonstrates the price 
track of a bond plus warrant in the absence of credit risk. 

The true price track of the convertible is found to fall below the 
fixed income value (Bcon < B) when the stock price falls to low levels 
due to the widening of credit risk spread when the company’s stock 
price falls. As the stock price falls to low levels there is an increas- 
ing correlation between the price of the convertible and the stock; 
however, as the price rises to very high values, the correlation be- 
comes insignificant (as the credit rating improves). A fall below the 
fixed income value is also seen in the pure bonds due to poorer credit 
rating. 

The following terms are used to describe various sections of the 
price track of the convertible: 


1. Distressed debt -In this region the convertible is on close to a de- 
fault event. If a default event occurs, a sum proportional to the 
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Fig. 20.1 The price track of a convertible bond 


recovery rate R is paid out to the holder of a convertible. The 
value of the convertible is highly sensitive to the credit risk spread 
(a parameter often referred to as omicron) in this region. 


2. Busted convertible-A term often used to describe a convertible that 
is out-of-the money but above the distressed zone. 


3. Hybrid zone-The convertible shows behaviour between the stock 
and a pure bond. 


4. Equity zone-The convertible price is more equity-like than debt. 
Credit risk factors become insignificant since the company’s credit 
rating is high due to the high stock value. 


The hybrid nature of convertibles is often exploited by arbitrageurs in 
the rather popular practice of convertible arbitrage. Convertible arbit- 
rage is especially successful at times of high volatility in stock price, 
producing high returns with relatively low risk. Delta, gamma and 
other more sophisticated hedging strategies are used to capture the 
low-risk profits (though not completely risk-free). Instances can occur 


528 J.R.M. Roman 


when the equity options imbedded are priced different to those of the 
equivalent pure options that may exist in the market for various reas- 
ons. These situations, when a relative price difference is observed, are 
exploited by arbitrageurs in a long-short trade. 


20.1.1 A Model for Convertibles 
The stock is modelled as a lognormal Brownian process 
dS = uSdt + oSdZ, 
and the interest rate as 
dr = u(r, t)dt + w(r, t)dZ> 


where Z; and Z2 are two independent Wiener processes with a cor- 
relation p. The drift u(r,t) and volatility w(r,t) is dependent on the 
interest rate model. We make the choice 


dr = (a, — by r)dt + wdZy 


where all parameters are time dependent. The value of the convertible 
V depends on the stock price, the interest rate and of time, V(S, 7,1), 
and given by the following PDE 


ƏV oo 282 a2V 3V ow 02V ƏV aV 
+ + pa Sw + +r +(u—Aw) rV =0 
ðt 2 as2 aSor 2 ar as ðr 


This is found by hedging the two processes against each other and 
with the introduction of the market price of risk. The value of the 
convertible must be: V(S,r,t) > nS where n is the number of stocks 
on exercise. We get 


ƏV o?S2?3?V 32V w283? V av aV 
— + ——~ + po Sw——— + — — + rS— + (a, —bir)— -rV =0 
ðt 2 as aSér 2 dr as or 


This equation can be solved with a finite difference method. 


21 


A New Framework 


21.1 Pricing Before and After the Crisis 


Ten years ago if you had suggested that a sophisticated investment 
bank did not know how to value a plain vanilla interest rate swap, 
people would have laughed at you. But that isn’t too far from the case 
today. 

We will now give an introduction to yield curve constructions and 
how this has been changed since after the financial crisis. 


21.1.1 Introduction 


Pricing complex interest rate derivatives requires modelling the fu- 
ture dynamics of the yield curve term structure. Most of the literature 
considers the existence of the current zero-coupon yield curve as 
a given, and its construction is often neglected, or even obscured, 
as it is considered to be more of an art than science. Actually, any 
yield curve term structure modelling approach will fail to produce 
good/reasonable prices if the current term structure is not correct. 
Financial institutions, software houses and practitioners have de- 
veloped their own proprietary methodologies in order to extract the 
current zero-coupon yield curve term structure from quoted prices 
on subsets of liquid market instruments. These can be divided into 
two groups: “best fit” and “exact fit” algorithms. “Best-fit” algorithms 
start by assuming a functional form for the term structure and calibrate 
its parameters such as to minimize the re-pricing error of the chosen 
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set of calibration instruments. An example of this is the Smith-Wilson 
approach to discounting used by regulators for insurance companies. 

In banks, “exact-fit” algorithms are often preferred in practice. Such 
algorithms fix the zero-coupon yield curve on a time grid of N points 
in order to exactly re-price N pre-selected market instruments, often 
referred to as benchmarks. The implementation of these algorithms is 
often incremental, extending the yield curve step by step in increas- 
ing order of maturity for the selected instruments, in a “bootstrap” 
approach. 

Intermediate yield curve values are then obtained by interpolation 
on the bootstrapping grid. Here different interpolation algorithms are 
available but little attention has been devoted in the literature to the 
fact that interpolation is often already used during bootstrapping, 
not just after that, and that the interaction between bootstrapping 
and interpolation can be subtle if not nasty (see e.g. Patrick Hagan, 
Interpolation method for curve constructions). 

While naive algorithms may fail to deal with market subtleties such 
as date conventions, the intra-day fixing of the first floating payment 
of a swap, the turn-of-year effect, the futures convexity adjustment, 
etc., even very sophisticated algorithms used in a naive way may fail to 
provide relevant estimates of forward Euribor rates in difficult market 
conditions, such as those observed since the summer of 2007 and the 
so-called subprime credit crunch crisis. Today using just one single 
curve is not enough to account for forward rates of different tenors, 
such as 1, 3, 6, 12 months, because of the large basis swap spreads 
presently quoted in the market. 

Prior to the credit crisis, (zero-coupon) yield curve modelling was 
reasonably well understood. The underlying fundamental principles 
had existed for over 15 years with steady evolutions in areas that 
were most relevant to options and complex products. Credit and li- 
quidity issues were ignored as their effects were minimal. Pricing a 
single-currency interest rate swap was relatively straightforward. A 
single “default free” zero-coupon yield curve was calibrated to liquid 
market products, and future cash flows of other instruments were 
discounted and evaluated using this single curve. There was little vari- 
ation between implementations, and results across the market were 
consistent. 

Following the credit crisis, yield curve modelling has undergone 
nothing short of a revolution. During the credit crisis, credit and liquid- 
ity problems appeared in several markets. This drove apart previously 
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closely-related interest rates. For example, Euribor basis swap spreads 
dramatically increased and the spreads between Euribor and Eonia 
overnight indexed swaps (OIS) diverged. In addition, the impact of 
counterparty credit on valuation and risk management dramatically 
increased. 

Existing modelling and infrastructure no longer seemed to work 
and a re-assessment of the basic principles has taken place. Currently 
a new interest rate modelling framework is evolving which is based 
on OIS discounting and the integration of credit value adjustments 
(CVAs). Pricing a single-currency interest rate swap ARS) now takes 
into account the difference between projected rates such as Euribor 
that include credit risk and the rates appropriate for discounting cash 
flows that are risk free or based on funding cost. CVAs take into ac- 
count the likelihood that the counterparty will default, along with 
the expected exposure given default, the volatility of these expected 
exposures, and wrong way risk.! 

One of the many consequences of the liquidity crisis that started 
in the second half of 2007 has been a strong increase in quoted 
basis spreads in the market between single-currency interest rate in- 
struments, in particular for swaps. This is characterized by different 
underlying rate tenors (e.g. Euribor3M, Euribor6M, etc.), reflecting 
the increased liquidity risk, and an increased preference of financial 
institutions for receiving payments with higher frequency (quarterly 
instead of semi-annually, for instance). Such asymmetry has induced 
the “segmentation” of interest rate markets into various sub-areas, with 
1M, 3M, 6M and 12M underlying rate tenors. Each area is character- 
ized, in principle, by its own internal dynamic, reflecting the different 
views and interests of the market participants. 

In order to price derivatives we must now know 


1. the underlying interest rate (Libor/Eonia, etc.); 


2. the matching market prices for plain vanilla derivatives, which are 
used to construct zero-coupon yield curves and the term structure 
of volatilities, as well as for calibration and hedging; 


3. the transaction mechanics: collateral and liquidity/funding issues; 


4. the counterparties: credit and default issues. 


1 Wrong-way risk is defined by the International swaps and Derivatives Association (ISDA) as 
the risk that occurs when “exposure to a counterparty is adversely correlated with the credit 
quality of that counterparty”. 
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21.1.2 After the Crises - How the Market Has Changed 


In Fig. 21.1, we show the market quotations on 31 September 2008, 
for six basis swap curves corresponding to the four Euribor tenors 1, 
3, 6 and 12 months. 

As one can see, the basis spreads are monotonically decreasing from 
over 100 to around four basis points. There is neither way nor any 
good reason to ignore such quotations in a market pricing framework 
of interest rate derivatives. Before the credit crunch of August 2007, 
the basis spreads were just a few basis points. 

The consequences of the credit crunch crisis that started in August 
2007 can also be seen in the historical series of the Euribor 6-month 
(6M) rate vs. the Eonia OIS rate as in Fig. 21.2. The change in the 
spread is essentially a consequence of the different credit and liquid- 
ity risk reflected by Euribor and Eonia rates. This divergence is not 
a consequence of any difference in counterparty risk between these 
contracts. Euribor and OIS rates are exchanged in the interbank mar- 
ket by risky counterparties, but they also depend on different fixing 
frequences of the underlying rates. 


100 
— - Euribor 1M vs 3M 
-- Euribor 1M vs 6M 
30 Tor — Euribor IM vs 12M 
; — Euribor 3M vs 6M 
60 Euribor 3M vs 12M 
===- Euribor 6M vs 12M 
40 
20 


Fig. 21.1 A typical overnight index swap 
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Fig. 21.2 The 6-month Euribor vs. Eonia overnight indexed swap rate 


In Fig. 21.3, we show the basis spread between some of the term 
structures of interest rates in June 2011 

The different influence of credit risk on Libor and overnight rates 
is shown in Fig. 21.4, where we compare the historical series for 
the Euribor-OIS spread above credit default swaps (CDS) spreads for 
some main banks in the Euribor Contribution Panel (Commerzbank, 
Deutsche Bank, Barclays, Santander, Royal Bank of Scotland and 
Credit Suisse). We observe that the Euribor-OIS basis explosion in Au- 
gust 2007 exactly matches the CDS explosion, corresponding to the 
generalized increase of the default risk seen in the interbank market. 

An effect of the credit crunch has been the great increase of col- 
lateral agreements (Credit Support Annex (CSA) agreements) in an 
attempt to reduce the counterparty risk of over-the-counter (OTC) 
derivatives. Nowadays most of the counterparties in the interbank mar- 
ket have mutual collateral agreements. In 2010, more than 70% of all 
OTC derivatives transactions were collateralized.? 

The main feature of the CSA is a margination mechanism similar 
to those adopted by central clearing houses for standard instruments 
exchange (e.g. futures). Shortly, at every margination date, the two 


2 International swaps and Derivatives Association, ISDA (2010). 
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Fig. 21.4 The 5 months Euribir-OIS spread and credit default spread for some main 
banks in Euorope during the financial crisis 
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counterparties check the value of the portfolio of mutual OTC trans- 
actions and regulate the margin, adding to or subtracting from the 
collateral account the corresponding mark to market variation with 
respect to the preceding margination date. 

We can also look at CSA as a hedging mechanism, where the 
collateral amount hedges the creditor against the event of default 
of the debtor. The most diffused CSA provides a daily margination 
mechanism and an overnight collateral rate. 

An important consequence of the diffusion of collateral agreements 
among the interbank counterparties is that we can consider the prices 
of derivatives quoted in the interbank market as counterparty risk-free 
OTC transactions. A second important consequence is that, by no ar- 
bitrage, the CSA margination rate and the discounting rate of future 
cash flows must match, hence the name “CSA discounting”. Since the 
OIS curve is what is usually used, the alternative name is “OIS dis- 
counting”. Such a discounting curve is also the best available proxy 
for a risk-free yield curve. 


21.1.2.1 Discount Rate — The Risk-Free Interest Rate, Then 
and Now 


Before the crisis, the risk-free rate was defined in either two ways: 


1. The rate at which the government borrows. This assumes that 
governments do not default and that their bond yields are not 
contaminated by liquidity or tax premiums. 


2. The (Libor) rate at which the big international (Libor rated) banks 
borrow at short maturities (e.g. 3 months). This assumes that 
the credit spread over a short horizon of a highly rated bank is 
practically zero. 


Now, many actors in the markets use the OIS/Eonia as the risk-free 
interest rate. This is obtained by compounding the overnight interest 
rate at which Libor rated banks borrow. One way to justify this new 
definition is to argue that zero credit risk at short horizons is valid for 1 
day but not for longer periods. Another justification is that OIS/Eonia is 
the rate that is paid on cash collateral and is therefore the correct dis- 
counting rate for collateralized forward contracts and swaps. The latter 
argument does not require that Libor rated banks have zero credit risk 
at even overnight horizons. 
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21.1.2.2 Reasons to Change the Idea of the Risk-Free Rates 


The argument that highly rated entities have zero credit risk over short 
time horizons requires perfect observability of the balance sheet (as- 
set values and debt). However, the crisis in 2008 showed that bank 
balance sheets are terribly opaque. Banks valued illiquid assets not by 
marking to market but by marking to model (see model risk) or some- 
times marking to their own belief. Similarly, the true liabilities of the 
bank were hidden using off balance sheets. This meant that the true 
distance to default could be much lower than the estimated distance 
to default based on observable parameters. 

The unreliability of accounting information creates a “jump to de- 
fault” risk and creates non-negligible credit spreads over short time 
horizons. So, even 1-month Libor can no longer be regarded as risk 
free. Whether overnight Libor can be regarded as risk free is an open 
question. 

We also have to remember, even if a bank is Libor rated today (i.e. it 
can borrow at Libor today), this does not guarantee that it can borrow 
at Libor in the next quarter because it may no longer be Libor rated 
then. The Libor rate for the next quarter will be the rate at which a 
bank can borrow if it is Libor rated on that date. This makes a big 
difference between 6-months Libor and the 3-months Libor rate. 

This implies that an interest rate swap whose floating leg is 3-months 
Libor is not the same as a swap whose floating leg is 6-months Libor. 
The floating leg payments on the first swap are expected to be lower 
and therefore the fixed leg should also be lower. A tenor swap in 
which both legs are floating - say 3-months Libor on one leg and 6- 
months Libor on the other leg - should include a tenor spread? on the 
3-months leg to ensure that the swap is fair at inception. We are assum- 
ing here that the interest rate swaps are free of credit (counterparty) 
risk as discussed later when we turn to collateralization. 

Similar explanations can be given for the cross-currency swap (CCS) 
spread. In a basis CCS, both legs are floating but in different currencies. 
For example, US dollar Libor on a dollar notional amount may be ex- 
changed for Japanese yen Libor on an equivalent yen notional with an 
exchange of principals at the end. If Libor is regarded as risk free, then 


3 This tenor spread is usually called basis spread, but I think tenor spread is a better name. In 
the following I will call this spread basis spread 
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the two legs are floating rate bonds that must be worth par and the 
swap should trade flat. In reality however, for several years, there has 
been a premium on the yen leg of this swap. Credit quality differences 
could be one explanation for this, though liquidity issues might also 
play an important role. 


21.1.3 A Multi-Curve Framework 


To set up this framework, we denote with My, x = {d, fi, ..., fa} a 
multiple distinct interest rate sub-market, characterized by the same 
currency and by distinct money market accounts Bx, such as 


t 
B,(t) = exp J ry(t’)dt' } , 
0 


where r(t) are the associated short rates. We also have multiple yield 
curves c in the form of a continuous term structure of discount 
factors, 


Cy = {T > px(to, T), T > to}, 


where fo is the reference date of the curves (e.g. settlement date, or 
today) and p,(t, T) the price at time t > fo of the M,-zero-coupon 
bond for maturity T, such that p,(7, T) = 1. In each sub-market My we 
postulate the usual no-arbitrage relation, 


Px(t, T2) = px(t, Tı)px(t, Ti, 79), t S Ti < To, 


where p(t, Tı, T2) denotes the M, forward discount factor from time 
Tı to time T2, prevailing at time ¢. The financial meaning is that, in 
each market M,, given a cash flow of one unit of currency, at time T>, 
its corresponding value at time t < To must be unique. This must be 
true, both if we discount in one single step from T, to t, using the 
discount factor p(t, T2), and if we discount in two steps, first from 7> 
to Tı, using the forward discount p,(t, Tı, T2) and then from T to t, 
using p(t, T1). 

We also define continuously compounded zero-coupon rates 
Zx(to, T) and simply compounded instantaneous forward rates fx(to, T) 
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such that 


Px(to, T) = exp {—zx(to, T) . Tc(to, T)} E exp ~ fi u)du ry 


In px(to, T) = —2;(t0, T) . Te(to, T) =— [ fs u)du, 


where te(T1, T2) = t(71, T2, dcc) and dce, the day count convention for 
the zero rate. From the relationships mentioned earlier it is immediate 
to observe that 


e zx(to, T) is the average of fx(to, u) over [to, T]; 

e if rates are non-negative, p(toọo, T) is a monotone non-increasing 
function of T such that O < p(to, T) < 1VT > tọ; the instantan- 
eous forward curve Cy is the most severe indicator of yield curve 
smoothness, since anything else is obtained through its integration, 
therefore being smoother by construction. 


Using the rate formulas we define two other curves associated with 
C,, a zero curve and an instantaneous forward rate curve 


CZ = {x(to, T), T 2 to} 
Cl = {x(to, T), T > to} 


where zy and fy are given from the previous equation as 


1 
to, T) = -———~ - Inpx(to, T 
Zx(to, T) LoT) npx(to, T) 


a) 
fix(to, u) = -zz Inpxto. Dli=r = ixllo, T) + telto, Da “cal, Dh. 


In the following we will denote with Cy the generic curve and we 
will specify the particular typology (discount, zero or forward curve) 
if necessary. 

Denoting with F(t; Tı, T2) the discretely compounded forward 
rate corresponding to the M, forward discount factor p,(t, Ti, T2), 
resetting at time Tı and covering the time interval [T1, T2], we have 


px(t,T>) _ l 


Px(t, Ti, T2) = = —— 
” P(t, Tı) 1+ F(t; T1, T2) - tx(T1, T2) 
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where 1,(71, T2) is the year fraction between times Tı and T2 with 
day-count dc,. We then obtain the familiar no-arbitrage expression 


1 | 1 1 _ Px(t, T1) — px(t, T2) 


F,(t; TıT2) = = . 
Tx(T], T2) Px(t, Ti, T2) Tx(T], T2) : Px(t, Tə) 


This can be also derived as the fair value condition at time t of the 
forward rate agreement (FRA) contract with pay-off at maturity T2 
given by 


FRA,(T2;T, T2, K, N) =N - (T1, T2) (Lx(11, T2)- K] 


L (Ti; T2) = ’ 
i Tx(T1, T2) - px(T1, T2) 


where N is the nominal amount, L,(7T,, T2) the Tı-spot Libor rate for 
maturity T) and K the (simple compounded) strike rate (sharing the 
same day-count convention for simplicity). Introducing expectations 
we have Vt < Ti < To: 


T 
FRA,(T9; T1, T2, K, N) = pet, T2) - Ege? [FRA(T2; T1, T2, K, N)], 
T: 
= N - palt, To) tT, Ta) [ER ILT, Ta) Kf, 


where the expectation is taken in the M, — T2-forward measure corres- 
ponding to the numeraire px(t, T2), at time t with respect to measure 
Q and filtration F;, encoding the market information available up to 
time ¢, and we have assumed the standard martingale property of the 
forward rates 


T: T: 
F(t; T1, T2) = EL” [Fy(11; T1, T] = EL” ILT, T2)] 


to hold in each interest rate market My. The previous assumptions im- 

ply that each sub-market My» is internally consistent and has the same 

properties as the “classical” interest rate market had before the crisis. 
The value of a FRA at time t can therefore be written as 


x(t, T1) — px(t, T2) 
i = wT T2) -K | - pxl(t, T2), 
Dx(t, T2) 


=N: [px(t,T1)- (1 + t(T1, T2)) - K - px(t, T2)] . 


FRA,(¢;71,T2,K,N)=N- | 
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Regarding swap rates, given two increasing dates vectors T = 
{To, .-., In}, S = {So, ..., Sm}, Tn = Sm > To = So È to, and an 
interest rate swap with a floating leg paying at times S;, j = 1,...,m, 
the Euribor rate with tenor [S;_1; Sj] fixed at time S;_;, plus a fixed leg 
paying a fixed rate at times 7;, i = 1,...,n, the corresponding simple 
compounded fair swap rate on curve Cy with day count convention 
dcs is given by 


m 

> plt, S)tr(S)-1, S)F x(t, S;-1, S)) 

a l t, Tı E C t, T, 
sarg px(t, To) — px(t, Tn) 


A(t, T) A(t, T) , 


where 


n 
Ax(t,T) = F px(t, Tes TTo 
i=l 


is the forward annuity on curve C, and we have defined ts(Ti-1, Ti) = 
t(Ti1, Tj, dcs). Notice that on the right-hand side of the previous equa- 
tion for S,, we have used the definition of the forward rate and the 
telescopic property of the summation. Actually the telescopic prop- 
erty would hold exactly only if the forward rates end dates equal the 
next forward rate start dates, with no period’s gaps or overlaps. 

This is not true in general, because start and end dates are adjusted 
with their business day convention, and the resulting periods do not 
concatenate exactly. Typically, such date mismatch does not exceed 
one business day (which sometimes can be three calendar days). In 
practice, on the one hand the error is small, of the order of 0.1 basis 
points; on the other hand nothing prevents us from using the exact 
dates and accrual periods. 


21.1.3.1 The Single-Curve Framework 


So how did the market practice for pricing and hedging interest 
rate derivatives change through the credit crunch crisis? Using the 
notation described earlier we start by considering a general single- 
currency interest rate derivative with m future coupons with pay-offs 
T = {T, ..., Mm}, with m; = mi(F,), generating m expected cash 
flows c = {c1, ..., Cm} on the future dates T = {7T), ..., Tm}, with 
ta Ti <e < TIm 
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The pre-crisis standard market practice was based on a single-curve 
procedure, well known to the financial world, that can be summarized 
as follows: 


1. Select one finite set of the most convenient (e.g. liquid) vanilla 
interest rate instruments traded in real time in the market with in- 
creasing maturities. For instance, a very common choice in the EUR 
market was a combination of short-term EUR deposit, medium-term 
futures on Euribor3M and medium-long-term swaps on Euribor6M. 


2. Build one yield curve, Cq using the selected instruments plus 
a set of bootstrapping rules (e.g. pillars, priorities, interpolation 
method, etc.) 


3. Compute the relevant forward rates using the same yield curve 
Cq as 


Patt, T;-1) — Patt, Ti) 
Ta(Ti-1, Ti) < palt, T 
where t < Tia < Tj, i€ {1, ..., m} 


Falt, Ti-1, Ti) = 


4. Compute cash flows c; as expectations at time t of the correspond- 
ing coupon pay-offs n;(F4) with respect to the 7;-forward measure 
associated to the numeraire pa(t, T;) from the same yield curve Cg, 


Tj 
c= c(t; Tun) = EY” Ea: 


5. Compute the relevant discount factors pg(t, Ti) from the same yield 
curve C4. 


6. Compute the derivative’s price at time t as the sum of the discoun- 
ted cash flows, 


m m 


Ti 
mi Tord = > patt,Ti)- e Teny =) pat, Ti) By ’[xi(Fa)l. 
i=l i=l 


7. Compute the delta sensitivity with respect to the yield curve Cg and 
hedge the resulting delta risk using the suggested amounts (hedge 
ratios) of the same set of vanillas. 


542 J.R.M. Roman 
21.1.3.2 The Multi-Curve Framework 


Unfortunately, the pre-crisis approach outlined earlier is no longer 
consistent, at least in its simple formulation, with the present market 
conditions. 


e First, it does not take into account the market information carried 
by basis swap spreads, now much larger than in the past and no 
longer negligible. 

e Second, it does not take into account that the interest rate market 
is segmented into sub-areas corresponding to instruments with dis- 
tinct underlying rate tenors, characterized, in principle, by different 
dynamics (e.g. short rate processes). Thus, pricing and hedging an 
interest rate derivative on a single yield curve mixing different un- 
derlying rate tenors can lead to “dirty” results, incorporating the 
different dynamics, and eventually the inconsistencies, of distinct 
market areas, making prices and hedge ratios less stable and more 
difficult to interpret. On the other side, the more the vanillas and 
the derivative share the same homogeneous underlying rate, the 
better should be the relative pricing and the hedging. 


e Third, by no arbitrage, discounting must be unique: two identical 
future cash flows of whatever origin must display the same present 
value (PV); hence we need a unique discounting curve. 


In principle, a consistent credit and liquidity theory would be required 
to account for the interest rate market segmentation. This would also 
explain the reason why the asymmetries cited earlier do not neces- 
sarily lead to arbitrage opportunities, once counterparty and liquidity 
risks are taken into account. Unfortunately such a framework is not 
easy to construct. In practice, interest rate derivatives with a given un- 
derlying rate tenor should be priced and hedged using vanilla interest 
rate market instruments with the same underlying rate tenor. 

The post-crisis market practice may be summarized in the following 
working procedure: 


1. Build one discounting curve C4 using the preferred selec- 
tion of vanilla interest rate market instruments and bootstrapping 
procedure; 

2. Build multiple distinct forwarding curves Cfi, ..., Cm us- 
ing the preferred selections of distinct sets of vanilla interest 
rate market instruments, each homogeneous in the underlying 
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Xibor rate tenor (typically with 1M, 3M, 6M and 12M tenors) and 
bootstrapping procedures. 


3. Compute the relevant forward rates with tenor f using the corres- 
ponding yield curve as 
pr(t, Ti-1) — pe(t, Ti) 
wl. Ti) - prt, T 
where t < T1 < T, ie{1,...,m}; 


Fy(t; Ti, Ti) = 


4. Compute cash flows c; as expectations at time t of the corres- 
ponding coupon pay-offs x;(Ff) with respect to the discounting 
T;-forward measure associated to the numeraire pg(t, Tj), as 


T; 
ci = c(t; Ti, Ti) = EY i [xF]; 
5. Compute the relevant discount factors pa(t, T;) from the discount- 


ing yield curve Cy; 


6. Compute the derivative’s price at time t as the sum of the discoun- 
ted cash flows, 


m m 
Ti 
(tT) = J peT - elt; Ti, i) = J palt, Ti) E [rE]; 
i=l i=1 
7. Compute the delta sensitivity with respect to the market pillars of 


each yield curve Cg, Cfi, . . ., Cm and hedge the resulting delta risk 
using the suggested amounts of the corresponding set of vanillas. 


The FRAs are now priced under the T forward measure associated to 
the numeraire p,(t, T2) as 


T: 
PRAM Ti, T2, K) = palt, T2) + Ta(T1, Ta) (ER ILT, T2)]- Kf , 
= px(t, T2) - t(T1, T2) [Fx(t; T1, T2)- K]. 


21.1.3.3 Pricing in Single- and Multi-Curve Framework 


Now, we will study and compare a 5.5Y maturity EUR floating swap 
leg on Euribor1M (not directly quoted in the market) in the two 
frameworks. 

In the single curve framework this is commonly priced us- 
ing discount factors and forward rates calculated on the same 
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depo-futures-swap curve cited earlier. The corresponding delta risk 
is hedged using the suggested amounts (hedge ratios) of 5Y and 6Y 
Euribor6M swaps. Notice that the expectation in step 3 is taken with 
respect to the pricing measure Z associated to the numeraire 
Pat, Ti) of the discounting curve. Any other equivalent measure 
associated to different numeraire may be used as well. 

In the multiple curve framework the forward rates are calculated 
on the Cy forwarding curve, bootstrapped using Euribor1M vanillas 
only, plus discount factors calculated on the discounting curve Ca. 
The delta sensitivity is computed by shocking one by one the market 
pillars of both Cım and C4 curves and the resulting delta risk is hedged 
using the suggested amounts (hedge ratios) of SY and 6Y Euribor1M 
swaps plus the suggested amounts of 5Y and 6Y instruments from the 
discounting curve C4. 

In the single curve framework approach, a unique yield curve is 
built and used to price and hedge any interest rate derivative on a given 
currency. This is equivalent to assuming that there exists a unique fun- 
damental underlying short rate process able to model and explain the 
whole term structure of interest rates of all tenors. It is also a relative 
pricing approach, because both the price and the hedge of a derivative 
are calculated relatively to a set of vanillas quoted in the markets. We 
notice also that it is not strictly guaranteed to be arbitrage-free, because 
discount factors and forward rates obtained from a given yield curve 
through interpolation are, in general, not necessarily consistent with 
those obtained by a no-arbitrage model; in practice bid-ask spreads and 
transaction costs hide any arbitrage possibilities. 

The multiple-curve framework is consistent with the present mar- 
ket situation, but it is also more demanding. First, the discounting 
curve clearly plays a special and fundamental role and must be built 
with particular care. This “pre-crisis” obvious step has become, in the 
present market situation, a very subtle and controversial point. In fact, 
while the forwarding curves construction is driven by the underlying 
rate homogeneity principle, for which there is (now) a general market 
consensus, there is no longer, at the moment, any general consensus 
for the discounting curve construction. 

At least two different practices can be encountered in the market: 


1. the old “pre-crisis” approach (e.g. the depo, futures/FRA and swap 
curve cited earlier) that can be justified with the principle of 
maximum liquidity, and 


21 A New Framework 545 


2. the OIS curve, based on the overnight rate (Eonia for EUR), con- 
sidered as the best proxy to a risk-free rate available on the market 
because of its 1-day tenor, justified with collateralized (riskless) 
counterparties. 


Second, building multiple curves requires multiple quotations: many 
more bootstrapping instruments must be considered (deposits, fu- 
tures, swaps, basis swaps, FRAs, etc., on different underlying rate 
tenors), which are available in the market with different degrees of 
liquidity and can display transitory inconsistencies. 

Third, non-trivial interpolation algorithms are crucial for producing 
smooth forward curves. 

Fourth, multiple bootstrapping instruments imply multiple sensitiv- 
ities, so hedging becomes more complicated. 

Last but not least, pricing libraries, platforms, reports, etc., must 
be extended, configured, tested and released to manage multiple and 
separated yield curves for forwarding and discounting, not a trivial task 
for quants, risk managers, developers and IT people. 

The static multiple-curve pricing and hedging methodology de- 
scribed earlier can be extended, in principle, by adopting multiple 
distinct models for the evolution of the underlying interest rates with 
tenors fi, ..., fn to calculate the future dynamics of the yield curves 
Cri, .-., Cj and the expected cash flows. The volatility/correlation 
dependencies carried by such models imply, in principle, bootstrap- 
ping multiple distinct variance/covariance matrices and hedging the 
corresponding sensitivities using volatility- and correlation-dependent 
vanilla market instruments. 


21.1.4 Bootstrapping with Multiple Curves 
21.1.4.1 Settings 


A yield curve is a complex object that results from many different fea- 
tures that concur to shape the curve. We have different types of yield 
curves, for example, the discount factor curve CP, the zero-coupon 
yield curve C‘, and the instantaneous forward rate curve C. 
Since the discount factor curve is observed to be monotonically 
decreasing, the zero rates are chosen to be continuous, as in 

f 


Px(to, T) = exp {—Zx(to, T) - Te(to, T)} = exp 4 - i fix(to, u)du 
0 
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The associated year fraction must be monotonically increasing with 
increasing time intervals and additive, such that 


tc(T1, T2) + tc(T2, T3) = tc(T, T3). 


The day count convention satisfying the aforementioned conditions 
that will be used is the common actual/365(fixed), such that 


tc(Ti, T2) := t[T T2; actual/365(fixed)] = (T2 — T1)/365. 


The forward rates are chosen to be simply compounded as with an 
associated year fraction (Euribor rates are quoted as actual/360 so that 


tF(Ti, T2) := t[Ti T2; actual/360(fixed)| = (T2 — T,)/360. 


The reference date, tọ can be, today, spot (two business days after 
today according to the chosen calendar) or, in principle, any business 
day after today. Once the yield curve at spot date is available, the cor- 
responding yield curve at today can be obtained using the discount 
between these two dates implied by O/N and T/N deposits. 

We choose a time grid of the yield curve as a predetermined vector 
of dates, defined by the set of maturities associated with the selec- 
ted bootstrapping instruments. The first point in the time grid is the 
reference date fo. 

Finally we have the bootstrapping instruments, quoted in the 
market, chosen as input for the bootstrapping procedure. We will use 
an algorithm that ensures exact re-pricing of the choosen input boot- 
strapping instruments. We also use an interpolation algorithm for 
calculating the yield curve outside the time grid points. Notice that 
interpolation is also used during the bootstrapping procedure. In prin- 
ciple, we can interpolate on discounts, zero rates or log discounts. We 
build curves for all needed currencies for their specific calendars 
used to determine holidays and business days. In some cases we also 
specify the bid, mid or ask price chosen for the market instruments, if 
quoted. 


21.1.4.2 Market Instruments 


In the current market situation, similar instruments may display very 
different price levels, liquidity, and may also give erratic forward rates. 
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Therefore, the first step for multiple yield curve construction is a very 
careful selection of the corresponding sets of bootstrapping instru- 
ments while they roughly cover different maturities and overlap in 
significant areas. For this reason we select those with more liquid ones 
with a tighter bid/ask spread. 

We start by examining these instruments in detail. In order to fix the 
data set once for all, we use the quotes on liquid instruments observed 
in the market. 

We begin with the interest rate deposits (depos) which are over- 
the-counter zero-coupon contracts that start at a reference date to 
(today or spot), and pay the interest accrued until maturity with a 
given rate fixed at fo. Let ROP °(to, T;) be the quoted rate associ- 
ated to the i:th deposit with maturity 7; and underlying rate tenor 
x = tọ months. The implied discount factor at time T; is given by the 
following relation: 


1 
ee a arr 
1+Ry (to, Ti): tr(to, Ti) 
The previous expression can be used to bootstrap the yield curve C, 
at point 7;. 

We continue with FRAs which are forward starting deposit contacts. 
For instance the 3x9 FRA is a 6 months deposit starting 3 months for- 
ward. In some markets FRAs are quoted between IMM days. In EUR, 
FRAs do concatenate exactly; for example, the 6x9 FRA starts when 
the preceding 3x6 FRA ends. The underlying forward rate fixes two 
working days before the forward start date. 

Market FRAs provide direct empirical evidence that a single curve 
cannot be used to estimate forward rates with different tenors. This 
can be seen if we observe the levels in the market. For instance, if 
the 1x4 FRA3M between March 18 and June 18, (TF:1x4 = 0.25556), is 
F me = 1.696% and the 4x7 FRA3M between June 18 and September18, 
(TF:4x7 = 025556), is F ie = 1.580%, we can compound these two rates 
to obtain the implied 1x7 FRA6M between March 18 to September 18, 
(TF: 1x7 = 0.50556), as 


Px(to, Ti) = to < Tj. 


k k 
po = (1 + Fi 4 TF;1x4) i (1 i F1 TF;4x1) =N = 1.641% 


while the market quote for the 1x7 FRAGM might be F m = 1.831%, 
that is, 19 basis point higher. The difference is a liquidity risk premium 
in the market. 
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Market FRAs can be used together with deposits to construction 
of the short-term structure of the yield curve Cy. If Fy(t; Tj-1, Tj) is 
the i:th Euribor forward rate resetting at time T;-ı with tenor x = 7; — 
T; months associated to the i:th FRA with maturity 7;, the implied 
discount factor at time 7; is obtained by 


Px(to, Ti-1) 
1+ Fy(to, Tii, t) © te(Ti-, Ti)’ 


Px(to, Ti) = to < Tia < Ti. 

Note that the closer the forward starting date of each FRA 7j_; gets to 
the current date fo the closer each FRA rate gets to the matching spot 
deposit rate 


lim F (to, EEP T) = RIP? (tọ, T;). 
Ti-1—> to 

Interest rate futures are the exchange-traded counterparts to the 
over-the-counter FRAs. While FRAs have the advantage of being more 
customizable, futures are highly standardized contracts which can be 
bought and sold at some exchanges. The most common contracts are 
traded at the Intenational Money Market in Chicago (so-called IMMFu- 
tures). They refer to Euribor3M or USD Libor3M. Some contracts are 
also traded in London at LIFFE. These standardized contracts expire 
the third Wednesday, every March, June, September and December 
(the IMM dates). The closing rates for these contracts will be fixed 
by the exchange on the third Wednesday of the maturity month, the 
last trading day being the preceding Monday (because of the 2 days of 
settlement). 

Notice that such a date grid is not always regular: if S; is the maturity 
date of the 7’th futures, then S; and T;, such that t(S;, Ti) = 3M, are the 
underlying FRA3M start and end dates, respectively, and, in general, 
T; # Siz. On some exchanges, futures are quoted in terms of prices 
instead of rates, for example, eurodollar futures. Then, the relation 
between rate and price is 


PM'(ty, Si, Ti) = 100 —RE(to, Si, T). 


Since exchange-traded futures have a daily marking-to-market mechan- 
ism they do not have exactly the same pay-off as FRAs. An investor 
who is long a futures contract will have a loss when the futures price 
increases (and the futures rate decreases) but he will finance such loss 
at a lower rate. On the other hand, when the futures price decreases 
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the daily profits will be reinvested at a higher rate. This means that the 
volatility of the forward rates and their correlation to the spot rates 
have to be accounted for; hence a convexity adjustment is needed to 
convert the rate RE “ implied in the futures price to its corresponding 
forward rate Fx, 


F (to, Si, Ti) = RE (to, Si, T) — Cyto, Si, Ti). 


In other words, the trivial unit discount factor implied by daily mar- 
gination according to the rules of the exchange introduces a pricing 
measure mismatch with respect to the corresponding FRA case that 
generates a volatility-correlation-dependent convexity adjustment. The 
calculation of the convexity adjustment thus requires a model for the 
evolution of the underlying interest rates. While advanced approaches 
are available in literature, a standard practitioner’s recipe is based ona 
simple short rate 1-factor Hull & White model. 

EUR futures on Euribor3M uses in their names, the letters “H”, “M”, 
“U” and “Z”, stand for March, June, September and December expiries, 
respectively. A future FUT3MH6 means a March 2016 expiry. 

Futures on an x-tenor Euribor can be used as bootstrapping instru- 
ments for the construction of short-medium term structure section of 
the yield curve. Be aware of the fact that futures contracts have ex- 
piration dates gradually shrinking to zero and generate rolling pillars 
that periodically jump and overlap the depo and FRA pillars. Hence 
some priority rule must be used in order to decide which instruments 
should be excluded from the bootstrapping procedure. 

Given the i’th futures market quote P (to, Si, Ti) with underlying 
FRA maturity 7;, the implied discount factor at T; is given by 


Px(to, Ti-1) 


Pxlto, T) = — SS 
eee 1+ {RFt(tg, Si, Ti) — Coto, Si, TD} - tr(Si, T) 


The aforementioned expression can be used to bootstrap the yield 
curve Cy at point T; once point S; is known. 

Interest rate swaps are OTC contracts in which two counter- 
parties agree to exchange fixed against floating rate cash flows. The 
EUR market quotes standard plain vanilla swaps starting at spot date 
with annual fixed leg vs. floating leg indexed to x-months Euribor rate 
paid with x-months frequency. Such swaps can be regarded as port- 
folios of FRA contracts (the first one being actually a deposit). The 
day count convention for the quoted (fair) swap rates is 30/360 (bond 
basis). 
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Swaps can be selected as bootstrapping instruments for the con- 
struction of the medium-long term structure section of the yield 
curve. 

By setting To = So = t = to and Tn = Sm = T; = S; in 


m 


D plt, S) TrCSj-1, 8) Ft, S;-1,8)) 


is 4 t, T T y t, T, 
SET _ Px(t, To) - px(t, Tn) 


Ax(t, T) Ax(t, T) 


3 


we obtain 


j 
X px(to, Sk) tF(Sk-1, Sk)Fx(to, Sk-1, Sk) 


k=l 
S t, T; = ’ 
x ? Ax(to, Tj) 
jl 
= | So prlto, Se)t(Sk-1, Se)Fx(to, Set, Sk) + Px(to, Sr-1) — Pu(to, Ti) 
k=l 
1 


x > 
Ax(to, Ti-1) + px(to, Ti) ts(Ti-1, Ti) 


where the last discount factor p,(to,7;) has been separated in the 
second line and ts5(7), T2) = t[T1, T2, 30/360(bond basis)|. A,x(t, T) 
is the annuity factor. This can be inverted to find p,(to, Ti) as 


jI 
Px(to, Ti) = bs Px(to, Sk)TH(Sk-1, Sk)F (to, Sk-1, Sk) 
k=l 


1 
+ px(to, Sk-1) — Sx(to, Ti-1)Ax(to, r] x 


1+ Sx(to, T))ts(Ti-1, Ti) 


This last formula can be used to bootstrap the yield curve at point T; = 
S; once the curve points at {7),..., Ti } and {S1, . . . , Sj-1 } are known. 
Since the fixed leg frequency is annual and the floating leg frequency is 
given by the underlying Euribor rate tenor, we have that {7,...,7;} C 
{S1,...,5; = Tj} for any given fixed leg date T;. Hence some points 
between p,(fo, Ti-1) and p,(to, Ti) may be unknown and one must 
resort to interpolation and, in general, to a numerical solution. 
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Interest rate (single-currency) basis swaps are floating vs. floating 
swaps admitting underlying rates with different tenors. The EUR mar- 
ket quotes standard plain vanilla basis swaps as portfolios of two swaps 
with the same fixed legs and floating legs paying Euribor xM and yM, 
for example, 3M vs. 6M, 1M vs. 6M, 6M vs. 12M, etc. Basis swaps are 
a fundamental element for long-term multi-curve bootstrapping, be- 
cause they allow one to imply levels for non-quoted swaps on Euribor 
1M, 3M and 12M, which can be selected as bootstrapping instruments 
for the corresponding yield curves construction. If A,6y(to, Tj) is the 
quoted basis spread for a basis swap receiving Euribor xM and paying 
Euribor6M plus spread for maturity T;, we simply have 


Sx(to, Ti) = Sem(to, Ti) + Axem (to, Ti). 


Note that the bootstrapping of yield curves requires extrapolation of 
basis swap quotations. 


21.1.4.3 Interpolation 


The choosen interpolation method determines how reasonable the 
yield curve will be. For instance, linear interpolation of discount 
factors is an obvious but extremely poor choice. Linear interpolation 
of zero rates or log discounts are popular choices leading to stable 
and fast bootstrapping procedures, but unfortunately they produce 
horrible forward curves, with a sag saw or piecewise-constant shape. 

In Fig. 21.5, we show examples of bad (but very popular!) inter- 
polation scheme, linear interpolation. The lower curve is the Swedish 
zero swap curve in April 2016 based on deposits (O/N, T/N, 1W, 1M, 
2M and 3M), IMM FRAs (with maturities between June 2016 and June 
2018) and swaps GY-LOY, 12Y, 15Y, 20Y and 30Y). The upper curve 
shows the forward curve. 

While the zero curve displays a smooth behaviour, the forward 
curve is non-smooth with oscillations that can exceed several basis 
points. Such discontinuities in the forward curves correspond to angle 
points (knees) in the zero curves, generated by linear interpolation 
that forces them to suddenly turn around a market point. 

It is recommended to only use the most liquid swaps, with ma- 
turities 3-10, 12, 15, 20, 25 and 30 years in the bootstrapping. 
The remaining less liquid quotations for 11, 13, 14, 16-19, 21-24, 
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Fig. 21.5 A bootstrap of SEK swap curve with linear interpolation. This shows the 
very bad shape of the forward curve 


26-29 years maturity should only be included in the linear interpol- 
ation schemes in order to reduce the amplitude of the forward curve 
oscillations. 

The choice of cubic interpolations is a very delicate issue. Simple 
splines suffer of well-documented problems such as spurious inflec- 
tion points, excessive convexity and lack of locality after input price 
perturbations. Some researchers found the classic Hyman monotonic 
cubic filter applied to spline interpolation of log discounts to be the 
easiest and perhaps the best approach. Its monotonicity ensures non- 
negative forward curves and actually removes most of the unpleasant 
waviness. Notice that the Hyman filter can be applied to any cubic 
interpolants. This helps to address the non-locality of spline using 
alternative more local cubic interpolations. 


21.1.4.4 No-Arbitrage and Forward Basis 


Now, we wish to understand the consequences of the assumptions 
given earlier in terms of no arbitrage. First, we notice that in the 


4 James M. Hyman. Accurate monotonicity preserving cubic interpolation. SIAM Journal on 
Scientific and Statistical Computing, 4(4):645-654, 1983. 
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multiple-curve framework the classic single-curve no-arbitrage rela- 
tions are broken. For instance, 


Patt, T2) = pa(t, T1) - pf(t, Ti, T2) fT, <7) 


1 
(t, Ti, 72) = : 
Pf 1+ Fp(t; Ti, T2) - (T1, T2) 


No arbitrage between distinct yield curves Cq and Cy can be immedi- 
ately recovered by taking into account the forward basis, the forward 
counterparty of the quoted market basis, as 


1 


O_O 
f 1+ Fa(t; T1, T2) - BA; Ti, T2) - ta(T1, T2) 


or through the equivalent simple transformation rule for forward rates 
Filt; Ti, T2) (11, T2) = Falt; Ti, T2) ta(T1, T2)BAsa(t; Ti, To). 


From this equation we can express the forward basis as a ratio 
between forward rates or, equivalently, in terms of discount factors 
from C4 and Cy curves as 


Filt; Ti, T2) (T1, T2) _ palt, T2) - prt, T1) — ppt, T2) 


BA f(t; Ti, Tz) = : = . 
Falt, Ti, T2) - ta(Tı, T2) — p(t, T2) - paft, Ti) — palt, T2) 


Obviously the following alternative additive definition is completely 
equivalent 
1 
1 + [Fa(t; T1, T2) + BA' falt; Ti, T2)| - ta(T1, T2) 
F(t; Ti, T2) « 71, T2) — Falt; T1, To) - ta 1, T2) 
Ta(Tı, T2) 

a 1 K Tı) palt, =| 

ta(T1, T2) Lpg(t,T2) palt, T2) 
= F(t, T1, T2) [BA a(t; T1, T2)— 1] 


p(t, Ti, Tz) = 


BA' a(t, T1, T2) = 


which is more useful for comparisons with the market basis 
spreads. Notice that if Ca = Cy we recover the single-curve case 
BAya(t, Ti, Tz) = 1, BA A(t, Tı, Tz) = 0. 
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The forward basis in the previous equations is a straightforward con- 
sequence of the aforementioned assumptions, essentially the existence 
of two yield curves and no arbitrage. Its advantage is that it allows for 
a direct computation of the forward basis between forward rates for 
any time interval [7;, T2], which is the relevant quantity for pricing 
and hedging interest rate derivatives. In practice its value depends on 
the market basis spread between the quotations of the two sets of 
vanilla instruments used in the bootstrapping of the two curves C4 
and Cy. 

The approach can be inverted to bootstrap a new yield curve from 
a given yield curve plus a given forward basis, using the following 
recursive relations: 


Pa Die” Baii Pd ate Pd 
i= i1 = Jai 
' Pf. i1 — Pfi + pfi- BAfa,i Pf i-1— pfi- BA’ fai + Tai i 
Pd,i Pd,i 
Pfi Pf,i-1 = Pf,i-1, 


~ Pait (Pais — Pai) BAfa,i -Pijit Pori + BA fai + tai 


where tx(Ti-1, Ti) = Tx, Px(t, Ti) = px i and BAfg(t, Ti1, Ti) = BAfa,i. 
Given the yield curve C, up to step px,i-1 plus the forward basis for the 
step i— 1 — i, the earlier equations can be used to obtain the next 
step Px,i- 

We now discuss a numerical example of the forward basis in a 
realistic market situation where we consider the four underlying 
interest rates J = {Lhm, Bm, 16m, lzm}, where I = Euribor index, 
and we bootstrap from market data five distinct yield curves C = 
{Ca, Cim, C3u, Cou, Ci2m}, using the first one for discounting and 
the others for forwarding. 

We build the discounting curve Cg following a “pre-crisis” tra- 
ditional recipe from the most liquid deposit, IMM futures/FRA on 
Euribor3M and swaps on Euribor6M. The other four forwarding curves 
are built from convenient selections of deposits, FRAs, futures, swaps 
and basis swaps with homogeneous underlying rate tenors; a smooth 
and robust algorithm (monotonic cubic spline on log discounts) is 
used for interpolations. Different choices (e.g. an Eonia discount- 
ing curve) as well as other technicalities of the bootstrapping are 
described by Ferdinando M. Ametrano and Marco Bianchett (see 
references). 
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The corresponding multiplicative forward basis curves can be 
calculated as 
F(t; Ti, T2): (T1, T2) _ palt, T2) - ppt, Ti) — ppt, T2) 


BAga(t; Ti, Tz) = = 
Falt; Ti, T2) - ta(Tı, T2) — pe(t, T2) - palt, T1) — palt, T2) 


and the additive forward basis are given by 


F(t; Ti, T2) - (11, T2) — Fa(t; T1, T2) + ta(T1, T2) 
Ta(Tı, T2) ' 


Eee: ee Tı) palt, ro 
Ta(Tı, T2) | pft, T2) palt, T2) 


= Falt, T1, T2) [BA y(t; T1, Tz) - 1]. 


BA' a(t, Ti, T2) = 


21.1.5 Modern Pricing 


According to Bianchetti and Morini (2010), the new market situation 
has induced a sort of “segmentation” of the interest rate market into 
sub-areas, mainly corresponding to instruments with 1, 3, 6 and 12 
months underlying rate tenors. These are characterized, in principle, 
by different internal dynamics, liquidity and credit risk premia, reflect- 
ing the different forwarding horizon views and interests of the market 
participants (“on average” i.e. views are weighted by how much Cap- 
ital each investor is willing and able to invest according to his/her 
opinion). 

In response to the crisis, the classical pricing framework, based 
on a single yield curve used to calculate forward rates and discount 
factors, has been abandoned, and a new modern pricing approach has 
been created by practitioners. This new methodology takes into ac- 
count the market segmentation as an empirical fact and incorporates 
the new interest rate dynamics into a multiple curve framework as 
follows. 


e Discounting curves: these are the yield curves used to discount 
futures cash flows. These curves must be constructed and selected 
so that they reflect the funding cost of the bank in combination 
with the actual nature of the specific contract that generates the 
cash flows. In particular: 
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o an OJS-based curve is used to discount cash flows generated 
by a contract under CSA with daily margination and overnight 
collateral rate; 


o a funding curve is used in case of contracts without CSA; 


o incase of non-standard CSA (e.g. different margination frequency, 
rate, threshold, etc.), an appropriate funding curve should in 
principle be selected, but we will not discuss this topic here. 


e We stress that the funding curve for non-CSA contracts is specific 
to each counterparty that will have its specific funding curve. This 
modern discounting methodology is called CSA-discounting. 


e Forwarding curves: these are the yield curves used to compute 
forward rates. As discussed before, the curve must be construc- 
ted and selected according to the tenor and typology of the rate 
underlying the actual contract to be priced. For instance, a swap 
floating leg indexed to Euribor6M requires a Euribor6M forwarding 
curve constructed from quoted instruments with Euribor6M as the 
underlying rate. 


During the pricing process it is important to generate all cash flows 
from the curve based on instruments with the same tenor and then 
discount with the matching discount curve. It is also important to use 
an exact fit bootstrap method and a good interpolation method. 


21.1.6 Pricing Under Collateralization 


The financial crisis in 2008 was a catalyst for significant changes in 
the market. Financial practitioners witnessed a tremendous increase 
in basis swap spreads, implying a divergence from implied rates and 
traded rates in interest rate markets. When collateral agreements in- 
creased in use and when counterparties began discounting at the 
overnight rates dictated by the credit support annex, the world had 
forever changed. 

As a result, many financial institutions are currently in the process 
of migrating to new market standards. But questions remain as to the 
potential impact on existing portfolios and how to effectively manage 
instruments with longer-dated maturities when spreads in Libor vs. 
OIS rates diverge. 

Collateral discounting and the impact of standardization in the mar- 
ket are adding a whole new level of complexity to derivative pricing 
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and risk management. Market participants are seeking a deeper under- 
standing when it comes to the potential consequences of moving to 
collateral discounting. 

With such collateral agreements, also called CSA agreements, we 
have to make changes in the valuation formulas. The swap market has 
already moved to the dual-curve approach, with the London Clear- 
ing House (LCH) using OIS discounting for clearing swaps, and the 
International swaps and Derivatives Association (ISDA) standardiza- 
tion on the imminent horizon. Other markets, including swaptions, 
caps/floors, exotics and equities, are still evolving. Now, receivers of 
fixed rates gain under OIS discounting, while payers lose when com- 
pared to the old ways. Taking into account recent studies of the swap 
market moving towards the multi-curve approach, we come to the 
conclusion that swaps can be significantly mispriced under the single 
curve framework. 

In the market, one can observe that spreads and OIS rate risks are 
close for par swaps, but differ significantly for out-of-the-money (OTM) 
swaps (vs. ATM swaps). Under the single curve approach, one neglects 
the risk that actually exists, while under the dual-curve approach one 
can estimate these risks. 

When we price under a collateral agreement we make two simplify- 
ing assumptions: 


1. We consider full collateralization (with zero thresholds) by cash. 


2. The collateral is adjusted continuously with zero minimum transfer 
amount (MTA). 


Most CSA contracts include a threshold and an MTA to avoid posting 
small amounts of cash each day. The CSA contracts can be one sided or 
bilateral. About 16% of the agreements are unilateral (one sided) while 
84% bilateral. 


21.1.6.1 Problems in the Old-Style Implementation 


We start by considering a plain vanilla tenor basis swap,’ typic- 
ally 3-months floating against 6-months floating rate as illustrated in 
Fig. 21.6. 


> Tt is also common that payment of short-tenor leg is compounded and paid at the same time 
with the other leg. 
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Libor (short tenor) 
+ spread 


Libor (long tenor) 


Fig. 21.6 A 3-month floating rate (the upper cash flows) against a 6-month float- 
ing rate (the lower cash flows). The arrow above the upper “wave” represents the 
spread over the floating rate 


In the old-style implementation, the spread between the rates in 
Fig. 21.6 was typically zero (or very near zero). However, since 2008 
the spread is quite significant and also volatile as seen in Fig. 21.7, 
where tenor swap spreads for maturities 1, 3, 5, 7, 10 and 20 years are 
shown. This spread represents the difference in risk; the long tenor 
has a higher risk than the short tenor. 


50 
45 — USD BASIS SWAP 3v6 1 YR 
——USD BASIS SWAP 3v6 3 YR 
40 ——USD BASIS SWAP 3v6 5 YR 
35 — USD BASIS SWAP 3v6 7 YR 
i — USD BASIS SWAP 3v6 10YR 
30 ——USD BASIS SWAP 3v6 20 YR 
25 F 
7 The upper curve is the 1 YR 
20 and lowest the 20 YR 
15 
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5 
0 = 
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Fig. 21.7 Historical data for USD 3-month vs. 6-month TS spread. The curves are 
given in the same order as the legends® 


6 Source Bloomberg. 
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USD Notional 
JPY Notional 


JPY Libor — 3 month + spread 


USD Libor — 3 month 
JPY Notional 


Fig. 21.8 A 3-month floating rate in JPY (with a constant spread) against a 3- 
month floating rate in USD. The arrow above the upper “wave” represents the 
JPY spread over USD 


USD Notional 


ý 


Similar problems are found if we consider cross-currency swaps 
(CCS or CIRS) as shown in Fig. 21.8. 

In the old-style implementation the spread in a CCS was zero. How- 
ever, since 2007 the spread is quite significant and volatile. Here we 
have seen a drastic change in recent years. This is shown in Fig. 21.9 
where we show the cross-currency USD/JPY spreads for maturities 1, 
3,5, 7 and 10 years. 

Before 2008, most banks did their funding on the interbank mar- 
kets; that is, the trades were made without collateral. Then we had the 
situation as in Fig. 21.10. 

According to the old view, the purchase of an OTC derivative was 
funded with an unsecured external loan. 

Here counterparty R (Red) buys a derivative of B (Blue) where the 
value of the derivative is given by some optional payment. To fund the 
trade, R takes a loan from another party or internally in his own bank.’ 
The interest A pays on this loan is (Supposed to be) given by Libor. We 
then have unsecured funding. A Libor rate is §, an unsecured offer rate 
in the interbank (deposit) market and Libor discounting is appropriate 
for such unsecured trades between financial firms with Libor credit 


7 If the trader is a person who funds the trade with his/her savings from his/her own money- 
market account, he/she will typically earn less interest on this amount than the running rate on 
his/her own money-market account. 
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The upper curve is the 10 YR 
and lowest the 1 YR 


2007-01-01 2007-07-01 2008-01-01 2008-07-01 2009-01-01 2009-07-01 2009-01-01 


Fig. 21.9 Historical data for USD/JPY cross-currency spread. The curves are given 
in the reverse order as the legends® 
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LIBOR 
Funding 


Cash 
Cash = PV 


Fig. 21.10 Funding via the nterbank market 


quality. In addition, Libor discounting makes the present value of the 
loan equal to zero. 

After the credit crisis in 2008, banks started requesting collateral. 
This modifies the pricing of OTC derivatives. The current view is to 
require collateral, so the funding is secured, and the picture looks like 
Fig. 21.11. 

In the new view, we finance the trade by “a loan” from the 
counterparty via collateral from counterparty B. 

Here, the outright cash flow, that is, the collateral is equal to the 
PV of the trade so no external funding is needed. And when the PV of 
the trade changes, R pays back or receives more collateral from B. This 
collateral is Gdeally) paid/received once every day. Also, the party who 


8 Source Bloomberg. 
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Option Payment 


Funding 


Collatral Rate 


Cash = PV 
Fig. 21.11 Funding via collateral 


holds the collateral must pay the other party interest. The size of this 
interest is specified in the collateral agreement. Since the collateral is 
posted on a daily basis, the relevant default risk is an overnight risk. 
Therefore, the collateral rate is usually the over-night (O/N) rate, so we 
use OIS to bootstrap the discount curve. When funding is provided by 
the collateral agreement Libor discounting is inappropriate. 


21.1.6.2 Collateral Agreements 


In Table 21.1 we provide some details for a given collateral agreement. 

In many CSA agreements you also have the opportunity to post (gov- 
ernment) bonds. When posting bonds there normally is a haircut. This 
haircut is a percentage that is subtracted from the market value of 
the asset that is being used as collateral. The size of the haircut re- 
flects the perceived risk associated with holding the asset. The higher 
the haircut, the safer the loan as regarded by the lender. For example, 


Table 21.1 A simplified collateral agreement 


Base currency USD 

Eligible currency USD, EUR, GBP 
Independent amount 5M 

Haircuts [Schedule] 

Threshold 50M 

Minimum transfer amount 500,000 

Rounding Nearest 100,000 USD 
Valuation agent Party A 

Valuation date Daily, NewYorkBusinessDay 
Notification time 2PM, NewYorkBusinessDay 
Interest rate OIS, EONIA, SONIA 


Daycount Act/360 
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United States Treasury bills, which are seen as fairly safe, might have a 
haircut of only 10%. For most other bonds the collateral is higher. 

In most CSA agreements you also have to consider the MTA and the 
threshold amount. The MTA is the smallest amount of currency that 
is allowed for transfer as collateral. For large banks, the MTA is usually 
in the USD 100,000 range, but can be lower. The threshold amount is 
the amount of unsecured credit risk that two counterparties are willing 
to accept before any demand for additional collateral demand will be 
made. The counterparties typically agree to a threshold amount prior 
to committing to the deal, but this is a source of on-going friction 
between OTC counterparties and their brokers. 

All these choices in the CSA agreements complicate the application 
of the theoretical priniciples. However, we will omit this here since 
it is a rather complex subject. But we need to remember this when 
dealing with the full valuation in a practical situation. For simplicity, 
we will pursue the analysis of CSAs with MTA and threshold amount 
set to zero and in addition assume that the posting frequency is 
instantaneous. Instantaneous posting enables us to approximate the 
collateral with an integral instead of a sum. Furthermore, we only use 
cash as collateral. 

Also note that, except for the MTA and thresholds, CSA typically 
does not involve hedging 100% of the credit exposure. As illustrated in 
Fig. 21.12, a 6-year swap (paying fixed at 2.5% on $1,000,000 notional) 


— Credit exposure ==- Swap value 


Fig. 21.12 A 6-year swap paying fixed rate at 2.5% on 1,000,000 notional. The 
collateral amount is the difference between the credit exposure and the swap 
value 
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cannot be 100% funded with collateral. The reason is that the value of 
the trade changes faster than the collateral frequency, and MTA and 
threshold considerations must also be taken into account, in addition 
to cash flow events. We must also consider the movements in the value 
of the collateral itself. 

Fig. 21.12 shows the total exposure of a swap over the lifetime of 
the trade, and a hypothetical collateral amount. Since imperfections 
exist, the total exposure of the trade is not perfectly hedged. The mis- 
match between the current exposure and the amount in the collateral 
account means funding is not purely OIS. This is the reason for funding 
value adjustments (FVA), which we will discuss later. 


21.1.6.3 Overnight Index Swaps 


For an OIS the floating rate is a daily compounded O/N rate and the 
market quotes; the fixed rates are called OIS rates” (see Fig. 21.13). 

Figure 21.13 shows the Libor vs. OIS spread in USD and JPY over a 
few years. According to the old view, an inappropriate kind of dis- 
counting was used for secured trades which mispriced future cash 
flows. This has as a significant impact on multi-currency trades where 
the change is given by 


Change ~ Notional x Duration x (Difference in discount rate) 


This is also inconsistent with CCS. The implied Foreign Exchange (FX) 
forwards are off market, leading to mispricing of foreign Libors. Long- 
dated FX products are most affected. We also get wrong forward Libors 
and overestimation of forward Libors with short tenors. 

Fig. 21.14 shows the 3-month Libor vs. OIS spread in USD and JPY 
over a few years. 


OIS rate 


Compounded O/N 


Fig. 21.13 A typical overnight index swap 


? Usually, there is only one payment for tenors (maturities) shorter than a year. 
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Fig. 21.14 A 3-month Libor rate vs. the overnight index swap spread in USD 
and JPY 


The bootstrapping of OIS curves is made by the same formula that 
we use when we bootstrap interest rate swaps 
D(t, To) -RGS Dpi AnDUt Tn) 


PE T+ROIS. A, 


Here we get the discount factors from the quoted OIS rates ROS with 
maturity Ty and day fraction A,„. The continuously compounded zero- 
coupon rate is then calculated by 


In (DM, Ty) 


dy 
365 


Z9" = -100 


21.1.7 Pricing with Collateral Agreements 


With no collateral (or collateral in domestic rate rf), the price of a 
contingent claim paying II(T) at maturity T is given by the Feynmann- 
Kaè representation 


T 
BA) 
diy — ro d md — FQ . y4 
11? (t) = E2! | exp J r“(u)du ¢ - OTF, | =E É XD I 7|, 


t 
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Here the expectation is taken with respect to the domestic martingale 
(risk-neutral) probability measure Q1 and the domestic money-market 
account 


t 
B“(t) = exp I r{(u)du 
0 


dB4(t) = r4(t)B4(t)dt 
B4(0) = 1 
SO 


T 


exp - f udu 


t 


_ BO 
“BY 


If there is a collateral agreement where we are supposed to get or pay 
cash collateral in a foreign currency f, then we have to add an extra 
term in the pricing formula for the aforementioned contingent claim. 
Say that we are in the money and get collateral from our counterparty. 
In this case we are free to invest this collateral in the market. However, 
we also have to pay interest to our counterparty since he/she still owns 
the collateral. This collateral rate is specified in the CSA agreement 
between the counterparties. If we denote the instantaneous return 
(or cost when it is negative) by holding the collateral in the foreign 
currency f at time t by y(t), we have 


YFA =r O-O. 


Here 7f (t) denotes the “risk-free” interest rate, assumed to be equal 
to the funding or repo rate for uncollateralized assets and cf (t) the 
collateral rate in the foreign currency f. In the case of cash collateral, 
the collateral rate c(t) is usually given by the overnight rate, that is, the 
OIS rate of the corresponding currency. When using treasury bonds 
as collateral, the collateral rate is usually the repo rate. For corporate 
bonds, we can use Libor plus a spread. 

The collateral rate is pre-specified in the collateral agreement. There- 
fore, we let r(t) accumulate in a money-market (bank) account and c(t) 
in a collateral account. When we get collateral, we invest this money 
in the market. 

The spot measure is taken to be that measure under which uncollat- 
eralized assets grow at r(t), and collateralized assets grow at c(t). Both 
uncollateralized p(t, T) and collateralized D(t, T) zero-coupon bonds 
pay 1 unit in the domestic currency at time T. 
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If we denote the PV of the derivative at time t by T14(t) Gn terms of 
the domestic currency d), the collateral amount posted from the coun- 
terparty in the foreign currency f is given by W(t) = Of (t), 
where fx% (t) is the foreign exchange rate at time t,representing the 
price of the unit amount of currency f in terms of currency d. These 
considerations lead to the following calculations for the collateralized 
derivative price: 


f 
mosen g + fet . EY jz (5 aie 
BÍ 


T 


We start by calculating the expectation 


T . 
f f djp at d 
= E24 ja, B, . Bs . B; [fx i y (4 i. 5) ds|F; 
af ? 

J By Bfe By fe 


where we have used B(s) = BU sf? (s)V; in order to change the 
measure from Qp to Q4. This can be simplified into 


f d d df d 
, rp II f BS II 
B Txs A Bi X; fxs 


0 
T 
1 Bd 
= —_ faa £. yf. fds |, 
pt" | ae ys- Tids 
0 
i T sS 
sT f - f iau -y - Mds]. 
Xt 


Using this formula we have 
T T 
11? = E" | exp - fru mip f exp - f iau -yf - Mds 
t t 


t 
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which is equivalent to 
t 


TI? = exp J r4du 
0 


T T 
x £2 | exp [die nt fo - | ridu yf das | 


mg mi , mni , 
= BY. | E% a f oli - f Zi d 
T S 0 S 


Now we see that 
t 


T 
ns T? r? r2 
poi +f Ja =e f T -vlas= me 
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since the expectation is a martingale under the domestic money- 
market account. Here we have also defined M(t). Remark that 
Ee [M(T)] = M(t). Next, use Ito’s lemma on 

$ 


1 1 
M = — -Tf + yt . T14ds. 


Bi BE 
0 
The result is 
aM aM aM 
dM, = —‘dt + — dB? + t ante 
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= exp - f ridu -ri . Mat + exp - f ridu all? 
0 0 
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or 
t 
exp J ridu $ dM; = (y; =r) fdr + arif 
0 


that is, 
t 
arig = (ot - 1!) . Idt + exp [au dM;. 
0 


We can now solve the previously mentioned Stochastic Differential 
Equation (SDE): 


dII = « - IIdt + dB 
T(t) = 2 
where we use 
aff d 
a= (0-0) 


and 


t 


dB = exp J rdu $ dM(t). 


0 
Set X = e ®™TI and use Ito 
ox ox 
dX = a tt za =œ -X -dt+a -X -dt+e®™dB = e™dB. 


If we integrate this we get 


A T 
X(T) —X(t) = [ets => eT) —e*' M(t) = |w 
t 


t 


and by taking expectation value we get 
T(t) = e" eA TT 


Finally we have the following theorem. 
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Theorem 21.1. Suppose that T14(T) is a derivative’s pay-off at time 
T in terms of the domestic currency d and that the foreign currency 
f is used as the collateral for the contract. Then, the value of the 
derivative at time t, 114(t) is given by 


T T 
MYA = E24 | exp } - J r"(s)ds $ - exp / yw (s)ds $ OKT) |, 


t t 
T 


= D4(t, T)E;" | exp J CORO) ds  T1“(T) 


t 


Here we also have defined the collateralized zero-coupon bond of 
currency d as 


T 


PGT) =E% exp {— J c4(s)ds 


t 


i oT , 


and the collateralized forward measure Tz, where the collateralized 
zero-coupon bond is used as the numeraire. 

In the case where the deal and collateral currencies are different, (d) 
and (f) respectively, we define the foreign collateralized zero-coupon 
bond D(d, f) by 


T T 
D (t, T) = EM exp 4 — | c4(s)ds -exp 4— J y (ods 


f t 


In particular, if c and yO = 4A - c4(t)) - A - cf) are 
independent, we have 


T 
ptt (t, 7) = D(t, T) - exp {— J yif (s)ds > , 


t 
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where 


S 


l 3 , 
y(t, s) = = In Ee exp 4 — ] y (sds 
s 


t 
We also have some corollaries to the Theorem 21.1. 


Corollary 21.1 When the foreign “risk-free” rate is the foreign col- 
lateral rate, that is, when rf (t) = d (t) and the collateral is posted in 
foreign currency, the PV of the derivative at time t becomes 


T 
M(t) = EP | exp 4- J r(s)ds } - TCT) | = p(t, DE (TY), 


f 


where p(t, T) is an uncollateralized zero-coupon bond and El 1 the 
expectation under the domestic forward measure with information 
given up to time t. 


Corollary 21.2 When the domestic “risk-free” rate is the domestic 
collateral rate, that is, when ri (t) = A(t) and the collateral is posted 
in domestic currency, again we have 


T 
T(t) = EÈ | exp } - J r“(s)ds $ - TI(T) | = p(t, TEs“ (TI(T)]. 


t 


Corollary 21.3 Ifthe foreign and the domestic economies are inter- 
changed so that the domestic risk-free rate is equal to the domestic 
collateral rate, that is, r(t) = c4(t), and the collateral is posted in 
domestic currency, then the previous equation is transformed into 


WO = pane! [ma], 


This result clearly shows the fact that the effective funding cost is 
given by the collateral rate, regardless of the risk-free rate of the 
corresponding currency. 
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It is also possible to construct a maturity 7 dependent measure 
associated with the collateralized zero-coupon bonds 


T 
D(t, T) = E24 | exp 4 - J cds 


t 
that is similar in some ways to forward measures. Define or (with 
expectation ET) by Q' = Z(T) - Q, where 
T 
exp }— f cHs)ds 
0 


a= DO, T) 


Then 


T 
Ez! ex l- f esas xn] 
0 


T 
Ee es l- f cous 
0 


T 
Ee |e l- f nel xn] 


DE, T) 


ETIX(T)] = 


If the spread y4(t) is deterministic, then Q’becomes the standard 
T-forward measure Q7 because 


0 1 


J ~ BUT) pXO0, TY 


T 
exp l- f [re (s) — y(s)] al 


Z(T) = 


T 
EL! |e |-/ [r4(s) — y%(s)| ds 


0 


which is the Radon-Nikodym derivative for QO" In particular, QT =Q" 
when yA) = 0 and the collateral rate is equal to the overnight rate, 
that is, c4(t) = r“(t). 
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From the definition of D(t, T), foranyO<s<t<T 


f T 
EQ exp 4— J c{(s)ds D(t, T) = EQd Ez exp 4 — | c4{(s)ds 
0 0 
T S 
a d _ d 
= E24 | exp {= | c“(s)ds > | = exp į- J c“(s)ds } D(s, T) 
0 0 


giving a result that will prove to be useful, for example, in modelling 
c(t) and D(t, T) as SDEs. 


Corollary 21.4 The value of the collateralized zero-coupon bond 
discounted by the collateral rate 


t 


exp 4 — J c{(s)ds $ Dit, T) 
0 


is a O¢-martingale, confirming thet the drift under On of D(t, T) 
re nd 
is c(t). 


In the literature OT is often referred to as the T-forward measure in- 
duced by D(t, T) as numeraire because it makes collateralized trades 
QO’ -martingales; that is, changing measures from of to oF allows 
Theorem 21.1 to be restated as the followng corollary. 


Corollary 21.5 When payment and pricing currencies are different 


T 
E? | exp -| o9- rola - TI(T) 


0 
T(t) - exp l-i [y4(s) -y (8)] ie 
0 


DET) 


and when payment and pricing currencies are the same 


T(t) 
D(t,T) 


E?” [I(T)] = 
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21.1.7.1 Pricing with Multiple Currency Collateralization 


With full collateralization in multiple currencies, the funding spread 
between currencies i and k is given by y = yi—yk = (ri —c!)—-(r* — c*). 
The pricing formula can then be written as (remember: y (t) = r(t) — 


c(t) 


T 
T(t) = E& | exp }- I max (Mw - viw) du} - 04T) |, 
jec 
t 


T 
= D4(t,T) - exp {— J max (yw) du \ . ET [nn]. 
0 


Here we use a blended cheapest-to-deliver (CTD) curve through the 
life of the trade constructed for each unique CSA agreement. If we 
have an agreement with six possible collateral currencies (USD, EUR, 
GBP, JPY, CHF and CAD), we need to bootstrap 29 different interest 
rate curves, as shown in Fig. 21.15. 

Also, note that we have the following spread relations: 


ye) = ye) Vt>0 
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Fig. 21.15 The complicated bootstrap process if 6 collateral currencies are used 
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21.1.7.2 What Currency Gives the CTD? 


The CTD currency is given by the one with the highest interest rate. 
The reason is that the standard contracts in the market always stip- 
ulate that it is the party who posts the collateral that can decide on 
the currency. By choosing the highest paying currency for the collat- 
eral he/she will get the highest rate of return. This contractual feature 
is closely related to the fact that the party who posts the collateral 
can at any time decide to get the collateral back and then post an- 
other collateral in another currency. The receiver of the collateral, on 
the other hand, has to pay interest in the currency chosen by his/her 
counterparty. 

In some agreements, the counterparties can also post bonds. When 
bonds are posted, a haircut is made depending on the credit quality of 
the bonds. 

The party holding the collateral can do what he/she likes with the 
cash or the bonds, like selling it in the market. But note that he/she 
must return exactly the same collateral when the other party wants it 
back and it can be difficult to buy back the same bond when market 
liquidity is low. 

When you have to post collateral yourself, you might not be able 
to deliver the CTD currency. Then you might have to deliver the next 
CTD currency or even another one. 

If you could decide to return the interest rate in any other currency 
than the one your counterparty chooses, then for you the currency 
with the lowest rate would be CTD. 

In Fig. 21.16, we illustrate the CTD curve for two different curren- 
cies, GBP (SONIA) and EUR (EONIA). 


21.1.7.3 Bootstrap Technique 


Before we end this section, we will describe how to bootstrap in some 
different currencies. The exact bootstrap procedure depends on what 
liquid instruments are available in each market. 

In general, curve stripping is based on market instruments such as 


e OIS 
e swaps (floating Xibor vs. fixed) 


e Basis swaps also known as tenor basis swaps (e.g. Xibor 3M vs. 
Xibor 6M) 
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Fig. 21.16 A CTD curve for two currencies, GBP (SONIA) and EUR (EONIA) 


e Cross-currency basis swaps (CCBS) (e.g. USD Libor 3M vs. GBP 
Libor 3M) 


The building blocks of these instruments are 

e Fixed cash flows: Fixed * YF * Notional 

e Libor payments: Libor * YF * Notional 

To price these we only need the elementary bits 

e Discount factors: DF = PV(1 unit of currency) 

e Forwards: FWD = PV(Libor)/DF 

Then, the PVs are given by 

e PV (Fixed cash flow) = Fixed * YF * Notional * DF 

e PVC€loating cash flow) = YF * Notional * FWD * DF, 


where YF is the year fraction. 
Once we know the discount factors for all maturities and the for- 
ward rates for all maturities and tenors that is, the dscount curve 
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t > DF(t) and the forward curves t — FWD*(t) for all tenors 5 we 
are able to price all linear instruments. 

In practice, the discount factors and the forward rates are stripped 
from market instruments for a set of maturities and tenors. For other 
maturities and tenors the values are obtained by interpolation. Here 
we have some practical issues: 

How is this interpolation performed? 

How are the curves represented? In terms of discount factors, 
forwards directly, etc.? 

Liquid OIS markets exist for the G5 currencies (USD, EUR, GBP, 
JPY and CHF) and AUD and CAD among others. If there is no OIS 
market or the market is not liquid enough we cannot find an OIS 
curve. Therefore we cannot strip the projection curves in the usual 
manner. One possible idea is then to turn to a cross-currency market 
which is liquid enough and try to simultaneously strip the pro- 
jection curve and the implied discounting curve from local swaps 
and CCBS. 

In the single-currency case, we first do the stripping of the OIS 
curve, typically done using OIS. OIS is a fixed/float interest rate swap 
with the floating leg based on the published overnight rate index. 


O/N rate + OIS swaps —> OIS curve used for discounting 


Next, we do the stripping of the projection curves (e.g. 1M, 3M, 6M, 
1Y curves) given the OIS curve. Here we use instruments indexed on 
Xibor: 


Cash deposits + FRA/futures + swapsprojection curve 


In the cross-currency case we first assume that the domestic and for- 
eign curves for all needed tenors have been already stripped. Then we 
strip the implied foreign basis curve 


FX forwards + CCBS>Foreign discount curve 


For the 3-month tenor, we now have Table 21.2. 

In the most liquid markets we have the following available curves 
(Table 21.3). 

For different currencies, the trading conventions and instrument 
dependencies are slightly different. 
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Table 21.2 The 3-month tenor in the bootstrap process 


Domestic floating leg Foreign floating leg 
3MDOM index 3MFOR index + spread 
Projection Curve DOM3M swap Curve Projection Curve FOR3M swap curve 
Discount curve DOM OIS curve Discount curve Implied DOM3M/FOR3M 
basis curve 


Table 21.3 Available curves for the most liquid markets 


Currency Overnight rate Standard curve Forward Basis curves 
curves 
USD Fed funds 3M USD Libor Muni swaps 1M vs. 3M 
effective rate 3M vs. 6M 
3M vs. 12M 
3M Prime/Libor 
Basis swap 


Muni vs 3M Libor 
T-bill vs 3M Libor 


EUR EONIA 6M Euribor 1M Euribor 3M vs. 6M 
3M Euribor 6M vs. 12M 
12M Euribor 

JPY MUTAN 6M JPY Libor 1M JPY Libour 1M vs. 3M 
3M JPY Libour 3M vs. 6M 

GBP SONIA 6M GBP Libor 1M GBP Libour 3Mvs. 1M 


3M GBP Libour 12M vs. 6M 
12M Euribor 12M GBP Libor 


CHF TOIS 6M CHF Libor 3M CHF Libour 12M vs. 6M 
1M CHF Libour 
CAD Bnak of Canada 6M CAD-BA 6M vs. 3M&3M 
Overnight Repo vs. 1M 
Rate 
(CORRA) 


21.1.8 Market Instruments 


We now start to present some formulas that are used for the curve 
generation. By a time schedule, we mean times t < Tg < T; <... < 
Ty and the corresponding year fractions { A^, PG where A,, is the year 


fraction from T,_; to Ta. We also denote by L,(T;_7) the Libor rate with 
tenor x that starts at T„—7 (the rate is usually fixed 2 days before T,_7). 


21.1.8.1 Overnight Index Swap 


If the collateral and the pay currency are the same, we have the fol- 
lowing definition of the forward overnight rate: For a given currency 
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i, collateral rate c; and times Tz and T2, the time t forward overnight 
rate is defined by 


O(t, Ti, T2) = 


1 e E | 
AT), T2) (DCG, T2) l 


An OIS in a currency i with collateral rate c; is a contract where you 
exchange a daily compounded floating overnight rate for a fixed rate 
at pre-specified points in time. The fixed rate OJSy is chosen so that 
the PV of the swap, according to market discounting conventions, 
becomes zero today 


Sm 


N M 
OlSy Y AnD(t, Tn) = Y E; | ex | = l doas 
(0) 


n=1 n=1 


Jin 
: (TI (1 + OW i m Uj-1,m, Uim)) . dim E ')}. 


jel 


where {7,} and {Sm} are the time schedules for the fixed leg and 
the floating leg respectively, {Uj} is the time schedule for the daily 
compounded overnight rate, U*;_7m is the fixing date for the rate in 
the period (Uj_7m, Ujm] and An, ôjm are the year fractions for the 
periods (Tn-1, Tn] and (Uj_1,m; Ujml. Here, Uo = Sm-1, UJm,m = 
Sm, To = So and Ty = Sy. By using the forward overnight rate and 
the approximation 


ok NN 
OF m Uj-iym Ujm) mo O(Uj-1,m: Uj-iym Ui); 
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In the last step, we were splitting up the integral and used the tower 
property of conditional expectation: if H C G then E[E[X|G]|H] 
= E[X|H] = E[E[X|H]|G]. 

Finally, we end up with the following approximation (where we 
ignore that the fixing is made in two trading days before the start of 
the time interval). 


N 
(OIS) — OlSy È AnD, Tn) = DO, To) — DO, Ty) 
n=1 
We can also get the previous equation from a simplification of a 
collateralized OIS 


N Th 
OISy Ý AnE® | exp 4 - J c(s)ds 
n=l 0 
N Th Tn 
DD exp 4 — l c(s)ds exp 4 — J c(s)ds è — 1 
n=l 0 Tn-1 


by using the collateralized zero-coupon bond. The collateralized zero- 
coupon bond price can then be bootstrapped as 


DU, To) = OS Y 7 AnD Th) 
1+ OISn g AN ` 


D(t, Ty) = 


If the collateral and the pay currency are not the same, we use that 
fair value of a payer OIS swap in currency i with collateral in currency 
j, collateral rate cj and quote OIS 


M Sm 
PV= Y bn + D(t, Sm) + exp } — / y¥(s)ds $ - O(t, Sm, Sm) 


m=1 t 


N T 
~ OISy > An - D(t, Ty - exp } = J yll(s)ds 


n=1 t 
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21.1.8.2 Forward Rate Agreements (FRA) 


Given times T; and T> the T; x T2 FRA traded in the market is a contract 
paying out to the buyer at time T; the amount 


A(T), T2) - (Lx(T1)- K) 
1- A(T), T2)-L(T1) ? 


where K is the pre-specified fixed rate of the FRA, x is the tenor starting 
at T; and ending in T2, and A(77, T2) is the year fraction for the tenor 
interval (T7, T2]. Under certain model assumptions and reasonable 
market conditions it can be shown that 


€ 


FRA) Kw Ep ? LEAT 


where the expectation is taken with respect to the forward measure at 
Tə discounted with the current (time 4 price of a zero-coupon bond 
with expiry T2 using the collateral rate in currency i (domestic). 


21.1.8.3 Interest Rate Swaps 


An IRS in currency i with collateral in currency j with collateral rate c; 
and quote IRS is given by 


N Tp. 
IRSy X} An D(t, Tn) - exp l- f sous] 
t 


n=1 


(IRS) 


M se Sm e 
= > bmD(t, Tm) i E; m [Ly (Sm-1)] exp l- f oas 
t 


m=1 


and a collateralized tenor swap (TS) by 


N . Tn |. 
D AnD(t, Tn) (E™ [LTn-1, Tn, t5)] + TSN} + exp l- (i voal 
A 
(TS) 
M ó Sm ne 
= D ômD(t, Sm) - ESm [Lx(Sm-1, Sm, TL)] - exp l- J vou] 
t 


m=1 
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Market quotes of collateralized OIS, IRS and TS and proper spline 
method allow us to determine 


{DO,T)}, 1E (L(Tn1, Tm, T)] + TSy} 


for all relevant T, Tm and tenor T. 


21.1.8.4 Cross-Currency Swaps 


When dealing with multiple currencies, we also need to consider CCS 
and FX forwards. The cross-currency formula for currencies i, j, with 
collateral in currency k, and collateral rate c and quote s is given by 


Ty 
N : ci 
E Sabet. Ta) -exp j= f Ods [E La Ts] 


n=l d 
To 
—D(c;, t, To) - exp - | oas 
t 
Ty 
+ D(cj, t, Ty) - exp - [ vio 
t 
(CCS) P ” 
= E mDj t Sn) -exp |= | Y Ods } ES L Spt, Sm) 


n=1 
t 


So 
—D(cj, t, So) - exp - | yon 
t 


SN 
+ D(cj, t, Sy) + exp - | oas 


t 
21.1.8.5 FX Forwards 


A collateralized FX forward rate at time ¢ for the currency pair (i,j) 
with collateral in currency k with collateral rate cg and maturity T is 
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defined as the solution ft? (t, T) of the following equation: 


T T 


0= fP, T)E® exp į — I ri(s)ds } - exp | yek(s)ds 


t t 


T T 
-FEDES exp - | ras < exp | oas > 
t t 


where an amount of currency i is exchanged for one unit of currency 
j, assuming a collateralization in currency k. Then we have 


[oo | Fone 
E;' | exp } f y/(s)ds 
pit, T) t 
iT) T i 
Pug) E" Ç [Fs 
t 


If the spread y is stochastic, the currency triangle, such as USD/JPY x 
EUR/USD = EUR/JPY, holds only for the same collateral currency. 
When y is deterministic, this becomes 


FMC, T) = fp 


T 


exp I y(t, s)ds 


t 


(i) aP GT id, D(,T) 
FDE T =O EDA OET 


which is independent of the choice of collateral currency. In a boot- 
strap all the previous except y'/ can be observed in the market. This 
can only be used for short maturities due to the lack of liquidity in the 
forward market for longer maturities. In the previous equation 


T 
7) B 
yit, T) = -zr in E" exp - | oas 


t 


21.1.9 Curve Calibration 


In this section, we describe our procedure for generating the curves 
needed. We fix t = 0 as today and we will write Eo['] = EL]. 
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Definition 21.1 By a collateral pair, we mean G, cj), where i is a 
currency and c; a collateral rate in currency i. 


Definition 21.2 Given a currency i and a collateral pair (j, cj), we 
call the mapping 


T 


Tr D(0,T) - exp [ "60s 
0 


a discount curve and denote it by D a0 T) 


Definition 21.3 Given a currency i and a set of collateral pairs 
[ (fj, i) } , we call the mapping 


j=] 
T 
Tr Da(0, T) - exp fe (y8(s)) ds 
j 
0 


a CTD curve and denote it by D j (0, T). 


Definition 21.4 Given a tenor x and a currency i we call the map- 


ping T > ET) [Lx (S; T)] a forward curve, where S, T are the start 
and maturity date, respectively. 


21.1.9.1 Single-Currency Collateral 


Given a currency i from the set of available currencies, a set of ten- 
ors (the number depending on i) and a collateral pair, our task is to 
generate the corresponding discount and forward curves. Due to li- 
quidity problems, we have to make some more assumptions regarding 
the quotes that we observe in the markets. Moreover, in some curren- 
cies the complete lack of OIS quotes forces us to make some sort of 
approximations. In this section we describe the currency-independent 
assumptions and we will come back to certain currency-dependent 
issues later. 


e ACCS quote observed in the market is interpreted to be of constant 
notional type. Mark-to-market (MtM) CIRS are not considered in the 
curve generation. 

e The CCS quote in Equation (CCS) is independent of the collateral 
pair (j,c j). 
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The last assumption makes it possible for us to compute for example 
EUR discounting under SEK collateral where the collateral rate is 
the Swedish repo rate (such a cross-currency quote is definitely not 
available in the markets). 


21.1.9.2 Stripping in EUR 


In EUR, the available instruments are OIS, IRS, FRA, CIRS and FX for- 
ward and the collateral pair is EUR, EONIA). The OIS and IRS quotes 
are assumed to be under (EUR, EONIA). 

We start to bootstrap the (EUR, EONIA) discounting curve using 
Equation (OIS) given earlier and then the forward rate curves for the 
tenors {1M, 3M, 6M} using Equations (FRA) and (IRS) 

These are the consistency conditions needed to get the market 
quotes of various swaps. We have denoted the market observed OIS 
rate, IRS rate and TS spread respectively as O/Sy, IRSy and TSy, 
where the subscripts represent the lengths of swaps. {T;};>0 are the 
reset/payment times of each instrument. We distinguish the day-count 
fraction of the fixed and floating legs by A and ô, respectively. These 
are not necessarily the same for different instruments. L(Tm-1; Tm, T) 
is the Libor with tenor t whose reset and payment times are Tm- 
and Tm, respectively. In Equation ARS), we have distinguished between 
two different tenors, ts and tz (> Ts), although we have used the same 
payment frequencies in the fixed and floating legs of the IRS. 

Ros is defined as 


-In (D(t, T)) 


Rors(t) = Tot 


For the forward Libor, the zero-rate curve R; is determined reclusively 
through the relation 


Te 1 e Pr (Tmi) Tmt) 
E ” p (t, Tmi, Tm I Om ( e`Rrı(Tm) (Tmt) 1) 


= 1 (= Tm-1) 1) 
Ôm D(t, Tm) 


= l (ee) 


bm Dt, Tn) 
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21.1.9.3 Stripping in USD 


In the United States, 30-day fed funds futures and 3-months OIS futures 
are commonly used, although there are only 36 standard 30-day FF 
futures contracts and eight standard 3M OIS futures contracts, based 
on www.cmegroup.com. Therefore, one can only bootstrap the OIS 
curve up to 2 or 3 years. 

In order to capture longer-term OIS discounting, one could then 
bootstrap the short end of OIS curve from fed funds OIS (and/or fed 
funds futures) and then bootstrap the long end through fed funds basis 
swap (FFBS) (typically combined with IRS to get Libor information). In 
the United States, the instruments are typically FFBS ranging from ly 
to 30y, which are averaged index basis swaps. This has become the 
common practice in the USD market. However, the actual implement- 
ation methodology could vary in several ways. We describe these in 
Method 1, Method 2 and Method 3 in the following sections. 


Method 1: Approximation Approach 


By ignoring discrepancy business day count adjustment and com- 
pounding crude adjustment (see Bloomberg (2012), Extending USD 
OIS curves using Fed Funds Basis swap quotes), we obtain 


OIS(t) approx. a 7 i) 


OIS(t)aaj © 4 - (1 r 


where 


ro- TBS. a 
4 


IRS(t 
ro=4-((1+ uaes an 4), 


where IRS(t) is the market quote of swap rate. 

FFBS(D is the market spread quote of FFBS. 

One can observe from market quotes that the implied OIS rate has 
approximately deterministic spread to IRS and such an OIS has embed- 
ded market information from FFBS. Moreover, the OIS is reasonably 
close to market OIS for long tenors, so using this adjusted OIS for 
longer (say, 30y) tenors is justified and computationally efficient. 

Thus this approximation method makes the following assumptions 
that practitioners must recognize. 
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Ignore conventions like business day count, calendar, roll conven- 
tion, basis, spot lag, schedule mismatch, etc. 

The compounding approximation assumes “flat curve with rate 
equal to the difference between Libor and Fed Funds Basis while 
ignoring weekends and holidays”. 

The error between adjusted OIS and OIS market quotes is a 
combination of approximation errors and liquidity. 


Thanks to the simplifying assumptions this method is very fast. 


Method 2: Brute Force Approach 


The brute force approach jointly solves both FFBS and IRS to par. We 
know analytics for both FFBS and the IRS, as shown in Table 21.4. 


Table 21.4 The fed funds basis swap and IRS 


Fed funds basis swap 


Pay Receive 

Libor leg Fed fund leg 

-3m Libor + weighted average Fed fund + basis spread 
IRS 

Receive Pay 

Libor leg Fixed leg 

+3m Libor — fixed swap rate 


Note: 


This is the most accurate approach because it solves for discount 
factors for both OIS discounting curve and Libor 3m forward curve 
jointly. 

Schedules among legs between FFBS and IRS may not align, which 
burdens the solver and is very sensitive to choice of interpolation 
method. 

Computationally expensive (especially the weighted average on 
daily forward fed funds effective rate feature). 


Method 3: Synthetic Approach 


The synthetic approach is to bootstrap the OIS curve by re-pricing 
FFBS to par given an IRS. The idea is represented as in Table 21.5. 


Through the synthetic construction in the previous graph, we can 


have a single instrument (fixed vs. floating (weighted average FF + 
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Table 21.5 Re-pricing fed fund basis swaps to par with IRS 


Synthetic fed fund swap 


Receive Pay 
Fed fund leg Fixed leg 
+ weighted average fed fund + basis spread - fixed swap rate 


basis spread)) per tenor, which would fit into most standard curve 
stripping frameworks as a regular bootstrapping procedure. 


Note: 


e This is faster than Method 2. 
e This is more accurate than Method 1. 
e Schedule alignment is still an issue. 


e Tenors have to be matched so as to create a synthetic FF average 
swap. 


21.1.10 The Bootstrap 


We now assume we have OIS, IRS, FRA, CIRS and FX forward and the 
collateral pair is (USD, FEDLO1) where FEDLO1 is the USD fed fund 
rate. The OIS and IRS quotes are assumed to be under (USD, FEDLO1). 
We start to bootstrap the (USD, FEDLO1) discounting curve as 
shown previously and then the forward rate for tenors {1M; 3M}. The 
6M forward rate curve can be built using the 6M-3M tenor basis swaps 
to create a synthetic 6M swap curve. 
Remarks: In the previous calculations for EUR and USD, we have 
assumed that all the instruments are collateralized by the cash of do- 
mestic currency, which is the same as the payment currency. You may 
worry about the possibility that the market quotes contain significant 
contributions from market participants who use a foreign currency as 
collateral. In fact, some of the major financial firms prefer USD cash 
collateral regardless of the payment currency of the contracts. This 
gives rise to additional factors in discounting as in 


T T 
nie) = E24 es l- f rod - exp i f rod ma) , 


T 
= D(t, T)E/4 [es l (y4(s) -y (8)) al ma) . 
0 
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Changes in IRS/TS quotes are very small and impossible to distinguish 
from the bid/offer spreads in normal circumstances, because the cor- 
rection appears both in the fixed and floating legs which keeps the 
market quotes almost unchanged. However, the PVs of off-the-market 
swaps can be significantly affected when the collateral currency is 
different. 


21.1.10.1 Stripping in SEK 


In SEK, the available instruments are also OIS (STINA swaps, which are 
not real OIS, but T/N), IRS, FRA, CIRS and FX forwards. The collateral 
pairs are (SEK, STIB1D) and (SEK, SWRRATE). Here STIB1D is the T/N 
rate (STINA) and SWRRATE the Swedish repo rate. The OIS quotes are 
assumed to be under (SEK, STIB1D) and the IRS quotes are assumed 
to be under (EUR, EONIA). 

The assumption about the IRS quote makes the bootstrap procedure 
a little bit different from earlier. We begin to bootstrap (SEK, ST IB1D) 
as usual for OIS. Then we need to add a cross-currency formula so 
we simultaneously can bootstrap the discounting curve for SEK with 
collateral pair (EUR, EONIA) and the SEK 3M forward curve. 

The (SEK, SWRRATE) is generated as follows: 


1. Bootstrap a curve from available SEK RIBA (Riksbanks futures) in- 
struments and compute the daily instantaneous forward rates up to 
the maturity of the last RIBA instrument. This will be the short end 
of the resulting curve. 


2. Compute the spread between the last instantaneous forward rate 
in the first step and the corresponding rate generated from (SEK, 
STIB1D). 


3. Extrapolate the spread and convert the resulting instantaneous 
forwards to a discount curve. 


The reason we also have to bootstrap the Swedish repo rate is that it 
is used as collateral rate in most CSA agreements. 


21.1.10.2 Different Problems with Bootstrapping 


The Choices of Instruments 


The first problem is related to the availability and choice of liquid in- 
struments. Let us assume that we have an OIS curve (obtained from the 
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previous method) and we need to find a forward curve (e.g. USD Libor 
3m curve). In the single-curve world, we typically bootstrap from com- 
binations of cash (1w up to 1y), futures or FRAs (typically most liquid 
up to 3y or 5y) and swaps (ly up to 50y). However, in the dual-curve 
world, each Libor has its own risk, which means we cannot use spot 
Libor cash rates like 1w, 2w, 3w, 1m, 2m, 4m, 5m, 6m, etc., to strip the 
3m Libor curve. Therefore, one possible choice could be to use FRAs 
for short term and swaps for long term. Note that the use of FRAs 
is subject to a few problems. In the dual-curve world, the Libor be- 
comes risky, which means that the old approach would use the wrong 
discounting curve. 

The problem here is that trading parties can only put daily collateral 
until the start time of the FRA period. Therefore, the period between 
the start time and the first maturity date of the FRA is only for theor- 
etical Libor tenor, and there won't be any collateral transactions. This 
simple approach would assume that Libor forward is risky when it is 
settled, though any such settlement amount would still be discoun- 
ted at the risk-free rate (e.g. continuously full collateralization with 
zero MTF). For futures, it is even trickier in view of the convexity 
adjustment. 


The Libor Rates 


In the single-currency world, the standard curve construction would 
be to make a discount curve and then an implied forward curve. This 
implies that the Libor forward depends on the choice of discounting. 
However, intuitively Libor is a floating index that should not depend 
on the choice of a discount curve. Different players may have different 
funding curves, and CSA may also require discounting other than OIS 
discounting which really refers to cash collateral within that single- 
currency market. 

Mathematically, the curve construction depends on the discount- 
ing method in our typical risk-neutral environment. This is because 
the probability measure changes when the discounting method is 
changed. Prior to the crisis, the “standard” discount curve was the 
“same” as the Libor forward curve in USD. Now, the “standard” dis- 
counting method has shifted into OIS discounting. OIS discounting 
and Libor discounting would produce nearly the same price for a 3m 
Libor forward. 
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Furthermore, Libor fixings too have shifted to OIS discounting, 
although dealers/brokers may not report their quotes consistently. 
As CSA regulation tends to standardize collateral agreements and as 
the LCHuses OIS discounting for clearing swaps, the swap quotes 
and Libor fixings would be more reliable than during the crisis 
period, when Libor was criticized due to various reasons including 
collateral/funding discounting and liquidity. 


Collateral 


Another concern is collateral currency. This is actually a fairly com- 
mon problem in the real world. Suppose a EUR-based bank is trading 
collateralized vanilla US IRS with a US bank. Assume further that most 
of the EUR-bank’s business is in the EUR and that the EUR-based bank 
owns a lot of EUR treasury bonds and cash in EUR. Therefore, it would 
be very natural for the EUR bank to put collateral in a EUR cash or EUR 
treasury bonds, rather than put collateral in USD even if the trade is a 
USD single currency. In this situation, we are faced with the dilemma 
of how to discount the value of the trade price and the risks. This is 
the “cheapest to deliver” problem. 


Cross-Currency swaps 


There are two popular types of CCS. 


e Float-float (commonly used for major currency pairs like USD/EUR, 
USD/JPY, EUR/JPY) 


e Fix-float (commonly used for minor currency pairs like USD/TWD). 


We will focus mostly on float-float CCS (CCBS). Using USD vs. JPY as 
a simple example to illustrate arbitrage-free relationships in interest 
rate and FX markets (“interest rate parity”), let us suppose a JPY in- 
vestor has a future cash flow of 100 million in USD ($100 MM) at future 
date T. She/he could either discount it using USD discounting back to 
today t = 0 and convert it into JPY using the spot FX rate (FX(O)), or 
she/he could convert the future cash flow at T through FX rate (FX(7)) 
and then discount back to today using JPY discounting. If there were 
no-arbitrage opportunities, these two methods would yield the same 
amount. Such parity relations can be visualized in Fig. 21.17. 
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100MM USD FX(T) 
DF(USD) DF(JPY) 
—_—_—_—_—_——>> 
FX(0) 


Fig. 21.17 A parity relation on cross-currency 


Domestic Yield Curve 
“Ss FX Forward 
Foreign Yield Curve 


Fig. 21.18 With two known the third can be solved 


Domestic Yield Curve 


Foreign Basis Curve 


FX Forward 


Fig. 21.19 Bootstrap of an implied foreign yield curve 


Therefore, in theory, we should see a perfect triangle relationship 
betweeen the domestic yield curve, the foreign yield curve and the FX 
forward. As we illustrate in Fig. 21.18, one can solve for one out of 
three given the other two. 

However, both interest rates (domestic and foreign) and FX forwards 
are independent actively traded markets. Therefore, the triangle won’t 
be “perfect” even in normal market situations due to differences in 
liquidity among those three markets. 

If one needs to value a cross-currency trade under the domestic 
measure, then a typical way would be to bootstrap the “implied” for- 
eign yield curve (the foreign basis curve) from a set of FX forwards 
given a domestic yield curve as in Fig. 21.19. 
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Note that such a foreign basis curve typically won’t be same as the 
foreign local yield curve (e.g. bootstrapped from cash, FRA/future, 
swap). This reflects the fact that the foreign investor’s (domestic 
currency’s) funding interest rate (discounting) at the local (foreign cur- 
rency) market typically will be different from that of a local investor. 
Such a basis curve would satisfy 


e Re-pricing each node of FxFwd back to par. 


So far, everything is still consistent. Such a relationship implies that 
the cross-currency basis is zero (which is what “perfect” means here). 
However, if we look at CCBS quotes, they are not zero! This im- 
plies that the CCBS actually breaks the previous interest rate parity 
due to many reasons, supply and demand factors, including liquid- 
ity premium, currency strength, currency country’s credit profile 
differentials, etc. 

Prior to the crisis, such bases were very small, and many practition- 
ers simply ignored them . However, after the crisis these bases became 
much larger (50-100 bps, compared to 1-5 bps pre-crisis). Therefore, 
such a basis can no longer be ignored in derivatives pricing and risk 
management. As with FX forward, one can also bootstrap an “implied” 
foreign basis curve from a set of CCBS given both domestic yield curve 
and foreign yield curve (which makes Libor forward consistent with 
local Libor forward, see Fig. 21.20). 

Such a foreign basis curve should satisfy 


e re-pricing each node of CCBS back to par; 
e matching local Libor forward. 


Domestic Yield Curve 
Foreign Yield Curve 


Cross Currency Basis Swap 


Fig. 21.20 Bootstrap an “implied” foreign basis curve from a set of cross-currency 
basis swaps 


Foreign Basis Curve 
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In short, we solve for the discount curve while keeping the prevailing 
index forward unchanged. Therefore, the Libor definition is changed 
slightly due to different discounting methods. From a practical point 
of view, this problem can often be ignored due to the tiny effect on 
forward change which means one can still assume that the basis swap’s 
Libor forward in foreign currency is the same as Libor forward in the 
domestic market. As in the single-currency world, even if we change 
from single curve Libor to a very different discount curve (e.g. the OIS 
curve) when pricing a vanilla IRS, the Libor forwards in the IRS would 
be almost the same (maximum about 1-3 bps for 30y). 

Furthermore, FX forwards may only last up to 5 years, depending 
on the currencies, while cross-currencies basis swaps would last up to 
30 years. Therefore, from a liquidity point of view, one may strip the 
foreign basis curve from a combination of FX forward (for short term) 
and CCBS (for long term) (See Fig. 21.21). 

OIS curve 


We now turn our attention from the traditional problem in the cross- 
currency world to the problems arising in a multi-curve environment. 

Recall that within the single-currency framework, the dual-curve lo- 
gic is to bootstrap the OIS curve, then solve for the index forward 
curve while keeping the OIS curve unchanged. In the cross-currency 
world, the logic is to solve for the implied foreign discounting curve 
while keeping the foreign Libor forward unchanged. Thus, in a CCBS 
there are four curves involved: 


1. Domestic discount curve (e.g. OIS curve). 
2. Domestic forward curve (e.g. 3m Libor curve). 


Domestic Yield Curve 


Foreign Yield Curve 
FX Forward 
Cross Currency Basis Swap 


Fig. 21.21 A foreign basis curve stripped from a combination of four sources 


vq Foreign Basis Curve 
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3. Foreign discount curve. 
4. Foreign forward curve (e.g. 3m Libor). 


Still the goal here is to solve for a “cross-currency implied” foreign basis 
discount curve (or say a foreign OIS basis curve) where discounting 
and forward projections are decoupled. 

The foreign OIS basis curve can be bootstrapped: 


1. From FX forward and the domestic OIS curve 
2. Using domestic and foreign OIS curves and CCY OIS basis swaps 
3. Using CCY Libor basis swaps 


In point 1 we assume that both domestic and foreign FX forward mar- 
kets move to OIS discounting. Point 3 is more complex due to the 
market conventions. The common practice would be the following 
steps (using USD 3m Libor vs. JPY 3m Libor as an example): 


1. Construct the domestic OIS discount curve. 
2. Construct the domestic index forward curve (e.g. USD 3m Libor). 


3. Construct the foreign OIS discount curve (e.g. JPY OIS curve from 
Tonar). 


4. Construct the foreign index forward curve (e.g. JPY 6m Libor). 


5. Construct the foreign index forward curve (e.g. JPY 3m Libor from 
Libor basis swap). 


6. Solve for implied foreign basis curve (or, say, foreign OIS basis 
curve) from CCBS given in 1, 2 and 5. 


Cheapest to Deliver 


CTD is driven by the collateral agreement and CSA standardization. 
In CSA, the two counterparties may have rights to choose collateral 
on the fly among a few predefined currencies. This will create a huge 
headache for modelling, hedging and risk management. 

Most dealers agree that the discount rate should - in theory - be 
based on the CTD collateral. Market practice is altogether different, 
with even the major dealers taking a variety of approaches to pricing 
trades based on multi-currency CSAs. 
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“Theoretically, it is not difficult to put together a model - it would 
be an extension of a stochastic basis model where you have more than 
one basis,” says Vladimir Piterbarg, global head of quantitative research 
at Barclays Capital in London. “The huge question is whether you are 
able to execute the hedging strategy required.” 

Such a disorder in CTD modelling and hedging would have much 
bigger model risks. CSA standardization and simplification become 
critical while the market is moving towards a consistently multi-curve 
approach. 


21.1.10.3 Collateral Choice 


To compute the integral in the mapping 
T 
T +> Da(0, T) - exp J max (yi) ds}, 
j 
0 


we assume for each i and f; that y(t) is a piecewise constant function 
of t. This implies that the integral can be written as 


T K Tk 
J max (y w) ds = a max J yh (s)ds 
j j 
0 k=1 Tk 
where 0 = Tọ < Ty < ... < Tg = T is the approximation 


scheme where all the yili(t) are piecewise constant. Moreover, from 
the definition of the discount curve we can see that 


Tk 
| yi (sds = lIn (D.:(0, Tx)) 
Tk 


-In (D4 (0, 74) = In (Da (0, Te-1)) + In (D4 (0, Ti) 


and hence the computation is now straight forward given the relevant 
discount curves. 
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21.1.10.4 Stripping in HUF as an Example of Illiquid 
Currencies 


For example, say we are looking for curves in HUF (Hungarian For- 
int) rates. Then we can use local swaps, that is, HUF3M vs. fixed and 
CCS HUF3M vs. EUR3M. At the same time, we have to consider the 
collateral assumptions behind this. 

We then use vanilla swap HUF3M vs. fixed, quarterly HUF3M (3M 
BUBOR) annually fixed and CCBS EUR3M vs. HUF3M where the EUR 
floating leg is quarterly EUR3M (EURIBOR3M) and HUF floating leg, 
quarterly HUF3M + spread (Tables 21.6 and 21.7). 

When we do simultaneously stripping of the discount curves and 
the projection curves we use the following stripping equation for the 
1Y point: 


Q-5-D4=6,-F,-D,+62-Fo-D2+63-F3-D34+64-F4-D4 
6): (F1 +s) - Di + 62+ (Fo + s) - D2 + 53+ (F3 + s) - D3 + ô4 - (F4 + s) - D4 


_ EURLeg 
— XO 
where 
D; = Discount factors 
F; = HUF3M forwards 
6; = Year fractions 


(O) = Spot EUR/HUF exchange rate 
Q = Quoted 1Y par swap rate 
S = Quoted 1Y EUR3M/HUF3M basis spread 


We also assume that the EUR curves are already stripped (e.g. from 
EONIA swaps and vanilla EUR3M swaps). So the unknowns in the 


Table 21.6 Vanilla swap HUF3M vs. fixed 


Discount curve Projecting curve 
Fixed leg HUF5< curve N/A 
Floating leg HUF5 curve HUF3M curve 


Table 21.7 Cross-currency basis swap EUR3M vs. HUF3M 


Discount curve Projecting curve 


EUR leg EONIA curve EUR3M curve 
HUF leg HUF5< curve HUF3M curve 
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previous formula are only D4 and D4. The intermediate D, and F, Gf 
unknown) are computed via interpolation. This can also be handled 
by a general solver. 

In a practical situation, things might need to be done differently. The 
payment dates of the local vanilla swaps and those of the CCBS might 
be misaligned (due to differing day count conventions), or simply the 
quoted maturities for one set of instruments might be different from 
the quoted maturities for the other set of instruments. Thus perform- 
ing a bootstrap might not be the best solution. Instead, we could use a 
global solver for all quoted instruments simultaneously. This is slower 
but produces more stable results. We could still perform intermediate 
passes using the quotes up to some fixed maturities in order to find 
good initial guesses for the later passes. 

Simultaneous stripping produces two curves, a discount curve 
HUF®® and a projection curve HUF3M. By §construction, using these 
two curves as the discount curve and the projection curve, respect- 
ively, we will price at par both the vanilla swaps HUF3M vs. fixed and 
the CCBS HUF3M vs. EUR3M. 

Now, we might ask how this curve stripping does fit into a general 
model for derivatives pricing. So far we have only considered “linear” 
instruments and defined formally discount factors and forwards. Can 
these be used to price something else but swaps? Is this backed by a 
theory where the curves get back their usual meaning? 

Consider an economy with two currencies: a domestic (Dom) and 
foreign one (For). Assume that collateral can be posted in any of the 
two currencies. The choice of the collateral currency holds for the 
whole lifetime of the derivative G.e. assuming there is no option to 
switch collateral), so the domestic collateral earns foreign collateral 
earnings cf. A pricing theory can be constructed rigorously with the 
help of the following replication argument. 1° 

Denote the price of a collateralized domestic derivative by V. Then, 
in the domestic collateral 


T 
V(t) = E; | exp - f cuas VD|F; 


t 


10 See Piterbarg (2010). 
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In the foreign currency we have 


T 
V(t) = E; | exp -J ca(s) + h(s)ds } - V(T)|F; 


t 


Similarly, a collateralized foreign derivative VÍ is given by 


T 
VO = E exp -| [cp(s) —h(s)] ds} - VDF, 


t 


and 


T 
Vin) = Ef exp - | oas VI(T)F, 


t 


The fact that the spread when computing from the foreign point of 
view is —h follows from the domestic-foreign “parity” condition. If we 
fix the collateral currency and then compute the price of a contin- 
gent claim through foreign or through domestic yields we must get 
the same result in order to exclude arbitrage (see Fig. 21.22). 

This implies that the drift of the FX rate X (in the domestic 
measure E) should be rg, = ca- cp +h. 


Time T 


Domestic Foreign 


Fig. 21.22 Foreign or domestic yields we must give the same result to exclude 
arbitrage 
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Table 21.8 Pricing in classical theory and collateral 


Domesticrate Foreign rate FX rate drift 
Classical theory lg rf 'd-rf 
Domestic collateral Ca Cig = C¢-—h rag = Cg- Cs +h 
Foreign collateral Cop = Ca +h Cy lop = Cg- Cf +h 


Table 21.9 Vanilla swap HUF6M vs. fixed 


Discount curve 


Fixed leg HUF®5¢ curve 
Floating leg HUF®< curve 


Projecting curve 


N/A 
HUF6M curve 


Table 21.10 Tenor basis swap HUF6M vs. HUF3M 


Discount curve 


Fixed leg HUFS< curve 


Projecting curve 


HUF3M curve 
HUF6M curve 


Floating leg HUF@®< curve 


Table 21.11 Cross-currency basis swap EUR3M vs. HUF3M 


Discounting curve Projecting curve 
EUR leg EONIA curve EUR3M curve 
HUF leg HUF“5¢ curve HUF6M curve 


Once the collateral currency has been chosen, pricing under a CSA 
is in some way the same as pricing in the classical theory, as long as 
the appropriate curves are used. 

Hence we have a pricing theory that is consistent and extends the 
formal swap pricing theory based on stripping where the stripping 
produces the initial term structures of the rates, that is, today’s values 
of the curves used for discounting and forwarding. As seen before, 
the choice of the collateral is reflected in the curves that are used for 
discounting. 

A problem will occur if the market does not directly quote 3M 
swaps. If we have quotes for 3M vs. 6M tenor basis swaps then we 
can simultaneously strip three curves (Tables 21.9, 21.10 and 21.11). 


e Local swaps: HUFOM vs. fixed 
e Local TS: HUFOM vs. HUF3M 
e CCS: HUF3M vs. EUR3M 
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In the aftermath of the financial crisis in 2007 and 2008, the market 
has turned towards OIS discounting as the new standard. The G5 cur- 
rencies and a few others have liquid OIS markets which provide the 
quotes from which the OIS curves in those currencies can be stripped. 
But for the other currencies the OIS markets are more often than 
not inexistent or illiquid. This renders it impossible to strip out any 
OIS curve. There is no magical solution, but one can contemplate 
using available cross-currency quotes to infer a number of curves. 
This requires the simultaneous stripping of single-currency and cross- 
currency instruments. The discounting method is intimately tied to the 
choice of the collateral; hence there is a trade-off: not having a curve 
at all vs. using curves based on a different collateral assumptions. 


21.1.10.5 The Future — Standard CSA (SCSA) 


The problem with the aforementioned method is the embedded op- 
tionality in CSA discounting to post collateral in many currencies 
and the choice of instrument to post. This leads to the following 
problems: 


e no price transparency 
e difficult to do unwinding and novation of trades 
e difficult to hedge 


Many market actors believe in using the Standard Credit Support An- 
nex (SCSA), that is, a single-currency CSA. Here one will probably use 
the following: 


e 17 currency silos 


e Emerging currencies will use the multi-currency USD silo. For these 
currencies, we have no domestic OIS market. Therefore, we will 
use a simultaneous calibration of USD-collateralized domestic IRS 
and CIRS. 


This would eliminate the currency switch options and we could use 
standardization on OIS in order to get better alignment with Central 
Clearing Partners. The institutions can make the move from CSA to 
SCSA but not vice versa. The switch to SCSA will not be mandatory. 
A greater collateralization will reduce the counterparty risk but also 
increase the funding costs. 
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21.1.11 General Pricing in the New Environment with 
Funding Value Adjustments 


We will now describe general pricing in the new environment, with 
collateral and funding. Here the FVA can be considered as the differ- 
ence between the PV of the trades between the “real-world” and the 
theoretical CSA: FVA) = PV (real-word CSA) - PV (theoretical CSA). 
This is the cost above the risk-free price to fund the trades in the pres- 
ence of a real-world CSA. We can also consider the FVA as the cost to 
fund the trades above that of the risk-free price 


FVA(t) = PV(including funding costs) — PV(risk-neutral). 


There is still, in 2014, an ongoing debate on how to calculate and 
report the FVA and how to pass this cost to the counterparty. An- 
other unresolved question is whether the FVA should be included in 
MtM prices or just be included in the unrealized profit/loss (P&L). In 
addition, it is not clear what the regulators will decide. 

Market practitioners recognize that collateral has been successfully 
introduced in order to mitigate CVA. However, since collateral affects 
the funding of the trade, it also affects profitability. There is no argu- 
ment that the FVA charge is real, and passed onto traders by funding 
desks. 

In the classical theory we have postulated the existence of a single 
rate r(t) where all tradable assets grow at this rate under a risk-neutral 
measure. The portfolio value v(t) satisfied the standard Black-Scholes 
partial definition equation (PDE) (or term-structure equation) 


Lit)v(t) = rv), 


where LC) is the evolution operator of the model drivers. All payments 
are discounted with the single rate r(t). The modern market on the 
other hand has multiple rates, for example: 


rc the collateral rate (which is almost risk free) 
rr the interest rate for asset secured borrowing (the repo rate) 
rr the interest rate for unsecured funding 


For example, a fully collateralized payment is discounted with rc. 
Now, we will study what happens if the collateralization is only partial 
and how this theory is constructed. The evolution of the aforemen- 
tioned rates is modelled by correlated stochastic processes and the 
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portfolio may consist of general instruments, both vanillas and exot- 
ics, with known future payouts as function of the rates. Furthermore, 
the collateral is a known function of the portfolio value. 

The FVA is defined as the difference between the modern price 
and the base one. The base price is often related to a fully collater- 
alized portfolio. This is the only difference between the modern and 
the base models; the payments/indexes coincide. Only modern price 
is important; the base can be chosen from pragmatically numerical 
reasons. 

Piterbarg!! calculated a modern portfolio value by replication. This 
replication gives a modified pricing PDE (with respect to the classical 
theory) and is a unique way to determine the price. However, at this 
stage, the default risk (of the bank or the counterparty) is not taken 
into account. 

Let V(t) be the portfolio price and C(t) the collateral. The collat- 
eral is a known function of the portfolio (e.g. C(t) = (V(t))*). Then by 
replication we have 


LØVE) = rcC@) + rr(VE) - CH) = rr VE) + CO - rr), 


where LŒ) is the evolution operator corresponding to the rate pro- 
cesses. We also assume that the equity grows with the repo rate rp(t) 
minus dividends. In the right-hand side we divide the portfolio value 
in two parts V = C+ (V - C) where the part under collateral, 
C, is discounted with the collateral rate rc and the residual part, not 
covered by collateral, V — C, is discounted with the funding rate rp. 
One can generalize the theory to different rates for positive/negative 
parts of the collateral and funding (see Pallavicini et. al.). This would 
reflect a distinction between borrowing and lending. 


11 Vladimir Piterbarg (2010), “Funding beyond discounting: collateral agreements and derivat- 
ives pricing”, RISK, Feb. 
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The solution to the PDE between payment/exercise dates is ob- 
tained via the Feynman-Kaé theorem 


T 


Vit) = EÈ | exp t — J rp(s)ds } - VD|F; 


t 
T T 


+E? [oe - [ress (rru) —rc(u)) C(u)dulF; 


t t 
or equivalently 


T 
V(t) = EÈ | exp - | rewas VD|F; 


t 
T T 


+E? ] exp 4 — l re(s)ds ¢ (rcu) — rr(u)) Vu) - C(u)) dul F, 


t t 


The previously mentioned PDE is non-linear (except when collateral is 
a linear function of the price V). We also have some special cases: 
For C = 0 we get a solution for a non-collateralized deal 


T 


V(t) = E£ | exp } — I rp(s)ds ¢ - V(T)|F; 


t 
For C = V we get a solution for the fully collateralized situation 


T 
V(t) = EL exp - f rcas -V(T)|F;, 


t 


Notice here that a fully collateralized security is discounted with the 
collateralized savings account, while the uncollateralized security is 
discounted with the pure funding account. Note that both condi- 
tional expectations are taken with respect to the same measure. Also, 
some conditions must be imposed on the coefficients in the evolution 
operator L(t) in order to avoid explosive solutions. 
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Of course, one can change the pricing measure using, for example, 
the Radon-Nikodym derivative 
M(t) = D(t,T)/Bc(), 
where D(t, T) is a collateralized zero-coupon bond 
T 
D(t, T) = EÈ | exp -J re(s)ds } |F; 
t 


The collateralized deal propagation leads to a standard formula 


T 
Vc(t) = E2 exp - f rewas -Vc(T)|F; 


t 


T 
1 
= M(t). ETC exp - f rewas -VC(T)- un 
t 
= _Btc Vc(T) 
ears Fa m | 


= D(, T) - E": |[VeDIFi] 
However, for the pure funding security, such a change of measure gives 


T 
Vr(t) = EL exp - f roas -Ve(T)|F; 


t 
T 


1 
= M(t)- ETc exp -J rr(s)ds ¢ - VF(T) - mn 


t 


and will not lead to a standard numeraire-deflated form. 


21.1.11.1 Portfolio FVA 


One way to calculate the portfolio FVA can be summarized as follows: 
First, compute future values (exposures) for each instrument in the 
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portfolio using the single-curve model and then aggregate of the in- 
struments at the portfolio (netting set) level and compute the FVA with 
the following approximate formula: 


T t 


FVA(th=E J (reg(t, V(t- r(t)) exp -f rept (s, V(s))ds ¢ - V(t)dt 
0 0 


Here rez(t, V(t) — r(t) is the effective funding spread over the OIS rate 
and V(t) the future exposure (MtM) of the portfolio at time t. The last 
integral represents an effective funding discount factor. The effective 
rate is given by 


reg (t, VO) - Vt) = rcC() + rr(V(t) - CO). 
This gives us the price, including the funding costs as 


T 
Vip) = EL exp -| ref (Ss, V(s)) ds ¢ - VTF; 


t 


The important point to note in the previous formulas is that V(t) is the 
uncollateralized exposure, and not the total exposure of the deriv- 
ative. The average is taken over all possible future paths where both 
the derivative and the collateral change due to market movements. 

Intuitively you can think of the FVA formula as averaging the amount 
of the derivative trade that cannot be funded by the collateral. Also, in 
order to ensure that you will be able to meet your future cash flows, 
you must have access to funding at your external or market funding 
rate, which will typically be higher than the OIS rate prescribed in 
the CSA. 

We now define the effective funding rate as 


ref D = ec) - re) + - Get) - rrO), 


where gc is the fraction of traded value in the collateral account and 
1 - g,the fraction of traded value with credit exposure. With 100% 
collateral øc = 1 and reg the collateral (OIS) rate. With 100% funding 
(c= 0 and reg the funding rate. 
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The effective rate is a path-dependent quantity, due to the collateral 
changing depending on the prevailing rate environment - or if a secur- 
ity is used as collateral - due to the market fluctuations of the collateral 
itself. 

Over the lifetime of the trade, you can calculate the “effective fund- 
ing rate” on each path, which is a weighted average of the collateral 
rate (OIS) and the firm’s market funding rate, weighted by the amount 
inside and outside the collateral account respectively. Again, the in- 
tuition is that the collateral account cannot be 100% relied upon to 
ensure that you could meet all future cash flow obligations, due to 
the imperfect nature of the collateral arrangement with CTD options, 
MTM, thresholds, etc. The amount needed above the amount present 
in the collateral account will have to be funded in the market. This 
effective rate gives the average funding rate taking into account both 
sources of funding. 


21.1.11.2 The Law of One Price 


If a corporate (price-taker) gets bids from a number of banks, he/she 
will certainly get different quotes. These quotes may differ up to 100 
bps or thereabouts. Therefore FVA is a real charge to traders reflect- 
ing their diverse funding costs. Traders are passed this charge from 
their treasury/funding desk. 

Traders will attempt to pass this charge onto the purchaser as an 
upfront fee as a spread over the floating rate. However, this will 
be difficult in competitive markets. The FVA will certainly affect the 
profitability of the trade and all traders must know this charge. 

In the interbank market, where the price is determined by the law 
of supply and demand, the prices on average should not include any 
charge for FVA. 

In 2013, four banks (RBS, Lloyds, Goldman and Barclays) reported 
funding valuation in their financial statements. 


22 


CVA and DVA 


22.1 Credit Value Adjustments and Funding 


For years, a practice in the financial industry has been to mark de- 
rivatives portfolios to market without considering counterparty risk. 
All cash flows were discounted using the LIBOR or another “risk-free” 
curve. However, the true portfolio value must incorporate the possib- 
ility of losses due to counterparty default. This observation has gained 
wider recognition following the high-profile defaults of 2008. The 
credit value adjustment (CVA) is by definition the difference between 
the risk-free portfolio and the true portfolio value which should take 
into account the possibility of counterparty defaults. In other words, 
CVA represents the monetized value of the CCR. 


22.1.1 Definitions of CVA and DVA 


When reporting the fair value of any derivative position, we also 
need to consider counterparty credit risk (CCR). This is done by an 
adjustment to the value, known as CVA. 

A pure definition can be written as 


CVA = Discounted expected exposure x Default probability 
x Loss given default. 


For symmetry reason, we also need to consider the bilateral nature 
of CCR. This means that an institution would calculate a CVA under 
the assumption that they, as well as their counterparty, may default. A 


© The Author(s) 2017 607 
J.R.M. Roman, Analytical Finance: Volume II, 
https://doi.org/10.1007/978-3-319-52584-6_22 


608 J.R.M. Roman 


defaulting institution “gains” on any outstanding liabilities that cannot 
be paid in full. This component is often referred to as debt value 
adjustment (DVA). 

The justification for using DVA is that a bank could buy back its debt 
cheaply to realize DVA gains, but no firm has actually done this. 


22.1.2 Standard Approach 


To study CVA and DVA we consider the simplest deal we can think of, 
a loan where the premium paid for an amount K at time T is equal to 
the risk-free price minus CVA. This is 


P=e"'K-CVAL, 
where 


CVA, = Efe"? K - lpqpery] =e" K-Olt(B) < T] = eK: i z go"), 


Here, we denote the lender with L. t(B) is time to default for the 
borrower (B), and x (B) the credit default swap (CDS) spread for 
the borrower. Here we used the definition of the default probability 


Ont) > D= E [irg] = aA 
> 
If we have a certain recovery rate R, the previous formula should be 
1- e7 BT 
Qp(t(B) < DE [tæn] = — p 
For the borrower we have 
P =—e”TK + DVAg. 


The present values of all cash flows are given by 


V =e"! K —CVA,—P; 
Vg = -e "TK + DVAp + P. 
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To get an agreement between the borrower and the lender we need 
Vr = Vg. We then have 

2P = 2e"' K — CVA, — DVAg. 
With DVA; = CVAg we get 
P=e”TK et È _ gren] = K . e B)T, 

This is the agreed price between the borrower and the lender. But 
something is missing here. We need to know about the liquidity. 
The lender needs to finance the amount until maturity at a funding 
spread S(L) but the borrower can reduce his funding cost with P. 
The borrower should therefore see a funding benefit, and the lender 
should see that the fair value of the claim reduced by the funding cost. 


Therefore, the effect of this funding cost seems to be missing in the 
earlier formula. 


22.1.3 Approach Including Liquidity 


One method is to introduce DVA as a liquidity cost by adjusting 
the discounting term and introduce the possibility of defaulting on 
the payoff. This means that we need to replace the term e”TK with 
errr: líz(B)>T} and get 


VL=E ieee I(n@=r) | =e peek Ty 1B)>7)] -P 


= e Tr D+ LT .K. e TT -P= eta (Ltn (B)+y (LT .K-P. 
(22.1) 


Here y (L) is the deterministic default intensity for the lender, L 
s(X) = W(X) + y(X); X = {B, L}. 
Similarly, for the borrower B we have 
Vg = -E| 0T. K. emr) +P 


2P a K. I(n@=r)| +P 
-+7 (B)+Y(B))T | g. o ABT y p = e 2r BHyBT . KK 4 p 
(22.2) 


= 
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When comparing this result with the previous section, it is conveni- 
ent to confine ourselves to the simplest situation where the lender is 
default-free with no liquidity spread, while the borrower is defaultable 
with the minimum liquidity spread allowed, the CDS spread. In this 
case we have (s(L) = 0, s(B) = z (B) > 0) and we get 


Vi = etm (BT K E P 


Vg = e VAT BT g + P. 


Here we observe two bizarre aspects. First, even in a situation where 
we have assumed no liquidity spread the two counterparties cannot 
agree on the simplest transaction with default risk. A day-one profit 
should be accounted for by borrowers in all transactions with CVA. 
This contradicts market reality. 

Secondly, the explicit inclusion of the DVA term results in a duplic- 
ation of the funding benefit for the party that assumes the liability. 
Against all evidence the formula implies that the funding benefit is re- 
munerated twice. If this were correct then a consistent accounting 
of liabilities at fair value would require pricing zero-coupon bonds 
by multiplying twice their risk-free present value by their survival 
probabilities. 


22.1.4 How to Make It Right 


To solve this in a right way, we do not calculate the liquidity by the 
adjusted discounting approach as in Equations (22.1) and (22.2). In- 
stead, we generate the liquidity costs and the benefits by modelling 
explicitly the funding strategy. Here we take into account how the 
companies capitalize and discount money with the risk-free rate r, and 
then add or subtract the actual credit and funding costs that arise in 
the deal. 

This allows us to introduce explicitly both credit and liquidity and to 
investigate more precisely, where credit/liquidity gains and losses are 
financially generated. We take into account that the aforementioned 
deal has two legs. 

If we consider the lender, one leg is the deal leg, with net present 
value 


E[-P+e""Tl]. 
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Here II is the payoff at T, including a potential default indicator. The 
other leg is the funding leg with the net present value 


E[+P-e"'F], 


where F is the funding payback at T, including a potential de- 
fault indicator. When there is no default risk or liquidity cost in- 
volved, this funding leg can be overlooked because it has a value of 
E [+P age "P| = Q. In the general case the total net present value is 


Vr =E|-P+e "T+ P-e F] =E[e"(H-F)]. 


Thus, the premium at time 0 cancels out with its funding and we are 
left with the discounting of a total payoff including the deal’s payoff 
and the liquidity payback. An analogous relationship applies for the 
borrower, as will be described later. 

When we continue, we work under the hypothesis that all liquidity 
management happens in the cash market. Then, the funding is made 
by issuing bonds and excess funds are used to reduce or to avoid in- 
creasing the stock of bonds. This is the most natural assumption since 
it is similar to the assumption that banks make in their internal liquid- 
ity management, namely, what the treasury desk assumes in charging 
or rewarding trading desks. 


22.1.4.1 Risky Funding with DVA for the Borrower 


The borrower has a liquidity advantage from receiving the premium P 
at time zero, as it allows him/her to reduce the funding requirement 
by an equivalent amount. The amount P generates a negative cash flow 
at T, when funding must be paid back. This is equal to 


P. eT e(BT i lrB)>T}. 


This future outflow equals P, capitalized at the funding cost, times 
a default indicator. We need to include a default indicator in case of 
default and zero recovery. During default, the borrower does not pay 
back the borrowed funding and there is no outflow. Thus reducing the 
funding by P corresponds to receiving at T a positive amount equal to 


P- eTe®T 1 uwr) = {s(B) = 2(B) + y(B)} 


=P. eTa PTD. 1 sr). (22.3) 
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Thus, the total payoff at Tis 


Peg eet gay. lirB>T} -K> lirB>T} 


Taking discounted expectation, we get 


Vg =e” T p.e" BT eY BI _ oT, g.e BI = p. eY BT _K. e0 BT 


Compared with Equation (22.2), (P — K - e+? ®)+7/(®))T) we have no 
unrealistic double-accounting of default probability. Also notice that 


Vp = 0 > PpaK- eT eT N ey BIT, (22.4) 


where Pz is the breakeven premium for the borrower, in the sense 
that the borrower will find this deal convenient as long as Vg > 0 => 
P > Pg. Compared with Equation (22.2) i.e. the standard DVA with 
liquidity where we observe the double counting 


Vg = 0 > Pg = K- ee Bl ey Br 
As before, assuming 


Pp = K - e B)T 


we may conclude that, in this case, the computation from the standard 
CVA is correct, because it is taking into account the probability of 
default in the valuation of the funding benefit, which removes any 
liquidity advantage for the borrower. 

Equation (22.4) shows what happens when, in addition, there is a 
pure liquidity basis component in the funding cost. On the other 
hand, charging liquidity costs with an adjusted funding spread cannot 
be naturally extended to the case where we want to observe explicitly 
the possibility of default events in our derivatives. 

In writing the payoff for the borrower we have not explicitly con- 
sidered the case in which the deal is interrupted by the default of the 
lender, since we can replace the deal with a new counterparty. This 
keeps Vg independent of the default time of the lender, consistently 
with the reality of bond and deposit markets. 
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22.1.4.2 Risky Funding with CVA for the Lender 


If the lender pays P at time 0, he incurs a liquidity cost. In fact, he 
needs to finance (borrow) P until T. At T, the lender will pay back 
the borrowed money with interest, but only if he has not defaulted. 
Otherwise, he gives back nothing, so the outflow for the lender is 


Pid! ease P ere OF 1 ay ny 


while he receives in the deal K - 1{r(B)>r}. The total payoff at T is 
therefore 


Piglet ey DT . lirD>T} + K - lirB>T} 


Taking discounted expectation we get 


Vi =-P. e TDT (LT oy DT 1 Ter KK 


{t(L)>T} + 


=—-P.eYOT 4 + eB RK, 


{t(L)>T} 
The condition that makes the deal convenient for the lender is 


Vr >= 0s P= PL 


PL=K- et eT BT ozy (DT =K. ertr By DT 


where Pz is the break-even premium. 

It is interesting to note that when the lender computes the value of 
the deal and takes into account all future cash flows as they are seen 
from the counterparties, the valuation does not include a charge to 
the borrower for the component x (L). This term, the cost of funding 
would be associated with his own risk of default. This term is cancelled 
by the fact that funding is not given back in case of default. 

In terms of relative valuation of a deal, this fact about the lender 
is exactly symmetric to the fact that for the borrower, the inclusion 
of the DVA eliminates the liquidity advantage associated with 7 (B). In 
terms of managing cash flows, instead, there is an important differ- 
ence between the parties, which is discussed next. For reaching an 
agreement in the market, we need 


Vz > 0,Vg > 0 
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which implies 
Pr >P>PR S&S 
Ket "BI vOl > ps K. eT e BI yer 


=> 
y (B) = y(L). 


This solves the problem, and shows that, if we only want to guaran- 
tee a positive expected return from the deal, the liquidity cost that 
needs to be charged to the counterparty of an uncollateralized deriv- 
ative transaction is just the liquidity basis, rather than the bond spread 
or the CDS spread. 

This is in line with what actually happened during the liquidity crisis 
in 2007-2009, when the bond-CDS basis exploded. Next we also show 
how the picture changes when we look at the possible realized cash 
flows (as opposed to the expected cash flows). 


22.1.4.3 Positive Recovery 


Now, we study what happens if we relax the assumption of zero 
recovery. The discounted payoff for the borrower is now 


et (P oe e” V ey BT -K) < lftB)>T} 


khpe (P eTe OTOT K) deinen 

=> 

gat (P etl g” BIT oy (B)T -K) . (L - Ltr@)<r}) + Rg ` lr@<7)}) > 

=> 

(P gr pyr -e"'K) -(1-(1-Ra) - 1{2(B)>T})) ; 
where the recovery is a fraction, Rg of the present value of the claims 
at the time of default of the borrower, consistent with standard de- 
rivative documentation. Notice that the borrower acts as a borrower 
both in the deal and in the funding leg, since we represented the latter 


as a reduction of the existing funding of the borrower. By taking the 
expectation at time 0 we obtain 


Ve = (1-(1-Rz)Sp(T))- (P . e” BT oY BIT _ eK) 
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where Sg(T) is the survival probability of B. Using 7 (B) = A(B)(1 — Rp), 
we can apply a first order approximation 


|e 7@T ~ (1 —Rz)- (1 - ey) 
Here A(B) is the deterministic default intensity. Then 


Vg x e TBI (P ; eT OT VOT eTe) = Pp. eV BT _ K. ta BT. 


Rather surprisingly, this is the same formula we found when we were 
studying the risky funding with DVA for the borrower! 

Similar arguments apply to the value of the claim for the lender, 
which acts as a lender in the deal and as a borrower in the funding leg. 
For the lender we have 


V; = -P . eT. Lans? ett (BT K 
is recovered as a first-order approximation of 


Vi =-(1-(1— Rp) SUT)) < P -PT eV OT 4 (1 -1 — Rp)Sp(T)) - EK. 


22.1.4.4 Can DVA Become a Funding Benefit? 


One of the most controversial aspects of DVA concerns the con- 
sequences for the accounting of liabilities on the balance sheet. In 
fact, DVA enables the borrower to condition future liabilities on sur- 
vival and create a benefit on default. However, liabilities are already 
reduced by risk of default in the case of bonds when banks use the 
Jair value option according to international accounting standard, and 
when banks mark the bond liabilities at historical cost. 

So, what is the meaning of DVA? Can we really observe a benefit in 
case of our own default? 

In order to answer this question we need to study what happens if 
the borrower pretends to be default-free, thereby ignoring DVA. 

The borrower can perform valuation for accounting purposes us- 
ing an accounting credit spread 2(B) that may be different from the 
market spread and an accounting liquidity basis y(B) possibly differ- 
ent from the market one, albeit with the constraint that their sum s(B) 
must match the market funding spread. In particular, when the party 
pretends to be default free, we have 2(B) = 0 and y(B) = s(B), and 
there are no more indicators for our own default in the payoffs. 
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Assume that the borrower pretends, for accounting purposes, to 
have zero default risk. The premium P paid by the lender gives the 
borrower a reduction of the funding payback at T corresponding to a 
cash flow P- ee" at T, where there is no default indicator because 
the borrower is treating itself as default-free. 

This cash flow must be compared with the payout of the deal at T, 
which is - K, again without indicator, i.e. without DVA. Thus the total 
payoff at T is 


pigt Oe —K, 
By discounting to time 0 we obtain an accounting value Vg such that 
Vg = P i eS PT -K . ert 


which yields an accounting break-even premium Pg for the borrower 
equal to the break-even in Equation (22.4), 


Pg=K-.- eT eT T oy B)T =K. etr (By (BIT 


where now x (B) and y(B) are those provided in the market. So in this 
case too the borrower B sees a funding benefit that actually takes into 
account its own market risk of default 7 (B), plus an additional liquidity 
basis y (B), thereby matching the premium computed by the lender 
that includes the CVA/DVA term. But now this term is accounted for 
as a funding benefit and not as a benefit coming from the reduction of 
future expected liabilities, thanks to the risk of default. 

This shows how the DVA term can be implemented. When a bank 
enters a deal in a borrower position, it is making funding for an amount 
as large as the premium. If this premium is used to reduce existing 
funding which is equally or more expensive, which in our setting 
means buying bonds or avoiding some issuance that would be neces- 
sary otherwise, this provides a real financial benefit that is enjoyed in 
the case of survival by a reduction of the payments at maturity. A bank 
can buy back its own bonds, which is like ‘selling protection on itself’ 
fully funded. When a sale of protection is funded, there is no counter- 
party risk and therefore no limit to whom can sell protection, which is 
different from the case of an unfunded CDS. In fact buying their own 
bonds is a standard and important activity of banks. The reduction is 
given by the difference in P - e’7e" K. 


22 CVA and DVA 617 


If the quantity of outstanding bonds is sufficiently high to allow the 
implementation of such a strategy, we have shown how the DVA term 
can be seen not as a default benefit, but rather a natural component 
of fair value whenever fair value mark-to-market takes into account 
counterparty risk and funding costs. 


22.1.4.5 The Accounting View for the Lender 


The previous results show that the borrower’s valuation does not 
change if he considers himself default free by using an accounting 
credit spread 7(B) = 0 and treating all the funding cost s(B) he sees 
in the market as a pure liquidity spread y(B) = s(B). Do we have a 
similar property also for the lender? Not at all. 

If the lender computes the break-even premium using an accounting 
credit spread z(L) and an accounting liquidity spread y(L) = s(L) - 
x (L) different from those provided by the market, he gets a different 
breakeven premium, because 


Py = Ke BH OT, 


Thus, the breakeven premium and the agreement that will be reached 
in the market depend crucially on y(L). In Fig. 22.1, for a sample 
deal, we show how Pzr varies when, holding s(L) fixed, we vary 
K(L) = y(L)/s(L), which we call the liquidity ratio of the lender. This 


Lender 
Liquidity ratio 


Borrower Liquidity ratio 


Fig. 22.1 Break-even premium for L, P, as a function of the liquidity cost 
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is not the only difference between the situation for the borrower and 
the lender. Notice that the borrower’s net payout at maturity T is given 
by 

PFO aT) =K liT] 
and is non-negative in all states of the world if provided we keep 
P > Pz, although the latter condition was designed only in order to 
guarantee that the expected payout is non-negative. For the lender 
instead the payout at maturity is given by 

=P. eO ODT. 1r) +K- lirt} 
The condition Vz > 0 does not imply the non-negativity of the earlier 
expression. In particular, we can have a negative carry even if we 
assume that both counterparties will survive until maturity. 

Figure 22.1 shows break-even premium for the lender Pz as a func- 
tion of the liquidity cost ratios «(L) and k(B) when s(L) = 005, 
s(B) = 01, T = 20, K = 100, r = 002. The xyplane crosses the z-axis 
at the break-even premium for the borrower Pg. A deal is possible 
only in the blue region. 

If we want to guarantee a non-negative carry at least when nobody 
defaults we need z(L) < z(B). Otherwise, the lender is exposed to 
liquidity shortages and a negative carry even if the deal is convenient 
for him. Liquidity shortages when no one defaults can be excluded by 
imposing for each deal z(L) < z(B), or, with a solution working for 
whatever deal with whatever counterparty, by working as if the lender 
were default-free. Only if the lender pretends for accounting purposes 
to be default-free will the condition for the convenience of the deal 
based on expected cash flows be 


P<K. e THTT — K. gtr B)+y (L+ (LT 


This clearly implies that 


=P. ettr L +y DT s lO>T} +K. lir(B)>T}- 


should be non-negative. On the other hand, the lender’s assumption 
to be default-free makes a market agreement more difficult, since 


K. ert (B)+y (BT <P<K. eta B)+y DLT 
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implies 
y(B) = y(L)+x(L) 


rather than y(B) > y(L). Under this assumption, uncollateralized pay- 
offs should be discounted at the full funding cost also in our simple 
setting. Let us consider a bank X that pretends to be default-free and 
thus works under the assumption «(B) = «(X) = 1 when the bank X 
is a net borrower and «(L) = «(X) = 1 when X is in a lender position. 
When the bank is in the borrower position, we have 


Pg = Px = Keg VO", 


When it is in a lender position with respect to a counterparty the 
breakeven premium will be given by 


Pr, = K . er ts)T = Pp = Py. 


and the discounting at the funding rate r + s(X) is recovered for both 
positive and negative exposures. 


22.1.5 Final Conclusions 


The previous discussion showed how a consistent framework for the 
joint pricing of liquidity costs and counterparty risk can be formulated. 
This was accomplished by explicitly modelling the funding com- 
ponents of a simplified derivative where both counterparties might 
default. We saw how bilateral counterparty risk adjustments (CVA and 
DVA) could be combined with liquidity/funding costs without any 
unrealistic double counting effects. 

We also found that DVA has a meaningful representation in terms of 
funding benefit for the borrower, so that a bank can take into account 
DVA and find an agreement with lenders computing CVA even when it 
neglects its own probability of default. On the other hand, the lender’s 
cost of funding includes a component that is associated with his own 
risk of default, but this component cancels out with his default prob- 
ability, so that only his liquidity spread (or equivalently his bond-CDS 
basis) contributes as a net funding cost to the value of the transaction. 

We also discussed how the situations of the borrower and the lender 
were different; in particular, the lender could have negative carry upon 
no default even if the value of the deal was positive for him. 
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Thus, while the debate appears to be focused on the impact of ac- 
counting choices on the valuation of liabilities, the previous discussion 
illustrated that it is rather on the valuation of assets that such choices 
make a difference. 


23 


Market Models 


23.1 The LIBOR Market Model 


One of the general disadvantages of short rate and HJM models is that 
they focus on unobservable instantaneous interest rates. The so-called 
market models, that was developed in the late 90’s! tries to overcome 
this problem by instead modelling observable market rates such as 
LIBOR and swap rates. These models can be calibrated to the market 
and have gained a widespread acceptance from practitioners. 

The first market models were developed in the HJM framework 
where the dynamics of instantaneous forward rates are used to determ- 
ine the dynamics of zero-coupon bonds. The dynamics of zero coupon 
bond prices were then used to determine the dynamics of LIBOR. Mar- 
ket models are therefore not inconsistent with HJM models discussed 
in Chapter 16. 


23.1.1 Introduction 


The term structure models studied in the previous chapters have 
involved assumptions about the evolution in one or more continu- 
ously compounded interest rates, either the short rate r(t) or, as in 
the Heath-Jarrow-Morton (HJM) framework, the instantaneous forward 
rates f(t, T). One drawback of these models is that they are expressed 


1 See Miltersen, Sandmann and Sondermann (1997), Brace, Gatarek and Musiela (1997), 
Jamshidian (1997) and Musiela and Rutkowski (1997). 
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in terms of interest rates that are not directly observed in the mar- 
ket. Another disadvantage in HJM is that all model parameters must 
be recalibrated at each point in time; there is no mechanism for se- 
quential updating. Also it is difficult to calibrate the model to actively 
traded prices in the market. However, many securities traded in money 
markets (e.g. caps, floors, swaps and swaptions) depend on discretely 
compounded interest rates such as spot LIBOR rates L(t, t + 5), forward 
LIBOR rates L(t, T, T + 5), spot swap rates lswap(t, ô) and forward swap 
rates Lswap(t, T, ô). For the pricing of these securities it seems more 
appropriate to apply models that are based on assumptions about the 
LIBOR rates directly or spot and forward swap rates. 

We use the term market models for models based on assumptions 
about discretely compounded interest rates. Market models take the 
currently observed term structure of interest rates as a given and 
are therefore to be classified as relative pricing or pure no-arbitrage 
models. Consequently, they offer no insights into the determination 
of the current interest rates. LIBOR market models (LMM) are based 
on assumptions about the evolution of forward LIBOR rates. Similarly, 
swap market models are based on assumptions about the evolution 
of forward swap rates. By construction, market models are not suit- 
able for the pricing of futures and options on government bonds and 
similar contracts that do not depend on money market interest rates. 

In the recent literature various market models have been developed, 
but most attention has been given to the so-called lognormal LMM. 
In such models the volatilities of a relevant selection of the for- 
ward LIBOR rates are assumed to be proportional to the level of 
the forward rate so that the distribution of the future forward LIBOR 
rates is lognormal. Lognormally distributed continuously compounded 
interest rates have unpleasant consequences, but Sandmann and Son- 
dermann (1997) showed that models with lognormally distributed, 
discretely compounded rates are not subject to the same problems. 
Later, we will demonstrate that a lognormal assumption on the dis- 
tribution of forward LIBOR rates implies that pricing formulas for 
caps and floors identical to Black’s pricing formulas can be derived. 
Hence, the lognormal market models provide some support for the 
widespread use of Black’s formula for fixed income securities. 

We have to be aware that lognormal models can’t handle negative 
rates. So we here suppose we have a market situation with strictly 
positive interest rates. 
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23.1.2 General LIBOR Market Models 


In this section we will introduce a general LMM, also referred to as 
the BGM/J model (Brace, Gatarek, Musiela and Jamshidian), describe 
some of the model’s basic properties, and discuss how derivative 
securities can be priced within the framework of the model. 


23.1.2.1 Model Description 


As was described in Section 4.1.14, a cap is a contract that protects 
a floating rate borrower against paying an interest rate higher than 
some given rate K, the so-called cap rate. We let T1, ..., Tn denote the 
payment dates and assume that 7; — 7j_; = ô for all i. In addition we 
define To = Tı — ô. At each time T;(i = 1,...,”) the cap gives a payoff 
of 


Ci(T;) = Nê max{L(T;, T;—-5)—K, 0} = Nô max{L(T;—6, T;—6, T;)—K, 0}, 


where N is the face value of the cap. A cap can be considered as a 
portfolio of caplets, namely one caplet for each payment date with 
payoffs described by the previous formula. 

The definition of the forward martingale measure in Chapter 17 
implies that the value of the aformentioned payoff can be found as 
the product of the expected payoff computed under the 7;-forward 
martingale measure and p(t, T;) the current discount factor for time T; 
payments, that is, 


Ci) = Nép(t, T)E2 [max{L(T;- 8, T;-8, T) -K, 0}]; t < T;-6. 


The price of a cap can therefore be determined as the sum of the value 
of the caplets 


n g'i 
C(t) =N8)_ plt, TEP [max{LT;-ô, Ti- ô, T)-K, 0}]; t< To. 
i=1 


For t > To the first-coming payment of the cap is known so that its 
present value is obtained by multiplication by the risk-less discount 
factor, while the remaining payoffs are valued as the previous one. 
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The price of the corresponding floor is 


n 
Ti 
F() = Nb) plt, TEP [max{K-LT;- ô, Ti- ô, Ti), 0}]; t< To. 
i=1 


In order to compute the cap price, we need to know of the distri- 
bution of L(Ti — ô, T; — ô, Ti) under the T;-forward martingale measure 
Q" for each i = 1,...,n. For this purpose it is natural to model the 
evolution of L(t, T; — ô, T;) under Q” . The following argument shows 
that under the Q" probability measure the drift rate of L(t, T; — ô, T;) is 
zero, that is, L(t, T;—ô, Ti) isa Q" martingale. The simple compounded 
forward rate at time t spanning the future period [T1, T2], L(t, T1, T2) is 
defined by 


pit, Tr) _ 1 


We rewrite this as 


1 t, T;-—6 
LG, t-A T= e a i 
ô \ p(t, Ti) 


where ô = T - Tı. The following diagram illustrates a set of forward 
rates spanning the set of dates T; 


Discount Bond as Numeraire 


Under the T;-forward martingale measure Q" the ratio between the 
price of any asset and the numeraire, that is, the zero-coupon bond 
price p(t, T;) is a martingale. In particular, the ratio p(t, Ti — ô)/p(t, Ti) 
is a Q” martingale so its expected change over any time interval is 
equal to zero under the Q” measure. From the previous formula it 
follows that the expected change (over any time interval) in the peri- 
odically compounded forward rate L(t, T;ô, Ti) also is zero under Q" . 
We summarize the result in the following theorem 
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Theorem 23.1. The forward rate L(t, T; — 6, Tj) is a Q" martingale. 


Consequently, a LMM is fully specified by the number of factors (i.e. 
the number of standard Brownian motions) that influence the forward 
rates and the forward rate volatility functions. For simplicity, we focus 
on the one-factor models 


dL(t, T;—ô, Ti) = B (t, T;— ô, Ti, L(t, Tj, Dlr >) dz(t, Ti), i= 1 ree A, 


where z(t, Ti) is a one-dimensional standard Brownian motion under 
the 7)-forward martingale measure Q. The fourth argument in the 
volatility function 6 indicates that the volatility at time t can depend on 
the current values of all the modelled forward rates. In the lognormal 
LMMs we will study later, one assumes that volatility of each forward 
rate is proportional to the current level of the same forward rate 


B (t, Ti- ô, Ti, L(t, Tj, 5)\r>1) = yt, Ti- 8, T)LG, Ti - 4, Ti) 


for some deterministic function y. However, for now we continue to 
discuss the more general specification. We see from the general cap 
pricing formula that the cap price also depends on the current dis- 
count factors p(t, T1), p(t, T2),...,p(t, Tn). These discount factors can 
be determined by p(t, To), and the current values of the modelled for- 
ward rates (i.e. L(t, To, T1), L(t, T1, T2)),..., L(t, Ty_-1, T,),. Similar to the 
HJM model, the LMMs take the currently observable values of these 
rates as given. 


23.1.2.2 The Dynamics of All Forward Rates Under the 
Same Probability Measure 


The basic specification of the LMM involves n different forward mar- 
tingale measures. In order to better understand the model and to 
simplify the computation of some derivative prices we will describe 
the evolution of the relevant forward rates under a common probabil- 
ity measure. As discussed below, Monte Carlo simulation is often used 
to compute prices of certain derivatives in LMMs. It is much simpler to 
simulate the evolution of the forward rates under a common probabil- 
ity measure than to simulate the evolution of each forward rate under 
a different martingale measure associated with the respective forward 
rate. One possibility is to choose one of the n different forward mar- 
tingale measures used in the assumption of the model. Note that the 
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T;-forward martingale measure only makes sense up to time 7;. There- 
fore, it is appropriate to use the forward martingale measure associated 
with the last payment date (i.e. the 7,,-forward martingale measure 
QT”), since this measure applies to the entire relevant time period. 
In this context Q7” is sometimes referred to as the terminal meas- 
ure. Another obvious candidate for the common probability measure 
is the spot martingale measure. Let us look at these two alternatives 
in more detail. 


23.1.2.3 The Terminal Probability Measure 


We wish to describe the evolution of all modelled forward rates un- 
der a common probability measure - here the 7,,-forward martingale 
measure. For that purpose we shall apply the following theorem 
that outlines how to shift between the different forward martingale 


measures of the LMM. 
Theorem 23.2. Assume that the evolution in the LIBOR forward 
rates L(t, Ti- ô, Ti) fori =1,...,n, where T; = Ti- + ô, is given by 


dL(t,T; -8,T;) = P (t, Ti— 5, Ti, L(t, Tj, Ò)r;>1) det, T), i=1...,n. 
Then the processes z(Ti — 5) and z(T;) are related as follows: 


14+ d6L¢, Ti— ô, Ti) 


dt. 


dz(t, Ti) = dz(t, Ti — ô) + 
Using this repeatedly, we get that 


wL 56 (t, Tj, Tii, L(t, Tj, ÒlT;>1) 


dz(t, Tn) = dz(t, Ti) + 2 1+ ôL(t, Tj, Tj) 7 


Consequently, for each i = 1,...,n, we can write the dynamics of 
L(t, T; — ô, T;) under the Q7"-measure as 


= 6B (t, T;, Tj+1, L(t, Tj, 5)|7,>1) 


dz(t, Tn) - 
ere 3 1+ SL, T; Tai) 


dt |, 
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E y 5B (t, Ti- ô, T;, L(t, Tk, 5)|7,>1) B (t, T}, Tj+1, L(t, Tk, 5)|7,>1) 


= dt 


j=i 
+ B(t, Ti—ô, Ti, L(t, Tr, ©lT >t) dz(t, Ty). 


Note that the drift may involve some or all of the other mod- 
elled forward rates. Therefore, the vector of all the forward rates 
(L(t, To, T1), ..., L(t, Ta-1, Th)) will follow an n-dimensional diffusion 
process so that an LMM can be represented as an n-factor diffusion 
model. Security prices are hence solutions to a partial differential 
equation (PDE), but in typical applications the dimension n (.e. the 
number of forward rates) is so big that neither explicit nor numer- 
ical solution of the PDE is feasible. For example, in order to price 
caps, floors, and swaptions that depend on 3-month interest rates and 
have maturities of up to 10 years, one must model 40 forward rates so 
that the model is a 40-factor diffusion model! However, Andersen and 
Andreasen (2000) introduce a trick that may reduce the computational 
complexity considerably. 

Next, let us consider an asset with a single payoff at some point 
in time T € [7 ,7,]. The payoff N(T) may in general depend on the 
value of all the modelled forward rates at and before time T. Let V(t) 
denote the value of this asset at time ¢ (measured in monetary units, 
e.g. dollars). From the definition of the 7,,-forward martingale measure 
QT it follows that 

F| . 


In particular, if T is one of the time points of the tenor structure, say 
T = Tk, we get 
r| | 


N(T) 


V(t) = p(t, Tr EL” ee ae 
(t) = p(t, Tn) | T) 


N(Tx) 
p(Tk, Tn) 


V(t) = plt, Tr) E2" | 


From 


1 t, T;=ô 
L(t, T;—ô, Ti) = E au) 1). 
ô \ p(t, Ti) 
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we have that 


1  p(Tk, Tk) PT, Tka) Pk, Tna) 
P(Tks Tn) — Pes Th+1) P(Tks Tra2) °° p(Tks Tn) 
= [L + dL(T,, Tk, Tis) + ôL(Tk, Teri, Tk+2)] --. [1 + S(T, Tn-1, Tn)] 
n-1 
= | [E +6L%, Tj, Tal 
j=k 


so that the price can be rewritten as 


n-1 


VQ) = p(t, TES” | NT) [ [L + 8LT Tj, Tix) Fi 
jk 


The right-hand side may be approximated using Monte Carlo simula- 
tions in which the evolution of the forward rates under Q7” is used, as 
outlined earlier. If the security matures at time T,, the price expression 
is simpler 


V(t) = plt, TEL” INT pF. 


In that case it suffices to simulate the evolution of the forward rates 
that determine the payoff of the security. 


The Spot LIBOR Martingale Measure 


The spot martingale measure Q, which we defined before, is associated 
with the use of a bank account earning the continuously compounded 
short rate as the numeraire. However, the LMM does not at all involve 
the short rate so the traditional spot martingale measure does not make 
sense in this context. The LIBOR market counterpart is a roll over 
strategy in the shortest zero-coupon bonds. To be more precise, the 
strategy is initiated at time Tọ by an investment of one dollar in the 
zero-coupon bond maturing at time 7), which allows for the purchase 
of 1/p(To, T1) units of the bond. At time Tı the payoff of 1/p(To, Tı) 
dollars is invested in the zero-coupon bond maturing at time T>, etc. 
Let us define 


I(t) = minți € {1,2,...,n} :T; > t} 


so that Ty denotes the next payment date after time t. In particular, 
I(T;) = i so that Tyr) = T;. At any time t > To the strategy consists of 
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holding 
E l l 
P(To, T) pTi, T2) plio, Tio) 


units of the zero-coupon bond maturing at time Tr. The value of this 
position is 


M(t) 


(t)-1 


A*(t) = p(t, Tr) MO = pt Tro) | | 
j=0 


1 
p(T;, Tj+1) 
I(t)-1 


= ptt, Tio) | | (1+ 8L(7;, Tj, Tj), 
j=0 


where the last equality follows from 


ie =8 T= (? ah 1). 


Since A*(f) is positive, it is a valid numeraire. The corresponding mar- 
tingale measure is called the spot LIBOR martingale measure and 
is denoted by Q* . 

Let us look at a security with a single payment at a time T € [To, Ty]. 
The payoff N(T) may depend on the values of all the modelled forward 
rates at and before time T. Let us by V(t) denote the dollar value of 
this asset at time ¢. From the definition of the spot LIBOR martingale 
measure Q* it follows that 


VO _ po" E fi 


A*(t) _ A*(T)|" | 
and hence 
TA 
-F2 
Vit) =E Fame 


From the calculation 


I(t)-1 
[I [1+ 6L(7;, 7;, T1) 
A*(t) _ p(t, Ty) j0 ae 
ANT): p(T, tian 
I [1 + 6L(7j, Tj, Tj+1)] 
J= 


630 J.R.M. Roman 


I(T)-1 
[| O+ 8T, Ti, Tor 
Jal) 


_ ptt, Tio) 
p(T, Tin) 


we get that the price can be rewritten as 


I(T)-1 


[|] [+ 8L; Ti, TaD 
JIO 


N(T) 


V(t) = plt, Tq )EZ | ————~ 
"i p(T, Tin) 


Fi 


In particular, if T is one of the dates in the tenor structure, say T = Tx, 
we get 


I(T)-1 
V(t) = pt, TE” | NT) || + 8L0j, Tj, TdT! 
JIO 


F, 


since I(Tx) = k and p(Tk, Tiero) = p(Tk, Tk) = 1. 

In order to compute (typically by simulation) the expected value 
on the right-hand side, we need to know the evolution of the forward 
rates L(t, Tj, Tj+1) under the spot LIBOR martingale measure O* . It can 
be shown that the process z* defined by 


dz*(t) =d" H- [o-a] 


is a standard Brownian motion under the probability measure Q*. As 
usual, o7 (t) denotes the volatility of the zero-coupon bond maturing 
at time T. Repeated use of the previous theorem yields 


i-1 
ôB (t, T;, T1, L(t, Tk, ô 
oO) -oA = ‘D 88 (t, Ti, Tint, LG Te Dlr) 


so that 


i-l 
ôB (t, Ti, Tai, L(t, Tk, OIT, 
dz*(t) = dat, T)- > B (t Tj, Tit, LO, Te, Òr) PA 


j=KÀ 1 + ôL(t, Tj, Tj+1) 


Substituting this relation into 


dL(t, T; - 8, T) = B (t, Ti— ô, Ti, L(t, Tj, 5i7>1) dat, T), i= 1...,n. 
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we can rewrite the dynamics of the forward rates under the spot 
LIBOR martingale measure as 


dL(t, Tj -— 8, T) = P (t, Ti— ô, Ti, L(t, Tr, &lr >) dz(t, Ti) 


= B (t, Ti— ô, Ti; L(t, Tk, 8)|r,>1) 
i-l 
OB (t, Ti, Tia, L(t, Ty, VT, 
x dz*(t) > B ( ptt ( k r>) 


dt |, 


HIO 


3 5B (t, T;— ô, Tj, Lit, Tk, Dlr >) B (t, T;, Tj+1, L(t, Tk, 5)|r,>1) 


dt 


JIO 
+£ (t, Ti—ô, Ti, L(t, Tr, 5)|7,>1) az*(. 


Note that the drift in the forward rates under the spot LIBOR martin- 
gale measure follows from the specification of the volatility function 6 
and the current forward rates. The relation between the drift and the 
volatility is the market model counterpart to the drift restriction of the 
HJM models. 


23.1.2.4 Consistent Pricing 


As indicated earlier, the model can be used for the pricing of all se- 
curities that only have payment dates in the set {7}, T2,..., Tn}, and 
where the size of the payment only depends on the modelled forward 
rates and no other random variables. This is true for caps and floors 
on 6-period interest rates of different maturities where the price can 
be computed as before. The model can also be used for the pricing of 
swaptions that expire on one of the dates To, 7),..., Tn-1, and where 
the underlying swap has payment dates in the set {T;,..., Tn} and is 
based on the 6-period interest rate. For European swaptions the price 
can be written as 

\(T)-1 

V(t) = pt, TOE | NT) [| (+ 8L, Tj, TDF 

JIO 
For Bermuda swaptions that can be exercised at a subset of the swap 
payment dates {T1,..., Tn}, one must maximize the right-hand side 
over all feasible exercise strategies. See Andersen (2000) for details 
and a description of a relatively simple Monte Carlo-based method for 
the approximation of Bermuda swaption prices. 
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The LMM is built on assumptions about the forward rates over the 
time intervals [To, T1], [T1, T2], . . -, [Tn-1, Ta]. However, these forward 
rates determine the forward rates for periods that are obtained by con- 
necting succeeding intervals. For example, the forward rate over the 
period [To, 72] is uniquely determined by the forward rates for the 
periods [Tọ, Tı] and [T1, T2] since 


p(t, To) 1) 8 (£ To) p(t, T1) 1) 
p(t, T2) T2- To \p(t, Ti) p(t, T2) ' 


= sill + 6L(t, To, TILL + 6L0, Ti, T2)]- 1), 


1 
L(t, To, T2) = 
(t, To, T2) zrl 


where 6 = Tı -To = T2 - T; as usual. Therefore, the distributions of the 
forward rates L(t, To, T1) and L(t, T1, T2) implied by the LMM determ- 
ine the distribution of the forward rate L(t, To, T2). A LMM based on 
3-month interest rates can hence also be used for the pricing of con- 
tracts that depend on 6-month interest rates, as long as the payment 
dates for these contracts are in the set {7o, T1, ..., Tn}. More generally, 
in the construction of a model, one is only allowed to make exogenous 
assumptions about the evolution of forward rates for non-overlapping 
periods. 


23.1.3 The Lognormal LIBOR Market Model 
23.1.3.1 Black’s Model 


The standard model for valuing OTC interest rate options, caps, floors 
and European swaptions is the Black model. The Black model is used 
by traders in the market to price these derivatives, and as will be 
seen later on, the analytical Black formulas will play a key role when 
calibrating the LMM. 

The basic assumptions under the Black model are the following: 


e The underlying forward rate or swap rate is a lognormally distrib- 
uted stochastic variable. 


e The volatility of the underlying is constant. 
e Prices are arbitrage free. 


e There is continuous trading in all instruments. 
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In Black’s world we denote the forward/futures price with expiry T 
on an underlying with expiry T* as ®(T, T*). The price is lognormally 
distributed with standard deviation o yT — t. It is further assumed that 


on expiry, the expected futures price is equal to the current futures 
price 


ELIOT, T*)|F;] = D, T*). 
For a European call option in Black’s model we have 
Ci =e" [@(t, T*) -N (di [t, O(t, T*)]) -K -N (a [t, OC, T*)])}, 
where 


q, - MOOG TYK) - (07/2)(T-t) _ 


2 di —ovT -t. 


ONT -t 


23.1.3.2 Bond Options 
The unique no-arbitrage value at time t of a forward contract with de- 


livery at time T of a zero-coupon bond maturing at time S at delivery 
price K is given by 


VG, T, $) = pt, S)-K - pt, T). 
The unique no-arbitrage forward price on the zero-coupon bond is 


p(t, S) 


F(t, T, S) = . 
“ j p(t, T) 


Next, we consider a forward contract on a coupon bond where we 
assumed to yield payments at Ti < T2 < ... < Tn in time where 
T < T,. We denote the coupon payments as c;,i = 1,2,...,n. At time t 
the value of the bond is therefore given by 


PO = $ cip(t, Ti). 


Ti>t 
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Let V?"(t, T) denote the time ż value of this forward contract. Then we 
have a no-arbitrage value if the forward given by 


VP", T) = Dept, T)-K- pt, T)=PQ- $, ciplt, T) -K - pt, T). 


Tj>t t<T;<T 


The no-arbitrage forward price is given by 


> r>rcip(t, Ti) _fO= Vrer<replt, Ti) 
= eral 1): 


T;>T 


F?”(t, T) = 


Consider between time ¢ and delivery time T, the two portfolios. 


1. A forward contract, K zero-coupon bonds maturing at T and for 
each T; with t < T; < T, ci zero-coupon bonds maturing at T; 


2. The underlying coupon bond. 
These portfolios have exactly the same payments. At time T, the first 


portfolio equals P(T) — K + K = P(T), which is identical to the value of 
the second portfolio. Therefore the absence of arbitrage implies 


VPG, T)+K -p T)+ Ý cipt, Ti) = PO. 
t<T;<T 


The expected payoff of the forward contract is then given by 


E? [max(p(T) — K, 0) | Fi] = F?” (t, T) -N (di [t, F?"t, T]) 
-K -N (dz [t, F?"(t, T)]) 


where 


In(F?"(t, T/K) — (0 2/2)(T -t 
gam (t, T) E )( leg eal Poa. 
(on —t 


If we multiply the expected payoff with the relevant discount factor 
(e.g. the zero-coupon bond price p(t, T)), we get Black’s formula for 
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a European call option on a bond 


CPG, T, K) = p(T) pa T) -N (d1 [t, F?” (t, T)) | 


=K - N (dz [t, F?” (t, T))) 


=|P©- X cipt, Ti) |N (di [t F"E, T]) 
t<T;<T 
-K - p(t, TN (dz |t, F", T)]). 


Similarly, for a put option we have 


P®”(t, T, K) = K - p(t, T)N (-dz |t, F?"(t, T)]) 


-| P©- Yo cipt, Ti) |N (~a: [t, FP"). 


t<T;<T 


23.1.3.3 Caps and Floors 
For a caplet, with a payoff given by 

7, = Nå max{L(T;, T; - 8) - K, 0} 
we obtain the Black’s price as 


LOT- DN (E Ut, L(t, T;-8, T)I) 


C (t) = Nôp(t, Ti) i -K-N (diit, L(t, Ti— ô, Til) 


where t < T; < 6 and 


ML, T;- 8, TD/K)- (0?2)\(T;-8-®) 
d} = L =d iV T; —6 -t. 
2 a e 


We have assumed that the forward rates F(t, T; — ô, T;) in a risk-free 
world follow the process 


dL(t, T; - ô, T;) = o; L(t, Ti— ô, T;)dV(t). 
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The price of the cap and the floor is given by 


CA = NBD pt, T) 


L T;— 8, TON (di [t, LG, T;— 5, T) 
i=l 


F(t) = N68) pt, T) 


i=l 


K - N (-di [t, L(t, T; - 8, T;)]) 
-L(t, T; — 8, Ti) - N (-di [t, L(t, 7; - 8, TDI) | 


23.1.3.4 LMM Model Description 


The traditional derivation of Black’s formula is based on some inappro- 
priate assumptions. 

First, the assumed lognormality of bond prices and interest rates 
is doubtful. For several reasons the price of a bond cannot follow 
a geometric Brownian motion throughout its life. We know that the 
price converges to the terminal payment of the bond as the maturity 
date approaches. Furthermore, the bond price is limited from above 
by the sum of the future bond payments under the appropriate as- 
sumption that all forward rates are non-negative. When the bond price 
approaches its upper limit or the maturity date approaches, the volatil- 
ity of the bond price has to go to zero. The volatility of the bond price 
will therefore depend on both the level of the price and the time to ma- 
turity. The lognormality assumption can at most be a locally relevant 
approximation to the true distribution. In addition, the forward price 
and the futures price on a bond are not necessarily equal when the 
interest rate uncertainty is taken into account. It is less clear whether 
it is reasonable to assume that future interest rates are lognormally dis- 
tributed, and that the expected changes in the forward rates and the 
forward swap rates are zero in a risk-neutral world. 

Second, the multiplication of the current discount factor and the 
risk-neutral expectation of the payoff do not lead to the correct 
price. In fact, this is true only if we take the expectation under the 
appropriate forward martingale measure instead of the risk-neutral 
measure. 

Third, simultaneous applications of Black’s formula to different de- 
rivative securities are mutually inconsistent. If, for example, we apply 
Black’s formula to the pricing of a European option on zero-coupon 
bond, we must assume that the price of the zero-coupon bond is 
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lognormally distributed. If we also apply Black’s formula for the pri- 
cing of a European option on a coupon bond, we must assume that the 
price of the coupon bond is lognormally distributed. Since the price of 
the coupon bond is a weighted average of the prices of zero-coupon 
bonds and a sum of lognormally distributed random variables is not 
lognormally distributed, these assumptions are inconsistent. Similarly, 
the swap rate is a linear combination of forward rates. When Black’s 
formula is applied for the pricing of caplets, it is implicitly assumed 
that the relevant forward rates are lognormally distributed. Then the 
swap rate will not be lognormally distributed, so that it is inconsistent 
to use Black’s formula for swaptions also. Furthermore, lognormality 
assumptions for both interest rates and bond prices are inconsistent. 

Several research papers suggest other models for bond option pri- 
cing that are also based on specific assumptions on the evolution of 
the price of the underlying bond. The most prominent examples are 
Ball and Torous (1983) and Schaefer and Schwartz (1987). A critical 
analysis of such models can be seen in Rady and Sandmann (1994). A 
problem in applying these models is that the assumptions about the 
price dynamics for different bonds may be inconsistent, and hence 
the option pricing formula obtained in the model will only be valid for 
options on one particular bond. 

To ensure consistent pricing of different fixed income securities 
we must model the evolution of the entire term structure of interest 
rates. In many of the consistent term structure models we shall 
discuss in the following sections, we will obtain relatively simple 
and internally consistent pricing formulas for many of the popular 
fixed income securities. As we shall see in this section, it is in fact 
possible to construct consistent term structure models in which 
Black’s formula is the correct pricing formula for some securities, but, 
even in those models, applications of Black’s formula for different 
classes of securities are inconsistent. 

The lognormal LMM provides a more reasonable framework in 
which the Black cap formula is valid. The model was originally de- 
veloped by Miltersen, Sandmann, and Sondermann (1997), while 
Brace, Gatarek, and Musiela (1997) sorted out some technical details 
and introduced an explicit, but approximate, expression for the prices 
of European swaptions in the lognormal LMM. Whereas Miltersen, 
Sandmann, and Sondermann derive the cap price formula using PDEs, 
we will follow the approach taken by Brace, Gatarek, and Musiela and 
use the forward martingale measure technique, since this simplifies 
the analysis considerably. 
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In the development of Black’s cap pricing formula, we assumed 
among other things that the forward rate L(t, T; — ô, Ti) was a martin- 
gale under the spot martingale measure Q and that the future value 
L(T; — ô, Ti — 6, Ti) was lognormally distributed under Q. However, as 
was seen in the previous theorem this forward rate is a martingale un- 
der the 7;-forward martingale measure Q7 and will therefore not be 
a martingale under the Q-measure. (Remember: an equivalent change 
of measure corresponds to changing the drift rate.) Looking at the 
general cap pricing formula 


n 
C(t) =N8 X plt, TEP" (max{L(Ti) - 8,7; -8,T;)-K,0}] t< To 
i=l 
it is clear that we can obtain a pricing formula of the same form as 
Black’s formula by assuming that L(7; — ô, Ti — ô, Ti) is lognormally dis- 
tributed under the 7;-forward martingale measure Q". This is exactly 
the assumption of the lognormal LMM 


dL(t, T; T ô, Tj) = Lit, Tj =. Ô, Ti)y (t, Ti an ô, Tj) det, Ti), i= 1, 2 eee’ Ma 


where y(t, T;ô, T;) is a bounded, deterministic function. Here we as- 
sume that the relevant forward rates are only affected by one Brownian 
motion, but below we shall briefly consider multi-factor lognormal 
LMMs. 

A familiar application of Itô’s lemma implies that 


1 
d [n L(t, Ti- 8, T)] = -> (7 (t, T; — ô, TD)? dt + y (t, T; — ô, T;) dz(t, T;) 


from which we see that 
T;-6 


1 
MET de= T= MAL TOT / (y (u, T; — ô, Ti)? du 
t 


T;-6 
+ J y (u, T; — ô, T;) dz(u, T;). 
t 


Because y is a deterministic function, it follows that 


Tj-6 Tj-6 
I Y (u, T; ai ô, Ti) dz(u, Ti) m N 0, J (y (u, T; = ô, T)? du 
t 


t 
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under the 7;-forward martingale measure. Hence, 


T,-6 
1 

ln L(t, Tj = ô, Ti) — 2 J (y (u, T; = ô, T)? du, 

In L(T;—ô, T;-6, T) = N Ts f 


J (y (u, T; — 6, TD) du 


so that 7;—64 is lognormally distributed under Q”. The following result 
should not now come as a surprise. 


Theorem 23.3. Under the assumption 
dL(t, T; = ô, Ti) = L(t, T; = ô, Ti)y (t, Ti = ô, Ti) dz(t, Ti), i a l; 2 ee’ Ma, 


the price of the caplet with payment date T; at any time t < Ti — ô is 
given by 


C'(t) = Nép(t, T) {L(t, Ti- 5, Ti) -N (di) -K -N (d)}, 


where 
. W(Lt,T;- 6, TDK 1 
i= (AAT i) ds vi(t, Ti- 6, Ti) 
vt, T; = Ô; T;) 2 
di = di — v(t, T; — ô, T;) 
and 


Ti-6 
v(t, T; — ô, Ti) = | (y (u, T; — ô, T;))? du. 
f 


Note that v,(t, T; — 8, T)? is the variance of In[L(T; — ô, T; — ô, T;)] un- 
der the T'-forward martingale measure given the information available 
at time t. The previous caplet price is identical to Black’s formula if 
we insert 


_ vit, Ti — 6, Ti) 


J/T;-6—-t 


Oi 
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An immediate consequence of the previous theorem is the following 
cap pricing formula in the lognormal one-factor LMM. 


Theorem 23.4. Under the assumption 
dL(t, T; — ô, Tj) = L(t, Ti — 6, T;)y (t, T; — ô, T;) dz(t, Ti), i=1,2...,n. 


the price of a cap at any time t < To is given as 


N 
CA) = Nb X` p(t, T;) {LCt Ti- 6, Ti) -N (di) -K -N (d)}, 
I=1 


where d A and di are as above. 


Fort > To the first upcoming payment of the cap is known and is 

therefore to be discounted with the relevant discount factor, while the 

remaining payments are to be valued as shown previously. 
Analogously, the price of a floor is 


N 
F(t) = N8 Y` p(t, Ti) {K -N (-d)) - LCt, Ti — 8, Ti) -N (~ai) }. 
I=1 


The deterministic function y(t, 7; — ô, Ti) remains to be specified. We 
will discuss this matter below. 

If the term structure is affected by n exogenous standard Brownian 
motions, the assumption on dL(t, Ti — ô, T;) earlier is replaced by 


n 
dL(t,T; - ô, Ti) = L(t, Ti- ô, Ti) X yj (t, Ti— 8, Ti) dzj(t, Ti), 
j=l 


where all y;(t, Ti—ô, T;) are bounded and deterministic functions. Again, 
the cap price is given by the previous formula with the small change 
that v(t, T; — 6, Tj) is to be computed as 


a Tes 


v(t, Ti— 6, Ti) = D | (Ti =8, T}? du. 
j=l ' 
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23.1.3.5 Pricing of Other Securities 


No exact explicit solution for European swaptions has been found in 
the lognormal LIBOR market setting. In particular, Black’s formula for 
swaptions is not correct under the assumption on dL(f, Ti-p—ô, Ti). The 
reason is that when the forward LIBOR rates have volatilities propor- 
tional to their level, the volatility of the forward swap rate will not 
be proportional to the level of the forward swap rate. The swaption 
price can be approximated by a Monte Carlo simulation, which is of- 
ten quite time-consuming. Brace, Gatarek, and Musiela (1997) derived 
the following Black-type approximation to the price of a European 
payer swaption with expiration date To and exercise rate K under the 
lognormal LMM assumptions 
N 
P(t) = Nô X pū, Ti) {L(t,T; -8,T;)-N (d) +K -N (d )} t< To, 
I=1 

where d; and d are quite complicated expressions involving the 
variances and covariance of the time Tọ values of the forward rates 
involved. These variances and covariance are determined by the y- 
function. This approximation delivers the price much faster than 
a Monte Carlo simulation. Brace, Gatarek, and Musiela provide nu- 
merical examples in which the price computed using the previous 
approximation is very close to the correct price (computed using 
Monte Carlo simulations). Of course, a similar approximation also 
applies to the European receiver swaption. 

Under the assumptions of the lognormal LMM Miltersen, Sand- 
mann, and Sondermann (1997) derived an explicit pricing formula for 
European options on zero-coupon bonds, but only for options expir- 
ing at one of the time points Tọ, 7), ..., Ta-1, and where the underlying 
zero-coupon bond matures at the following date in this sequence. In 
other words, the time distance between the maturity of the option 
and the maturity of the underlying zero-coupon bond must be equal 
to ô. Representing the exercise price by K, the pricing formula for a 
European call option is 


C(t, K, T;-6, Ti) = (1 - K) p(t, TN (e) -K [p(t T; - 8) — p(t, T)] N (å) , 
where 

ea 1 ij (1. — K)p(, Ti) 
vi, Ti-8,T)  \ K- [p(t, Ti- 8) - p(t, Ti) 
e e — v(t, T; — 5, T;) 


1 
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and 


TS 
t 


or 


ji T;-6 


v(t, T; -— ô, Tj) = D | Ti -8, 19) du 
j=l t 


The price of the corresponding European put option follows from the 
put-call parity. 


23.1.3.6 Further Remarks 


De Jong, Driessen, and Pelsser (2001) investigated the extent to which 
different lognormal LIBOR and Sap market models can explain em- 
pirical data consisting of forward LIBOR interest rates, forward swap 
rates, and prices of caplets and European swaptions. The observations 
are from the US market in 1995 and 1996. For the lognormal one- 
factor LMM they find that it is empirically more appropriate to use 
a y-function that is exponentially decreasing in the time-to-maturity 
T; — 6 —t of the forward rates, 


y(t, Ti = ô, T;) = ye 


than to use a constant, y(t, Ti — ô, Ti) = y. This is related to the well- 
documented mean reversion of interest rates that makes long-term 
interest rates relatively less volatile than shorter-term interest rates. 
They also calibrate two similar model specifications perfectly to ob- 
served caplet prices, but find that in general the prices of swaptions 
in these models are further from the market prices than are the prices 
in the previous time homogeneous models. In all cases the swaption 
prices computed using one of these lognormal LMM exceed the mar- 
ket prices; that is, the lognormal LMMs overestimate the swaption 
prices. All their specifications of the lognormal one-factor LMM give 
a relatively inaccurate description of market data and are rejected by 
statistical tests. De Jong, Driessen and Pelsser also show that two- 
factor lognormal LMMs are not significantly better than the one-factor 
models and conclude that the lognormality assumption is probably 
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inappropriate. Finally, they present similar results for lognormal swap 
market models and find that these models are even worse than the 
lognormal LMMs when it comes to fitting the data. 


23.1.4 Calibrating the LIBOR Market Model 


In this section we describe how to calibrate LMM to market data. A 
basic assumption is that we have chosen the forward-rate-based LMM. 
This is the natural approach when pricing caps and floors both not 
when pricing swaptions. 

Let the tenor structure be 0 = To < Ti < ... < T1 < Tn and i an 
integer ranging over the reset dates of the rates (e.g. 1 <i< n). 

We define n(t) as the unique index such that Tyg is the next tenor 
date after t. The (one factor) model is given by the following SDE for 
the underlying rates (swap or forward) 


a = Mi (f(t), t) dt + of Ddz(0), 
L 
where 
Si = forward/swap rate at time i 
Hi = drift term 
Oj = volatility of rate i 
z(t) = isa Wiener process 


The solution to this SDE is 
T i T 
F(T) = f;(T) exp I (u (u) — 5) du+ J si (u) dz (u) 
0 0 


The drift terms depend on the choice of numeraire and can be de- 
termined by applying the assumption of no arbitrage. Suppose we 
have forward rates as the underlying rates and choose p(To, Tı) as the 
numeraire. Then the drift terms become 


i 
okfk(Tk+1 — Tk) Pik 
i= 5 + i 


A IAAT - Te)” 
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Determining the time-dependent forward rate volatilities is equivalent 
to calibrating the model. How the calibration is performed is explained 
in a section below. 

Although it is not necessary it will always be assumed here that 
the (instantaneous) volatilities o; of the rates are deterministic (not 
stochastic) functions of time. 

A one-factor model means that all the forward rates are perfectly 
instantaneously correlated. In this case, a single Wiener process is suf- 
ficient to evolve the rates. This is not often a reasonable assumption, 
and eliminates one of the advantages of employing the LMM. An m- 
factor model is one where m-independent Wiener processes are used 
to evolve the rates. In this case the equation becomes 


m 


a = ui fO, N dt+ $ oidd), 1 <k < m. 
' k=l 
This is solved for 
T l Po 
SiT) = fiT) exp / E (u) - 50700) | du + J 2 Oik (u) dz (u) 
0 0 k=1 


The loadings o;g(u) can be interpreted as the sensitivities at time 
u of the ith forward rate to the kth shock provided by the Wiener 
process zg. They must satisfy 


of (t)= > oh. 


k=1 


All that remains before we can start to analyse how the rates will evolve 
is to specify the instantaneous volatilities o;’s and their loadings 0;,x’s. 
This can be done in many different ways. One choice is presented 
below. 


23.1.4.1 Volatility Calibration 


Volatility calibration deals with the determination of the o;’s (the in- 
stantaneous volatility of the forward rate with reset at T;). This is done 
differently for caps and swaptions. Since the cap volatility calibration is 
a first step in the swaption volatility calibration, we start by describing 
the cap volatility calibration. 
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23.1.4.2 Cap Volatility Calibration 


In the Black model for caplets it is assumed that the underlying 
rate has a lognormal distribution with variance equal to T- oF ack 
where T is the reset date of the underlying forward rate. In the 
LMM, this lognormal assumption is also made for each rate separ- 
ately. The instantaneous volatility at reset for each rate is related to 
the aforementioned expression in the following way: 


T; 


J o; (t)dt = Topia 

0 
In other words, the instantaneous volatility at reset for each underlying 
rate is equal to the implied Black volatility, which can be read from the 
market. Although not necessary, we make the assumption that the ops 
are deterministic functions of time only. There are (infinitely) many 
solutions to these equations, and our goal is to pick one that fits our 
needs. We follow the approach suggested by Rebonato (2002). Let 


o(t)=(a+bihe"+d 
and 
di (t) = kio (Ti- t). 


This form is flexible and can by varying the constants a, b, c and d take 
many different shapes. The constants k; are rate specific and are used 
to assure that the caplet prices are exactly recovered. How the k;’s are 
set should become clear below. 


1 Find values on the constants a, b,c and d such that the forward rates 


dfi = ui (f(t), t) dt + 0;(t)dz(t) 


L 


fit as close as possible. We use both the Broyden-Fletcher-Goldfarb- 
Shannon dimensional variable metric method and the Levenberg- 
Marquardt method? in parallel to solve this problem and pick the 
best solution. 


2 The Broyden-Fletcher-Goldfarb-Shannon dimensional variable metric sometimes gives adjust- 
ing factors with the property that the first ones (i.e. i = 1,2,...) have a higher deviation from 
unity than the rest. The Levenberg-Marquardt method yields adjusting factors with a more 
constant deviation from unity. 
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2 Set values on the k;’s by computing 


2 . 
O Black Ti 


Ti 
fo2(@dt 
0 


The second step ensures equality for the forward rates; that is, the 
instantaneous volatility and the implied Black volatility are set to be 
equal at each reset. This completes the volatility calibration for caps. 

Before we end the cap calibration section, we shall discuss the issue 
of deciding the implied Black volatility, for example, OBlack, in the case 
when we are calibrating to an exotic cap with path-dependent strikes. 
For example, consider a ratchet cap where each caplet has a strike 
given by K; = fi-1+X where K; and X denote the strike for the i:th caplet 
and a spread respectively. Recall that opjg-, of a caplet is a function of 
the maturity of the caplet and the strike. 

The fact that opjacx depends on the strike gives us some problems. 
To see this, note that in order to get OBlack of a caplet we must know 
its strike. But if the strike is path-dependent, as it is in a ratchet cap, 
we cannot know the strike beforehand. To solve this problem the 
following approach is taken: 


1. Make good guesses on the start values of the strikes. 


2. Get the Oblack of the caplets by using the strikes from the first step. 

3. Evolve a small sample of the rates (e.g. 1024 Monte-Carlo simula- 
tions). Then compute the average rates for each caplet. 

4. Compute new strikes by using the average rates. 


5. Go to the second step with the newly computed strikes and repeat 
until some desirable convergence criterion is achieved. 


Empirical results show that this scheme always (although not proven) 
converges and gives good estimates on the strikes. Let us now turn to 
the volatility calibration for swaptions. 


23.1.4.3 swaption Volatility Calibration 


We concentrate on the volatility calibration of a Bermudan swaption. 
The volatility calibration of a European swaption is a special case of 
this discussion. The basic set-up of the calibration is that we want to 
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recover the implied Black volatilities of a set of co-terminal swap rates. 
Let T1, 72,..., T, be the expiry dates of the co-terminal swaptions and 
SR; denote the i:th swap rate (e.g. the swap rate for the swaption with 
expiry Ti). Recall that 


SR; = 2 wifi(t), 


j=i 


where 


BaTh- T) 
jt 14 j+ J 
Wj = 7 


X Bra (Tk — Th) 
= 


Furthermore, denote with Oixn(t) the instantaneous volatility of SR;, 
for example, the instantaneous volatility of the swap rate with expiry 
at T; and maturity at T,. It can be shown that 


Six (D? = $ 9 OKOA, 


j=i k=i 
where 


wk Ow OFO 
2° 
(È > wnt) 


m=i 


Gix(t) = 


The last equation is an approximation. One of the main results is that 


n n 
Tix (D? © X Y > Ge(To) pjxoi(Oox()- 

j=i k=i 
This is the key point in the calibration and one ought to under- 
stand how this greatly simplifies our task Git is recommended to read 
Rebonato (2002), “Modern Pricing of Interest-Rate Derivatives”, Prin- 
ceton University Press for a discussion regarding this issue). In order 
to recover the swaption prices, we must have that 


T; 
oBlack (27; = | Oixn(u)? du. 
0 


Our approach to achieving this is the following. 
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Find values of the constants a, b,c and d such that the equations for 
the forward rates fit as close as possible. We use both the Broyden- 
Fletcher-Goldfarb-Shannon dimensional variable metric method and 
the Levenberg-Marquardt method (for a description of these methods, 
refer to Press et. al, (2002) “Numerical Recipes in C”) in parallel to 
solve this problem and pick the best solution. Note that this is the 
same first step as in the cap calibration. 

Consider the last swaption in the set of co-terminal swaptions (e.g. 
with expiry 7;,). This is simply a floating rate exchanged for a fixed rate 
(i.e. a caplet). Set 


o Black AT, 


Onxn 
kn = 


i Onxn(u)*du 
which makes 
o;(t) = kjo (T; -t) 


for i = n an equality. 
Move a step back and consider the swaption with expiry 7,_;. Our 
goal is to set the value for k,_; such that 


Th-1 
gP! k 2 
On- o Tn-1 a On-1xn(u) du. 
0 


Since we have already solved for k, and are using the approximation of 


T; 


2 2 
fo (t)dt = TiO Black 
0 


kn-1 is the only unknown variable. Since k,_; appears in squared form, 
we need to solve a second-degree equation. Although straightforward, 
the algebra becomes quite messy. 

Next we continue to compute the values of the remaining k;’s in a 
similar fashion as in the previous step. Doing this yields values for the 
k;’s such that 


oi (t) = kio (T; - t) 


is fulfilled for 1 < i < n. 
Our approach to calibration of the instantaneous swap volatility can- 
not, as far as we know, be found in any published/academic/scientific 
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article. This is a brief motivation why it seems to be a desirable proced- 
ure. In all but very rare cases, the first step gives values on a, b,c and 
d such that the instantaneous volatilities of the forward rates are time 
homogeneous. A problem might be that the constants a, b,c and d are 
not designed to recover the swaption prices. However, they should 
not be too far off since we are considering the same rates, namely the 
forward rates. Finally, the remaining steps make sure that we recover 
the swaption prices exactly. 


23.1.4.4 Principal Component Analysis 


We use principal component analysis to reduce the number of driving 
factors needed when valuing plain-vanilla caps, European swaptions, 
and Bermudan swaptions. In this section we describe how this is done. 
Note that we do not use principal component analysis for the valu- 
ation of path-dependent caps such as ratchet, sticky, momentum, flexi, 
and chooser. In this case, we use as many factors as there are rates. 
More information on why we do this can be found in Rebonato (2002) 
“Modern Pricing of Interest-Rate Derivatives”, Princeton University 
Press. 

Consider a cap with n caplets with resets at T1, T2, ..., T, respect- 
ively. Each caplet has an associated forward rate f; for 1 <i < n. We 
describe the principal component analysis in a simple case, namely 
when the rates f1, f2,...,fn are evolved to the first reset date Tı. The 
complete picture can then be understood from this discussion. Now 
suppose we have an m-factor model where m < n and let 


Tı 
Cov; = I pijoi(u)o;(u)du 
0 


be the n xn terminal covariance matrix where 1 < i S nand I <j <n. 
Note that Cov;; is symmetric. 

Use the Jacobian transformations of a symmetric matrix method to 
find the n eigenvalues of Cov; and the corresponding normalized ei- 
genvectors as described by Rebonato (2002, section 11.1). Denote the 
vector of eigenvalues with e; = [e] e2 ...e,]/ and the corresponding 
normalized eigenvectors with v;. Furthermore, let vj; = [v1 v2... Vn]. 
Sort e; in decreasing order and change the vectors in vj; accordingly. 
Compute 


Bi = [vevi Verv EE, €mVm | > 
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where | < k < m. Compute 


Cov]] Covn2 Covan T 
Levy | 2 erz deny 
l=1 =1 l=1 


and 


S14/e1V11 SI E2012 . . . $1,/€mVim 
/ S2./@1 021 S24/ €2V22 . . . S24/€mV2m 
ik™ |. 


Sna/ €1Vn1 Snay/ C2Vn2- ++ Sna/ €mUnm 


Finally, compute the model covariance Cov; = =B, x Br. In particular 


note that Cov; = Cov; for all i’s, that is, Ge variance of each rate is 
unchanged, and that B, is the principal component matrix. To put this 
in the context of the section “Libor Market Model”, the equation for 
the drift terms 


i 
okfk(Tk+1 — Tk) pi 
aS kfk(Tk+1 — Tk) pik 


O "E +See Te) 


is rewritten as 


SiTe -T) / 
Lj = —______"_ Cov. 
Deen l+f(Tui-T) © 


and in equation 


a = uif (t), tdt + X oin(Odzx(0), l<k<m 
i k=1 


1 
Oik correspond to B}. 


m 


23.1.4.5 Correlation Matrix 


The n x n correlation matrix 


P11 P12 +++ Pin 
p21 22... 


Pnl ++ Pnn 
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is user defined. After the user has input P;;, we use the heuristic of ap- 
plying principal component analysis to it as described in the previous 
section before the calculations start. In other words, Pj is exposed to 
the same transformation as Cov;; was above. 


23.1.5 Evolving the Forward Rates 


We use two different approaches to evolve the forward rates - the 
short step and the long step. The short step approach evolves the 
“living” forward rates to each reset. To exemplify, suppose that the re- 
sets of the rates are 7), T2,..., Tn. Then fi, fo, ...,f, are evolved to T], 
jo, --+->fn are evolved to T2, and so on. We do not attempt to describe 
the technical details of the approach here. However, the implications 
of the results are that we can evolve a forward rate to an appropriate 
point in time in one step, without hardly any loss of accuracy when 
compared to the short step approach. 

We use both approaches when evolving the forward rates. In par- 
ticular, we use the long step when valuing a cap (of any kind). This 
allows us to evolve each forward rate to its reset date in one step. 
When valuing a swaption (of any kind) we use the short step for all 
but the first time sensitive date (i.e. the first exercise date in the case 
of a Bermudan swaption) - to the first time sensitive date when we 
use the long step. 


23.1.6 Pricing of Bermudan Swaptions 


In this section we make a theoretical elaboration of the pricing proced- 
ure for a Bermudan swaption. This pricing procedure is the only one 
we need to describe; a European swaption is, as we will see, a spe- 
cial case of this discussion, and the pricing procedure of a cap follows 
from its definition. 

A Bermudan swaption contract denoted by X-non-call-Y gives the 
holder the right to enter into a swap at a pre-specified strike rate “K” 
on a number of exercise opportunities. The first exercise opportun- 
ity in this case would be Y years after inception. The swap that can 
be entered into always has the same terminal maturity date, X. A Ber- 
mudan swaption entitling the holder the right to enter into a swap 
in which they pay the fixed rate is referred to as a Payer’s otherwise 
Receiver's. 
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Therefore, as the owner of a Bermudan swaption one is in a position 
where at each exercise date it is necessary to use one’s own judgement 
to determine if exercise is optimal or not. This makes the pricing of 
Bermudan swaptions a little more complicated than that of a European 
swaption. We have implemented Andersen’s strategy that can be found 
in Andersen (1999) “A simple approach to the pricing of Bermudan 
swaptions in the multi-factor LIBOR Market Model”, The Journal of 
Computational Finance, 3(2): 5-32, 1999/2000. Below we describe 
the approach and how we use it. 

When pricing a Bermudan swaptions the most important question 
is how to determine the free exercise boundary. In other words, given 
that the world is in a certain state at one of the exercise dates, under 
what circumstances should one exercise the option? 

Let 


Sse = The European payer’s swaption maturing at time T, and with 
a last cash flow at date Ty. 


Ssx,e = The Bermudan swaption with lockout date (first exercise 
opportunity) 7;, last exercise date T, and final swap maturity Te. 


The decision whether or not to exercise a Bermudan swaption at a 
date T; will in general depend on the state of all forward rates F7(7;). 
To simplify matters somewhat one could make the assumption that the 
strategy depends only on the intrinsic values of the underlying swap. 
Let /(7;), be the indicator function that equals one if exercise is optimal 
at dates T; and zero otherwise. It is hence assumed that 


I(T) =f (Sie, H (Ti), 


wheref is a specified Boolean function with a possibly time-dependent 
parameter H. The relationship 


1 S(T) >0 
0 otherwise 


I(T) = | 


must of course be fulfilled. For the other exercise opportunities the 
following form of I(T;) is assumed: 


1 Sie (Ti) >H (Ti) 
0 otherwise ; 


11) =| 


In this strategy the option is exercised if the intrinsic value of the 
underlying swap is above the barrier H. 
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The next step is to determine the value of the function H at the dates 
T;,..., Tx-H(T). The function H is characterized as being the function 
that maximizes the value of the Bermudan swaption. 

A brute-force way of determining the values of the function H is 
hence to solve the multi-dimensional optimization problem. This could 
be done as follows. Store a Monte Carlo simulation in memory. Simul- 
taneously find values on the function H at the dates 7;,..., 7, such 
that value of the Bermudan swaption is maximized. This would be an 
iterative process where in each iteration, first, new values on H at the 
dates T;,..., 7 are set, then, the value of the Bermudan swaption is 
calculated. Whichever choices on H at the dates 7,,..., 7), that give 
the highest value on the Bermudan swaption are picked. Obviously, 
this way of finding H is tremendously slow and not applicable in prac- 
tice. Fortunately, there is an another way of finding H, as proposed by 
Andersen (1999), which is much more efficient. It can be described as 
follows: 


1. Set H(T,) = 0. 


2. Compute an appropriate Monte Carlo simulation of the forward 
rates and store it in memory. 


3. Consider the Bermudan swaption Sn-1,x,e. At time T;,_; the exercise 
strategy must be the same as for Ssxe since at this time there is 
no difference between the options. H(T,) = 0 (ordinary European 
option) is known and determining the value of H(T;,_1) is hence a 
one-dimensional optimization problem. We solve this optimization 
problem with the Golden Section Search in One Dimension (see 
section 10.1 of Press et. al, (2002) “Numerical Recipes in C” for a 
description of this algorithm). 


4. Repeat in turn the previous step for Sy-2, xe, Sy-3,x,e. +--+» Ss+l, xes Ss,x,e- 


When doing this, it is sufficient to store a single Monte Carlo session in 
memory and to re-use it over and over again. Having determined the 
exercise boundary in this way another Monte Carlo simulation is run 
to calculate the price. This Monte Carlo approach to determine the 
free exercise boundary produces a lower bound on Bermudan swap- 
tion prices that can be shown to be very tight for many realistic term 
structures. 

Finally, note that valuing a European swaption is a special case of 
valuing a Bermudan swaption; a European swaption is a Bermudan 
swaption with one exercise event, Ty, and it follows that A(T;) = 0. 


24 


A Model for Exotic Instruments 


24.1 Managing Exotics 


In the following section, I refer to articles by Patrick Hagan. The adjus- 
tors described here are called “Hagan adjustors”. This is a method of 
turning bad prices into good prices, 

We will study the need of pricing and trading an exotic derivative, 
but because of limitations in our pricing systems, we cannot calibrate 
on the “natural set” of hedging instruments. Instead, we have to cal- 
ibrate on some other set of vanilla instruments, which provide only a 
poor cash-flow replication of the exotic. Consequently, our prices are 
questionable, and if we are bold enough to trade on these prices, our 
hedges will be unstable, chewing up any profit as bid-ask spread. Here 
we discuss how to get out of these jams by using “adjusters”, a tech- 
nique for re-expressing the Vega risks of an exotic derivative in terms 
of its “natural hedging instruments”. This helps to prevent unstable 
hedges and exotic deal mismanagement, and, as a side benefit, leads 
to significantly better pricing of the exotic. First, let us briefly discuss 
how we get in these jams. 

During normal times, the pricing of fixed income derivatives de- 
pends on two key markets. First we have the swap market (or delta 
market), from where we get the yield curve. swap desks maintain 
current yield curves by continually stripping and re-stripping a set of 
liquid swaps, futures, and deposit rates throughout the day. This curve 
determines all current swap rates, FRA rates, forward swap rates, etc. 
The yield curve also shows how to hedge all interest rate risks by 
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buying and selling the same swaps, futures, and deposit rates used 
in the stripping process. 

Second, we have the plain vanilla option market (or Vega market) 
for European swaptions, caps, and floors. Prices of these options are 
quoted in terms of the volatility o, which is inserted into Black’s 1976 
formula to determine the dollar price of the option. 

European swaptions are defined by three numbers: 


1. the exercise date, 
2. the tenor dength) and 
3. the strike (fixed rate) of the swap received upon exercise. 


Keeping track of this market requires maintaining a volatility cube, 
which contains the volatilities o as a function of the three coordinates. 
However, the vast majority of swaptions are struck ATM, i.e. at strikes 
equalling the current swap rate of the underlying forward swap, so 
desks normally track this market by maintaining a volatility matrix con- 
taining the volatilities of ATM swaptions, and a set of auxiliary “smile” 
matrices showing how much to add/subtract to the volatility for strikes 
50 bps, 100 bps, etc., above or below the current swap rate. 

Alternatively, some swap desks determine the adjustment by using 
a smile model, such as the SABR or Heston models. In any case, desks 
are reasonably confident that they can trade the vanilla instruments at 
the indicated prices. 

Now consider the typical management of an exotic interest rate 
derivative, such as a Bermudan swap or a callable range note. During 
the nightly mark-to-market, the deal will be priced by 


e selecting an interest rate model, such as Hull-White or Black- 
Karasinski, 

e selecting a set of vanilla swaptions and/or caplets as the calibration 
instruments, 


e calibrating the interest rate model so that the model reproduces 
the market prices of these instruments, either exactly or in a least 
squares sense, and 


e using the calibrated model to find the value of the exotic via finite 
difference methods, trees, or Monte Carlo. 
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The exotic’s Vega risks will then be obtained by 


e bumping each volatility in the matrix (or cube) one at a time, 


e re-calibrating the model and re-pricing the exotic derivative for each 
bump and 


e subtracting to obtain the difference in value for the bumped case 
versus the base (market) case. 


This results in a matrix of Vega risks where each cell represents the 
deal’s dollar gain or loss when the volatility of that particular swaption 
changes. These Vega risks are then hedged by buying or selling enough 
of each underlying swaption so that the total Vega risks are zero. Of 
course the desk first adds up the Vega exposure of all deals, and only 
hedges the net exposure. 

Calibration is the only step in this procedure, which incorporates 
information about marketing the volatilities. 

Under the typical nightly procedure the exotic derivative will only 
have Vega risks to the set of vanilla swaptions and/or caplets used 
in calibration. So regardless of the actual nature of the exotic derivat- 
ive, the Vega hedges will be trying to mimic the exotic derivative as a 
linear combination of the calibration instruments. If the calibration in- 
struments are “natural hedging instruments” which are “similar” to the 
exotic, then the hedges probably provide a faithful representation of 
the exotic. If the calibration instruments are dissimilar to the exotic, 
having the wrong expiries, tenors, or strikes, then the Vega hedges 
will probably be a poor representation of the exotic. This often causes 
the hedges to be unstable, which gets expensive as bid-ask spread is 
continually chewed up in re-hedging the exotic. 

For example, consider a cancellable 10-year receiver swap struck 
at 7.50%, where the first call date is in 3 years (LONC3@7.50). Surely 
the natural hedging instruments for this Bermudan are the diagonal 
swaptions: the 3y into 7y struck at 7.50%, the 4y into 6y struck at 
7.50%, ..., and the 9y into 1y struck at 7.50%, since a dynamic combin- 
ation of these instruments should be capable of accurately replicating 
the exotic. Indeed, if we do not calibrate on these swaptions, then our 
calibrated model would not produce the correct market prices of these 
swaptions, and if our prices for the 3y into 7y, the 4y into 6y,..., are 
incorrect, we don’t have a prayer of pricing and hedging the callable 
swap correctly. 


658 J.R.M. Roman 


When feasible, best practice is to use auto-calibration for managing 
exotic books. For each exotic derivative on the books, auto-calibration 
first selects the “natural hedging instruments” of the exotic, usually 
based on some simple scheme of matching the expiries, tenors, and 
effective strikes of the exotic. It then re-calibrates the model to match 
these instruments to their market values, and then values the exotic. 
Auto-calibration then picks the next deal out of the book, selects a 
new set of natural hedging instruments, re-calibrates the model, and 
re-prices the exotic, and so on. 

There are a variety of reasons why auto-calibration may not be feas- 
ible. If one’s interest rate model is too complex, perhaps a several 
factor affair, one may not have the computational resources to allow 
frequent calibration. Or if one’s calibration software is too “fractious,” 
one may not have the patience to calibrate the model very often. In 
such cases one would generally calibrate to all swaptions in the volat- 
ility matrix in a least squares sense and the calibration would only 
include ATM swaptions. Alternatively, an interest rate model may be 
more easily calibrated on some instruments than others. For example, 
a multi-factor Brace-Gatarek-Musiela (BGM) model is much easier to 
calibrate to caplets than to swaptions. Finally, one’s software may not 
be set up to calibrate on the “natural hedging instruments”. 

A callable range note provides an example. Consider a regular (non- 
callable) 10-year range note, which pays a coupon of, say, $1 each day 
Libor sets between 2.50% and 6.00%. Apart from minor date differ- 
ences, the range note is equivalent to being long one digital call at 
2.50% and short a digital call at 6.00% for each day over the next 10 
years. Since digital calls can be written in terms of ordinary calls, a 
range note is very, very close to being a vanilla instrument, and can be 
priced exactly from the swaption volatility matrix (or cube). To price a 
callable range note, one would like to calibrate on the underlying daily 
range notes, for if we don’t price the underlying range notes correctly, 
how could we trust our price for the callable range note? Yet many 
systems are not set up to calibrate on range notes. 


24.1.1 At-The-Money Volatility Matrix 


European swaptions are defined by the time-to-exercise (row), and 
length (column) and fixed rate (strike) of the swap received upon ex- 
ercise. A volatility matrix (as opposed to a volatility cube) contains the 
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Table 24.1 A volatility matrix 


a(in %) 3m ly 2y 3y ee 10y 

1m 5.25 12.25 13.50 14.125 Tae 14.25 
3m 7.55 13.00 14.125 14.375 ee 14.50 
6m 11.44 14.25 14.875 15.00 oe 14.75 
ly 16.20 16.75 16.375 16.125 oe 15.50 
2y 19.25 17.75 17.125 17.00 ee 15.75 
10y 14.00 13.50 13.00 12.50 eee 11.00 


volatilities of ATM swaptions, swaptions whose fixed rates are equal to 
the current forward swap rate of the underlying swap. Linear interpol- 
ation is used for the volatilities in between grid points. The 3m column 
is the caplet column. Such a volalility matrix is shown in Table 24.1. 


24.1.2 Migration of Risk 


We will now describe a method for moving the Vega risk, either all of 
it, or as much as possible, to the natural hedging instruments. Suppose 
we have an exotic derivative v which has h1, h2,...,h, as its natural 
hedging instruments. For example, for the 1ONC3 Bermudan struck at 
7.50%, the natural hedging instruments are just the 3y into 7y swap- 
tion struck at 7.50%, the 4y into 6y at 7.50%, ..., and the 9y into 1y at 
7.50%. Suppose that for “operational reasons”, one could not calibrate 


on hi, h2,..., hm, but instead were forced to calibrate on the swap- 
tions and/or caplets $1, S2,..., Sn. Let these instruments have market 
volatilities 0],02,..., On. 


Then after calibrating the model, all prices obtained from the model 
are functions of these volatilities. So let 


d d 
Vine = VACO, 02,..., On) 


be the value of the exotic derivative v obtained from the model. 
Suppose we use the model to price the natural hedging instruments 
hy, h2,...,hm. Let 

He (01, On.. .,0n) k-1,2,...,m 


be the value of these instruments according to the calibrated model. 
Finally, let 


HY k=1,2,...,m 
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be the market price of the natural hedging instruments. Let us cre- 
ate an imaginary portfolio consisting of the exotic derivative and its 
natural hedging instruments, 


where the amounts bg of the hedging instruments will be selected 
shortly. Using the calibrated model to price this portfolio yields 


m 
O -VA (o, On... On) = a2: (Ol, On On): 
k=1 


According to the calibrated model, this portfolio has the Vega risks 


arn aymi «A 3HP 
= bk—+ 
doj = 49; 32 


to the calibration instruments. 

In the next section we will show how to choose the amounts bg 
so as to eliminate the Vega risks, either completely or as completely 
as possible. For the moment just suppose we have chosen the portfo- 
lio weights by. We add and subtract this portfolio of natural hedging 
instruments to write the exotic derivative v as 


m m 
p= v=o bih + X brh 
k=1 k=1 


We now use the calibrated model to value the instruments in the first 
set of braces, and use the market prices to evaluate the instruments in 
the second set of braces. This yields the adjusted price 


m m 
yd = Jy mod _ y by gon m > bH?” 
k=1 k=1 


m 
=V mod + 3 by (ee _ a ). 
k=1 


This procedure is generally known as “applying an adjuster”. The 
terms inside {} are evaluated using the calibrated model, so they 
only have Vega risk to the volatilities of the calibration instruments 
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01, 02,...,0,.With the weights bg chosen to eliminate these risks as 
nearly as possible, the adjusted price v4 has little or no Vega risk 
to the calibration instruments. Instead, the Vega risks of the adjusted 
price come from the last term, 


m 


5 by nee 
k=1 


which only contains the market prices of the natural hedging instru- 
ments. So, as claimed, the adjuster has moved the Vega risks from the 
calibration instruments to the natural hedging instruments. In fact, to 
hedge these risks one must take the opposite position 


— 3 by hr 
kel 


in the natural hedging instruments of the exotic. For the 10NC3 
Bermudan struck at 7.50%, for example, the resulting hedge is a com- 
bination of the 3y into 7y, the 4y into 6y, ..., and the 9y into 1y 
swaptions, all struck at 7.50%, regardless which set of instruments 
were used originally to calibrate the model. 

The last part of the previous equation gives a different view. It shows 
the adjusted price as being the model price corrected for the differ- 
ence between the market price and the model price of the natural 
hedging instruments. 


24.1.3 Choosing the Portfolio Weights 


We wish to choose the hedging portfolio weights by so as to minimize 
the model’s Vega risks. This is an exercise in linear algebra. Define the 
matrix M and vectors Uj by 


ð Hod əri 
Mk= : Us 


00; 1 90; 


and let b be the vector of positions (b1, b2,... ‘Day so that the Vega 
risks to the calibration instruments are 


U-Mb. 
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There are three cases to consider. First suppose that there are fewer 
hedging instruments than model calibration instruments. One can- 
not expect to eliminate n risks with m < n hedging instruments, so 
one cannot eliminate all the Vega risks in this case. Instead one can 
minimize the sum of squares of the Vega risks: 


b=(M'M)!'MU (ifm <n) 


If there are exactly as many hedging instruments as calibration instru- 
ments, then we can expect to completely eliminate the risk entirely 
by choosing 


b=M!U (ifm=n) 


Finally, if there are as many hedging instruments than model cal- 
ibration instruments, then we can select the smallest hedge which 
completely eliminates the Vega risks to the calibration instruments. 
This yields 


b=M'(mm')'U (ifm>n) 


24.1.3.1 Examples 


Consider once more the cancellable 10-year receiver swap struck at 
7.50%, where the first call date is in 3 years. This derivative is normally 
booked as a straight 10-year swap, with a Bermudan option to enter 
into the opposite swap. Here we just price the Bermudan option, the 
option to enter a payer swaption at 7.50% on any coupon date starting 
on the third anniversary of the deal. For the purposes of this example, 
we assume a flat 5% yield curve, and use the Hull-White model with 
the USD volatility matrix from March 1999. 

Clearly the natural hedging instruments are the 3y into 7y swaption 
struck at 7.50%, the 4y into 6y swaption at 7.50%, ..., and the 9y into 
ly swaption at 7.50%. Suppose we calibrate the Hull-White model to 
these “natural hedging instruments” and then use the calibrated model 
to price the Bermudan. This leads to a price of 


V = 200.18 bps. 


This represents the best price available within the one factor Hull- 
White framework. Suppose we calibrate to the same “diagonal” swap- 
tions as before, but instead of calibrating to swaptions struck at 7.50%, 
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we calibrate to swaptions struck ATM, at 5.00%. This yields a much 
lower price, 


vied — 163.31 bps. 


If we add in the adjustor, we obtain the price 


m 
yd + Sb (Ap HE ) = 163.31 + 39.18 = 202.49 bps 
k=1 


a great improvement. 

Alternatively, suppose we calibrate the Hull-White model to the 
caplets starting at 3 years, at 3.25 years, at 3.5 years, ..., and at 9.75 
years, with all caplets struck at 7.50%. Now we have the correct strike, 
but the wrong tenors. The calibrated model yields the price 


vied — 196.82 bps. 


If we add in the adjustor, we obtain a price of 


m 


ymed 4S by (He -Hpo ) = 196.82 + 9.12 = 199.94 bps, 
k=1 


again a distinct improvement. 


24.1.4 Nothing Is Free 


At first glance, it appears that using an adjuster greatly increases the 
computational load. After all, to determine the adjustment requires 
computing the exotic derivative’s Vega risk avmt/ag; to all calibra- 
tion instruments. These risks are usually found via finite differences, 
so evaluating these risks would seem to require model calibrations 
in n + 1 separate scenarios (base case, and each oj bumped separ- 
ately). However, these Vega risks are needed for hedging purposes, 
and are nearly always computed as part of the nightly batch, even if 
one is not applying an adjustor. So computing the Vega matrix is usu- 
ally free. The computational load does increase modestly, because for 
each natural hedging instrument, one has to calculate the model price 
H™°4,, and its Vega derivatives 0H 4/303. This requires calculating 
the model price of m vanilla instruments n + 1 times. This is the same 
load as calculating the calibration error in each of the n+ 1 scenarios, 
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clearly much, much faster than actually calibrating the model in each 
of the n+ 1 scenarios. 


24.1.5 The SABR Volatility Model 


The SABR model is a stochastic volatility model, which attempts to 
capture the volatility smile in derivatives markets. The name stands for 
Stochastic Alpha, Beta, Rho, referring to the parameters of the model. 


24.1.5.1 Dynamics 


The SABR model describes a single forward F, such as a LIBOR for- 
ward rate, a forward swap rate, or a forward stock price. The volatility 
of the forward F is described by a parameter o. SABR is a dynamic 
model in which both F and o are represented by stochastic state 
variables whose time evolution is given by the following system of 
SDEs 


dF; = oF dw, 
do; = ao,dZ; 


with the prescribed time zero (currently observed) values Fo and oo. 
Here, W, and Z; are two correlated Wiener processes with correla- 
tion coefficient -1 < p < 1. The constant parameters f, œ satisfy the 
conditions 0 < 6 < la > 0. 

The aforementioned dynamics is a stochastic version of the constant 
elasticity of variance (CEV) model with the skewness parameter p: in 
fact, it reduces to the CEV model if œ = 0 The parameter a is often 
referred to as the volvol, and its meaning is that of the log-normal 
volatility of the volatility parameter o. 


24.1.6 Asymptotic Solution 


We consider a European option (say, a call) on the forward price 
F struck at K, which expires T years from now. The value of this 
option is equal to the suitably discounted expected value of the pay- 
off max(Fr — K,0) under the probability distribution of the stochastic 
process for F;. 
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Except for the special cases of 8 = 0 and 6 = 1, no closed form 
expression for this probability distribution is known. The general case 
can be solved approximately by an asymptotic expansion in the para- 
meter £ = Ta*. Under typical market conditions, this parameter is 
small and the approximate solution is actually quite accurate. Also sig- 
nificantly, this solution has a rather simple functional form; it is easy to 
implement in computer code and lends itself well to risk management 
of large portfolios of options in real time. 

It is convenient to express the solution in terms of the implied 
volatility of the option. Namely, we force the SABR model price of 
the option into the form of the Black model valuation formula. Then 
the implied volatility, which is the value of the log-normal volatility 
parameter in Black’s model that forces it to match the SABR price, is 
approximately given by 


In(Folk) | [2v2—vi + Enia ( 20CFmid) \” 
D(s¢) 24 


pyi ooC(Fmia) 2-37 
+ + : 
4 a 24 


Oimpl = & 
p a 


where, for clarity, we have set C(F) = F B The value Finiqd denotes a 
conveniently chosen midpoint between Fo and K (such as the geo- 
metric average ./FoK or the arithmetic average (Fo + K)/2). We have 
also set 


Fo 
_ a dx a 1-6 pl-p 
“lo = Siw E) 
K 
and 
_CFmia) B 
C(F, mid) F mid ? 
_C"Fnia) __ A-2) 
OT CEmd Foa 


mid 
The function D(¢) entering the previous formula is given by 


Atire) 


vo=n( 
L=p 
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Alternatively, one can express the SABR price in terms of Black’s nor- 
mal model. Then the implied normal volatility can be asymptotically 
computed using the following expression: 


r Fo- K 272 - vi ooC(F mid) 
Oimpl = D(c) l+ 24 a 


4 PY1 o0C(Finia) 7 2= 3p? g 
4. g 24 


It is worth noting that the normal SABR implied volatility is generally 
somewhat more accurate than the log-normal implied volatility. 


24.1.7 Conversion Between Log Normal 
and Normal Volatility 


We have seen in Section 4.1.15.2 how to convert between normal and 
log-normal volatility for ATM swaptions. This is when the strike rate K 
and the forward rate F are equal. We will now give a general formula of 
how to convert between volatilities. The subsequent formulas follow 
the articles by Hagan and Woodward (1998) and Hagan! where they 
used a singular perturbation expansion. 

The Black’s log-normal model is 


dF = ogFdW, F(0) =f, 


where f is today’s forward swap/caplet rate and where og is the 
implied Black dog-normal) volatility. The Black’s normal model is 


dF = 0;,dW, F(0) =f, 


where oy is the “normal” or “absolute” or the annualized “basis point” 
volatility. For a swaption with strike (fixed rate) K, the normal volat- 
ility oy (which gives the same price of the option) as the log normal 
volatility og is 


f-k : 
OB . ? 
2 
mK) 14 (1- aby [ngo] ) oft + soit? 


ON = 


1 Hagan, Volatility Conversion Calculators 
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where t is the time to the exercise date in years. When f- > K, 
the aforementioned formula goes to a “O over 0”. To avoid this 
complication, we should use the alternative formula 


1+ 4 [n/K] 
4 
on =0Byf K 1-2 I 42 
1 + 5306T + 376008 T 


when 


aa “ < 0.001. 
To calculate the Black log-normal volatility from the normal volatil- 
ity, we apply the previous formulas in a Newton-Raphson solver (or 
similar). 

Since Black’s log-normal model, in a risk-neutral world, can be 
written as 


dF 

— = ogdW 

p” 
and the normal Black as 

dF = oydW 


we say that og is a relative volatility while oy is an absolute volatility. 


24.1.8 Conversion Between Normal and CEV Volatility 
The CEV model is 
dF = aFfdW, F(0)=f. 


To convert the CEV volatility into a normal (absolute) volatility, one 
can use 


iat =K) d 

Bp _ Kl 2.846? . 

CJP KP 1 AREMO eD 02, 
1-48 A 
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When f is very near K, or when £ is very near 1, the previous formula 
is singular. To avoid this complication, we should use the alternative 
formula 


1 2 

= olf. p2. Habo 
on =a(f - K) EET 
1 24 


9. 2 
1S gI? B-P) 


_B)2 KI 
C8 Tinga? rT? 


1+ at 


when 


(1-8) — < 0.001. 


To calculate the CEV volatility from the normal volatility, we apply the 
aforementioned formulas in a Newton-Raphson solver (or similar). 


25 


Modern Term Structure Theory 


25.1 Term Structure Theory 


From the bootstrapped Swap curve, see Section 6.1.4 we have a given 
yield curve. This curve is used for discounting all cash flows. We define 
the discount function D(t) as 


D(t) = exp - f fo, T^dT' 
0 


To price exotic instruments we use the following schema: 


1. To examine the risk, select an arbitrage free model. 
2. Calibrate the model by 
a. selecting some vanilla hedging instruments 
b. matching the current discount curve D(T) 
c. matching the vanilla prices (caps, floors, swaptions) by e.g. 
least square method. 


3. Price the exotics via the model and interpolation of vanilla instru- 
ments. 


You can use a global calibration process with known instruments or a 
local calibration with the most likely vanilla instruments. 
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25.1.1 The Three Elements 


We use three elements as described next: 


1. A numeraire N(t). Use a positive value of a coupon free security, 
e.g. the money market account 


t 


N(t) = exp J r(t’)dt’ 
0 


or zero-coupon bonds 
N(t) = p(t, T) 


2. The valuation in a risk neutral world is then given by 


E V(T) 
V(t) = NOE Fad 


3. We also use a random evolution of the interest rates. 
There exist three approaches to model the interest rates: 


1. Heath-Jarrow-Morton (HJM) and BGM (Brace Gatarek Musiela) 
sometimes also called the Libor market model (LMM), 


2. A short rate model, such as Ho-Lee, Hull-White CHW), Black- 
Karasinsky (BK) etc. 


3. A Markov model. 


25.1.2 The BGM Model (Brace Gatarek Musiela) 


The BGM model is a financial model of interest rates. It is used 
for pricing interest rate derivatives, especially exotic derivatives like 
Bermudan swaptions, ratchet caps and floors, target redemption 
notes, auto-caps, zero-coupon swaptions, constant maturity swaps and 
spread options, among many others. 

The quantities that are modelled in BGM are a set of forward rates, 
which have the advantage of being directly observable in the mar- 
ket, and whose volatilities are naturally linked to traded contracts. 
Each forward rate is modelled by a log-normal process under its for- 
ward neutral martingale measure, i.e. a Black model leading to a Black 
formula for interest rate caps. This formula is the market standard to 
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quote cap prices in terms of implied volatilities, hence the term “mar- 
ket model”. The LMM may be interpreted as a collection of forward 
LIBOR dynamics for different forward rates with spanning tenors and 
maturities, each forward rate being consistent with a Black interest 
rate caplet formula for its canonical maturity. One can describe the 
dynamics of the different rates under a common pricing measure, for 
example, the forward-neutral measure for a preferred single maturity, 
in which case forward rates will not be log-normal under the unique 
measure in general, leading to a need for numerical methods such as 
Monte Carlo simulations. 

The LMM models a set of n forward rates L;, as log-normal processes 


dLi(t) _ 
L 


iLi}, thdt+o(tdW; i=1,...,n. 


Here, L; denotes the forward rate for the period [T;, Ti+1]. For each 

single forward rate, the model corresponds to the Black model. The 

novelty is that, in contrast to the Black model, the LMM describes the 

dynamic of a whole family of forward rates under a common measure. 
The valuation is given by 


ma Cag K RT) | 


= R(T) = 
Zo(t) Zo(T) 
dR; = y(t, R)dt + a(t, RAW; k=1,2..., 


where V(t, F) will satisfy 


av 
2 DDI = 


The function 


j j i 
ak di ak 
yo ee =} 9 pinaias 
; 1+ akrk - 1 +airi 1+ ars 
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which gives the BGM drift condition 


M(t, t) = Vinee 


Drar 
j=l J 


Compare this to the HJM drift condition. If we let Rm be the interest 
rate in the interval [Tm-1, Tm] and Zm = Z(t, Tm), then we can use the 
numeraire 


Zm(t) = Zo(t) I ET 


We now remember that a change of numeraire only changes the drift. 
We now have the valuation formula 


WE _ 5 R(T)) Ruis 7| 
Zm(t) Zm(T) 


dR; = vit, R)dt + a(t, R)dWe k=1,2... 
If we replace Zp with Zm in the preceding calculation 


m 


v(t, 7) = kajak-—— 
2. D Piaj 1 +ajrj 
j=l j=l 
for the numeraire Z,,(t). If we set k = m, we have no drift, i.e. 


dRm = Ant, Rm)dWin 


i.e. Rn is a martingale when Zm(t) is numeraire. 


25.1.3 A Caplet in the BGM Framework 


In this section we will study a caplet with strike Ry, a tenor a, and a 
FRA fair rate Rm over the interval [Tm-7, Tm]. Its payoff is given by 


payoff = Om(Rm(t) — Re) Zn) t < Tm- < Thm. 
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If we use 
V.(0, R? VAT, Rn (T 
ORI) _ p [VT Rn egy = po 
Zm(0) Zn(T) 
dRm = ax(t, Rm)dWm 
we get 


Ve(0, RG) = &mDTmdE | (Rn) — Rf)” Rn) = Ry] 


m 


If we use Black’s log-normal model: am(t, Rm) = OmRm we get 


Ve = tmin BS (1,89, Rj, Om) - 


m? 


If we use a local volatility model: am(t, Rm) = Am(Rm)Rm we get 


Ve = omD(Tn)BS (t, R9» Rf, om) 


l 0 
ET 5 ( M+R) U+... 


Using a SABR (Stochastic Alpha, Beta, Rbo) model we have 


’ 


Rg = vet, R)dt + wagR dW, k=1,2... 
dæ = yædW 


where the correlation matrix is given by 


dR, = v(t, R)dt + apRpdWe 
dW;W; = pjgdt = Corr fdR;, DRx} dt 
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and the volatilities a;(t, Rx) are determined by the volatility and smile 
of caplet k. The correlation matrix pj then determines swaption volat- 
ilities. High correlations give high swaption prices. The correlation 


matrix can be determined by 


1. historical correlations 


2. using factor analysis to write pj; = Maik + 5q2;G2K+ ... where 


only two or three factors, q1, q2,..., are significant 


3. adjusting factors to match swaption volatility as closely as possible. 
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The strength of the BGM model is that you can understand intuit- 
ively the relationships between caplet volatilities, correlations and the 
swaption volatility matrix. 

The weaknesses are as follows: 


e High dimensions, one per Rg. So, a 10-year deal has for example 40 
rates and 40 dimensions! 


e Limit valuation methods to Monte Carlo. 


e Monte Carlo is very slow, is very noisy, gives bad hedges and is very 
challenging for multiple exercise deals. 


So do not use Monte Carlo when there is an alternative. 


25.1.4 Short Rate Models 


The roll-over numeraire N(t) is the money market given by 


t 


N(t) = exp [rear 
0 


The valuation in a risk neutral world is 


t 


V(t, ro) = E / r(t’)dr- V(T, r(T))|r(t) = ro 
0 


Generalized Cross-Currency Interest Rate (CIR) models 
dr = [0(t)—K(t)- r] dt + o (trf aW, 


where £ is a skew, 6 = 1/2 gives CIR and f = 0 HW. The generalized 
BK models can be written by 


dr = [0(t) —k(t)- Y]dt+o(t)dW, 


where Y is Gaussian. Y = logr gives BK, Y = r the HW and Y = r!’ 
the generalized BK. 
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We can also use a multi-factor generalization. The calibration is made 
on spot prices of zero-coupon bonds 


t 


p(t, Yo, T) = E J r(t’)dt’ - 1|Y(t) = Yo 
0 


With a given x(t) and o(t) adjust 0(t) to match the yield curve; 
D(T) = p(O, Yo; T) or use the Jamshidean’s forward induction scheme 
and adjust «k (t) and o (t) to match caps/floors/swaptions. 

Short rate models are easy to work with. An n factor model has n 
dimensions. But the most serious valuations are done with either one 
or two factors. For n = 1 or n = 2 you can use tree models, lattice 
methods or finite difference methods to value deals with a fast and 
accurate result. The calibration process is often straightforward. 


26 


Pricing Exotic Instruments 


26.1 Practical Pricing of Exotics 


In this chapter we will give an introduction of valuation of exotic in- 
terest rate derivatives in a Gaussian framework and how to calibrate 
such models to market data of plain vanilla instruments. 


26.1.1 Discount Factors, Zeroes and FRAs 


Suppose at date t, one agrees to loan out $1 at date T, and get repaid 
the next day 


1 paid at T, 
1+f(t,T)AT received at T + AT 


By definition, the fair interest rate to charge is f(t, T) = instantaneous 
forward rate for date T as seen at date t. 

Now suppose at date t one agrees to loan out $1 on Ty, with the 
money repaid on Tena. Economically this is equivalent to loaning out 
$1 on Ts, getting repaid $1 plus interest the next day, re-loaning out 
the $1 plus interest, getting repaid $1 plus interest plus interest on the 
interest, .... Clearly, if one agrees at date ¢ to loan out 


1 paid at Tst 
the agreement should specify getting repaid 
Tend ¢ 7 7 
eltg FET MT paid at t 
1 received atT ond 
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for the deal to be fair. Alternatively, we can rephrase this as 


Tend 4 7 
e fia TET aT paid at t 


1 received at Teng. 


This type of single payment deal is equivalent to an FRA (forward rate 
agreement). Suppose we imagine that we are at date t, and we ask how 
much I would need to pay immediately to receive $1 at date T . Clearly 
the fair amount is 


Tena 7 7 . 
e~ hi AGT AT paid atTsr 
1 received at Tend. 
By definition, 
pe, Ty =e Morar 
i.e. the value at ¢ of $1 paid at T is the value of a zero-coupon bond for 
maturity T on date 0. Today is always t = 0 in our notation. Discount 
factors are today’s values of the zero-coupon bonds 
T JTI 
D(T) = pO, T) = ho SOD, 


where f (0, T) is today’s instantaneous forward rate curve. 
26.1.2 Swaps 


We start by studying the fixed leg. Consider a swap with start date fo, 
fixed leg pay dates t1, t2, . . . , tn, and fixed rate RF”. The fixed leg makes 
the payments: 


ajRM* paid at t,i = 1,2,...,n— 1 
l+a,R™ paid at ty, 
where 
aj = cvg(ti-1, ti,B) 


is the coverage (day count fraction) for interval i computed according 
to the appropriate day count basis 6. Discounting the future payments 
with some given market discount factors p(t, ti) the current value of 
the fixed leg on any given day t will be 


n 
Pilt) = R™ 2 aip(t, ti) + P(t, tn). 
= 
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Then we study the floating leg of the swap. Floating legs usually have 
a different frequency than the fixed legs, so let the reset and pay dates 
of the floating leg be denoted by 


to = T0,T1,--->Tm = Ín. 
We suppose that the pay dates for the floating leg are 10,71, ..., Tm- 
These will typically not be the same dates as the pay dates t1, f2,..., tn, 


for the fixed leg. Indeed, the number of pay dates for the fixed and 

floating leg may differ, often m > n. We do however require that the 

first and last pay dates for both legs coincide, i.e. to = To, and Tm = th... 
The floating leg pays 


airj paid at tj, j = 1,2,...,m— 1 
1+, rm paid at Tm = tn, 
where 
a’ , = cvg(ti-, Ti, B’) 


is the coverage for interval j computed according to the appropriate 
day count basis p’. Here r; is generally the Libor or Euribor floating rate 
for interval j. This rate is set on the fixing date; for most floating legs, 
the fixing (reset) date is 2 London business days before the interval 
starts On Tj-1. 

We may consider the floating leg as a sequence of FRAs in which 
one lends 1 unit of the currency at t;_; and receives | + rj units back 
at Tj. On the later date 1 unit of currency is lent again so the net cash 
flow of the series of FRAs will exactly duplicate the cash flow of the 
floating leg, i.e. a;r) for j < m and the final payment 1 + œ/,rm , the 
notional does not relent on date Tm . Taking a closer look at one of 
the FRAs we note that the current value of the lent amount must be 
equal to the current value of the repaid amount on each fixing date 


z “ar Using the relevant market discount factors we have 
PaF, 1) = (1 + afr) Bt, 7). 


From this formula we can solve for the discounted value of the interest 
payment a'r; on each fixing day 


vier”) = pee", Tj-1) prt", Tj). 


If the value of the floating rate payment is the difference between two 
freely tradable securities (two zero-coupon bonds) at the fixing time, 
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then the value must equal this difference for all earlier times as well. 


So in principle, the value of the j” floating interest rate payment is 


vere =p(t, Tj-1) —p(t, Tj) fort ar” 


for any date ¢, at least until the rate is fixed. At any date f, the for- 
ward fair or true rate r” “e (t) is defined so that the value of the interest 
payment exactly equals the theoretical value 


ari” “€(t) ‘p(t, tj) = theoretical value of interest rate payment 
= Plt, Tj-1) — P(t, Tj) 
So 
P(t, Tj-1) — p(t, Tj) 


r O = 


In practice, floating rates are not set at the fair rate; they are set at 
the fair rate plus a small offset s;, the forward basis spread, due to 
credit considerations and supply and demand. The value of the (for- 
ward) basis spread depends on which index is used for the floating 
rate (3M USD Libor, 1M fed funds, 6M Euribor, etc.), and on the start- 
ing date. Taking the spread into consideration the value of the floating 
rate payment paid at qj is 


Vit) = ÊC, t1) - ÊC, G) + a's {P(t 1). 
By definition, the forward rate for the floating rate is defined by 


at’ F Ped cr, Tj) = value of interest ratepayment 
so 
F A(t, Ti] —p(t, T; 
OOL _ Pl j l) P(t, Tj) +sj. 
j j 7 : 
a P(t, Tj) 


26.1.3 Basis Spread 


Basis spread curves are obtained by stripping basis swaps. One can 
show that forward basis spreads are martingales with respect to the 
appropriate forward measures. Since they are very small, usually just 
1-2 bps, and since they seldom vary, one always assumes they are con- 
stant. That is, one assumes that the gamma of the forward spread is 
inconsequential. 
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Summing these payments together, the value of the floating leg is 


m 


Varlt) =C, to) + X a’ spUt, 1). 
j=l 


This is true regardless of the model used. The value of the receiver 
swap (receive the fixed leg, pay the floating leg) is 


m m 


Vrec(t) = R™ X aplt, ti) +P, tn) Plt, to) — X > a splt, 1j). 
i=l j=l 


The value of the payer swap (pay the fixed leg, receive the floating leg) 
is 


m n 


Vpay(t) = Vrec(t) =P, to) + X a splt, tj) —R™ X apt, ti) —BU, tn). 
jel i=l 


26.1.3.1 Handling the Basis Spread 


Basis spreads are a nuisance. They are too big to be neglected (except 
for USD 3m Libor), yet small enough to be nearly irrelevant. One way 
of handling them is to treat them as another, very small, fixed leg of 
the swap. There is really nothing wrong with this approach, although 
one usually has twice as many fixed leg pay dates. 

However, here we will use a second, common approach in which 
each interval’s fixed rate is adjusted to account for the value of the 
basis spreads. If the basis spread is 0.625 bps in an interval, then we 
subtract 0.625 bps from the fixed rate instead of adding it to the float- 
ing leg. More precisely, including this adjustment, today’s value of the 
swap is, 


n m 
Vrec(O) = R™ $` ajDUti) + Dn) — Dito) — È a;s;jDCT)), 
i=1 j=l 


= JU œ (R — Si) Di) + Dt») — D(t). 
i=1 


Here S; is the basis spread expressed with the same frequency and day 
count basis as the fixed leg. If the floating leg frequency is the same or 


682 J.R.M. Roman 


higher than the fixed leg frequency, then 
‘> a’ iP i) 
Sel 
aj D(ti) 


where j € I; represents the floating leg intervals which are part of 


the i” fixed leg interval. That is, the floating leg intervals whose the- 


oretical dates p are contained in the i” fixed leg theoretical interval 


? 


pE < t < th. If the fixed leg frequency is shorter than the floating 
leg frequency (this is rare), then the same S; is used for all fixed leg 
intervals forming part of each floating leg interval. So, 

a’ ;5jD(%) 


EDY 


i=l 


where i € J; represents the fixed leg intervals i with ari < A < r”. 


j 
Either way, we may approximate the swap values as 


A n A rs A 
Vrec(t) = X a; (R™ — Si) p(t, ti) +Ê(C, tn) —B(t, to) 
i=1 


n 
Ppay(®) = P(t, to) ~P(t, tn) - 2a; (RF — Si) ÊC, ti) 

i= 
for all dates t, where the strike RY © and effective spread S; are known 
constants. We are neglecting any evolution of the basis spreads and any 
minor differences due to the differences between the legs’ day count 
bases and frequencies. We will use this approach throughout. Compu- 
tationally, it would be just as easy to modify the code to add another 
fixed leg, but this would make the formulas messier and debugging 
more difficult. 


26.1.3.2 Swap Rate and Level 


At any time 7, the swap rate R™ (t) is defined to be the break even rate, 
the value of R’ ™ which would make the swap value equal to zero. 
Clearly, 


n 
P(t, to) P(t, tn) + X aS P(E, ti) 
i= 


R” (H = 


2 


LO 
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where the level L(t) (also known as the PVO1, the DVO1, or the 
annuity) is 


L(t) = X apt, ti). 


i=] 


We can rewrite the swap values in terms of the swap rate and level as 
Vrec(t) = |R" -ROLO Pat = [RO - RM LO. 


In particular, today’s swap rate and level are 


Do) = D(t) + oS:D(t) 


QO i=] 
R= 70 


n 
=$ aD(ti) 
i=l 


and the swap values are 
Virec(t) = [r = R’| EP Viet) [r° = Rf] L. 


Swaptions. A swaption is a European option on a swap. Consider a 
receiver swaption with notification date rf. If one exercises on this 
date, one obtains the receiver swap. Clearly 

A opt , qt 

ye (1) = [r _ R "| LG) 


réc 


is the value of the receiver swaption on the exercise date. Swaption 
prices are almost always quoted in terms of Black’s model. To intro- 
duce this model, suppose we choose the level L(t) as our numeraire. 
dt is just the sum of a bunch of zero-coupon bonds, and hence is a 
tradable instrument.) The function L(t) is sometimes called the for- 
ward annuity. There exists a probability measure in which the value 
of all tradable instruments CGincluding the swaption) divided by the 
numeraire is a martingale. So 


r Vee) 
opt = 
VP!) = LAE l LT) rd : 


= t~* 


If we evaluate the expected value at T , we see that 


P(t) = LINE [[R® z ea) | F,| 
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Moreover, the swap rate 


R” _ P(t, topl, tn) + 2i = Ina;SÊ(t, ti) 


is clearly a tradable market instrument (a bunch of zero-coupon bonds) 
divided by the numeraire. So the swap rate is also a martingale. By the 
martingale representation model, then, we conclude that 


dR” = A(t, *)dW, 


where dW is Brownian motion, and A(t, x) is some measureable coef- 
ficient. Fundamental theory can take us no further. We now have to 
model A(t, x). Black proposed that A(t, x) = oR”, so that the swap rate 
is log normal 


dR” = oR”aW. 


Finding the expected value previously under this model yields Black’s 
formula 


Vre () = {R™N(d1)— R” ON} LO 


fix 
log Ro Æ jo? (tex —t) 


o V tey —t 


di2 = 


Today’s market price of the swaption is 


amkt 


pnt) = [R Na) 7 Rnd) ip 


R 1 12 
log RO + 79 tex 


dy 
; oN tey 
Do -D +n ajS;D; n 
0 _ 0 n i=] iSi i _ D. 
R? = D Lo = Dist aD. 


Here D; = D(t;) are today’s discount factors. A payer swaption is a 
European option to pay the fixed leg and receiver the floating leg. The 
value of the payer swaption is obtained by reversing RF” and R° in the 
aforementioned formulas. 

If one analyses Black’s formula, one discovers that the receiver and 
payer swaption values are both increasing functions of the volatility o. 
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Instead of quoting swaption prices in terms of dollar values, one can 
just as well quote the price in terms of the value of o that needs to be 
inserted into Black’s formula to obtain the market price. This value of 
the volatility is known as the implied volatility. 


26.1.4 Caplets and Floorlets 


Consider a floorlet for the interval tọ to t1. The floating rate r for the 
interval is set on the fixing date tey two (London) business days before 
the interval starts at to, and the floorlet pays the difference between 
the strike (fixed rate) and the floating rate at the end of its period, 
provided this difference is positive: 


a(Rjix—r)* paidat t7. 


Here a@ is the coverage (day count fraction) of the interval [To, T1]. As 
before, the value of the floating rate payment is 


Pltes To) -P (tex, T1) + OS P(tex, 1); 


on the fixing date, where sı is the (forward) basis spread for the 
interval. The floorlets payoff is 
re 


V foortet(t) = [a (Rf = 51) Blter T1) +P(tex, T1) —P( tex, 0) | 


This is the same payoff as a 1-period receiver swaption. Similarly, caplet 
payoffs are identical to the payoffs of 1-period payer swaptions. 

The analysis of caplets and floorlets parallels the analysis for swap- 
tions exactly. We define the forward or FRA rate as 


n 
P(t, To) —P(t, Tn) + D> ais p(t, Ti) 
Hi 


ap(t, T1) 


and choose the zero-coupon bond p(t, t1) as our numeraire. The value 
of the floorlet is 


n A + 
Pjoorleiltex) - Plte THE | [RM - R^) 1 Fi], 
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where the forward FRA rate is a martingale under this measure. 
Modelling this rate as log normal 


dR" = = o RP Maw 
again yields Black’s formula, 


Vitoorter(t) = {RN (dy) — RFRA (t)N(d2)} aP, t1) 


In (RO /RFRA(t)) + 507 (tex —1) 


O /tex —t 


di2 = 


Today’s market price of the floorlet is 


amkt 


V foorler(0) = {RN (a9) — R°N(dS)} aD(t1) 
In (R/R?) + 507 (tex) 
o tex 
Do -—Dy+as,D, 
aD, ` 


0. _ 
di2 a 


Ro = 


The value of the caplet is obtained by reversing R and R° in the 
previous formulas 


Vcaplet 


Vink (0) = | RONA) - RNN) } aDC) 


= Virie() ~ [R - R°} aDC). 
Note that the caplet and floorlet values are special cases of the payer 
and receiver swaptions with n = 1. As before, the implied volatility 
o is the value of the volatility which makes the previously mentioned 
formulas match the actual market values of the floorlet and caplet. 


26.1.5 Linear Gaussian Models 


A modern interest rate model consists of 3 parts: a numeraire, a set of 
random evolution equations in the forward risk neutral world, and the 
matching martingale pricing formula. The one-factor Linear Gaussian 
Models (LGM) has a single state variable, X. It starts at today = 0 and 
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satisfies 
| dX(t) = a(t)dW(t) 
X(0) = 


This is the evolution under the risk neutral measure induced by the 
numeraire, which will be named shortly. Clearly X(t) is Gaussian with 
the transition density 


y(t,x|T,X)dX = prob [x <X(T) < X+dXK() = x| 


given by 
1 2 
(t,x|T, X) = ————e -4 (Xx)? /AE 
Nl ree 
1 iz 
N(t — HOH (Oo) 
(t, x) Do’ 
i T 
¢(T)= | eae: At =¢(T)-f( = | eon, 
2 t 


where N (t, x) is the chosen numeraire. Note that in this case, the value 
of the numeraire is 1 today: N(0, 0) = 1. The last part of the model is 
the martingale valuation formula. Suppose at time ¢ the economy is in 
state X(t) = x. If V(t, x) is the value of any freely tradeable security, 
then (t, x) is a martingale. So in the LGM model we get the valuation 
formula 


Va, x) =ÑG, oe fee N R09 = |- N@x (VOX) iaag 
N(T,X 


/2n KE J N(T,X) 


If the security has intermediate cash payments, then we need to 
modify this formula appropriately. The prices in LGM model can be 
written in terms of the relative prices 


2. 


va, lE» 
Mx) 


Since the value of the numeraire is 1 today, values and relative values 
are equal today, i.e. vO, 0) = V(0,0) As we shall see, we only have to 
calculate the relative prices V(t, x) and never have to calculate the full 
prices V(t, x), This simplifies our formulas substantially. In terms of the 
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relative prices V(t, x), the LGM model is 


V(t, x) V(T, Xj 2 IAS ax, 


Jae | 


The value of a zero-coupon bond is 


Plt, X, T) 1 / 1 -1 (X-x)/2A¢ 
p(t, x, T) =~ = - e2 dX for anyT >t. 
M(t,x) S20 AT J N(T,X) 


Substituting for the numeraire and carrying out the integration yields 
the zero-coupon price 


p(t, x, T) = DT) HO 20s, 


At t = 0, the state variable is x = 0, by definition. Since p(0, 0, T) = D(T), 
the LGM model automatically matches today’s discount curve D(T). 

At t, x the instantaneous forward rate for maturity T, namely f (t, x, T), 
is defined via 


A T 7 7 
P(t, x, T) = p(t, x, DÂ, x, T) = eh FEXTAT, 


Similarly, the discount factor can be written in terms of today’s 
instantaneous forward rate fo(T) as 


D(T) = E So fo yar" 
This shows that for the LGM model, 
2 
Fx, T) = fT) + H' Tx + [AD] E). 


The last term [H’ (TEA is a small convexity correction; although it 
is needed for pricing, it does not affect the qualitative behaviour of the 
model. The other terms show that at any date t, the forward curve is 
made up of today’s forward curve fo(7) plus an amount x of the curve 
H'(T). The amount x of the shift is a Gaussian random variable with 
mean zero and variance ¢(f). 

The curve H’(T) is a model parameter; as we shall see, it replaces 
the mean reversion coefficient x(t) in the Hull-White CHW) model. The 
other model parameter is the variance ¢(t). It takes the place of the 
volatility o (t). As always, model parameters have to be set a priori dur- 
ing the calibration procedure by combining both theoretical reasoning 
(guessing) and calibration of vanilla instruments. 
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26.1.6 Hull-White 


We continue to study LGM. The HW model is given by 
dt = [O(t) —k(t)- r]dt+a(t)dW 
T 
p(t, ro; T) = E [ode ‘1 |r) = ro 


t 


For given x(t) and o(f) we adjust @(t) to match the yield curve; 
D(T) = p(0, ro; T). We can transform HW into a LGM model with 


t T T 
X(t) = exp [code rO -r0- f 00 - exp f eoar dt 
0 t 0 
Using It6 we can write 
dX = a(t)dW 
X(0) = 0, 


where X is a Gaussian martingale. In the LGM model we get the 
valuation formula 


V(t, x) V(T,X) V(T,X) 
=E xXQ=x|= p(T, Xvert, t, x)dX, 
N(t,x) N(T,X) N(T,X) 
where 
y(t, x] T,X) = —— eieo 
i ? ~V 2T AG 
N(t,x) = l HOROO 
f D(t) 
T T 
¿(T) = J a(t)dt; Ag=e(T)-¢(N= J a*(t)dt. 
0 t 


The model parameters are ¢(T), the accumulated volatility and H(T), 
the mean revision. For a zero-coupon bond (V(T, X) = 1) we have 


PIGXT) = Di) e BOROA 
N(t, x) 
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Table 26.1 The Hull-White and the LGM model functions 


Hull-White LGM 
a(t) = H’(t)/s'(t) rate volatility c accumulated volatility 
K(t) = -H“(t)/H'‘(t) reversion H(t) reversion 


A(t) = £6 (t) +k (t) - folt) + [H’(t)]?<(t) 


But 


T 
S F(tx, TOAT" 
p(t, x, T) = et 


gives 
F(t, x, T) = fo(T) + H'(T) -x+ H'(T)H(T) < o(0. 


So in terms of LGM we can write the HW model as 
dX = a (dW, X(0) =0 
f(t, x, T) =fo(T) + H'(T) -x+ H'(T) - HT): ®© 
r(t, x) = f(t, x, ) = PA + H'A) -x+ H'O - HO - cA). 


With the money market numeraire we have 


=F oa 
Vit,rop) =E |e? V(T, r(T))|r@ = ro 


dr = [O(t) —k(t)- r] dt + a(t)dw. 


To compare the models we use Table 26.1, where ¢(t) and H(t) are 
directly measured! 

The calibrated LGM prices are shown in Table 26.3 and their pre- 
dicted implied volatility in Table 26.4. The error in the At-The-Money 
volatility with diagonal = 2% are shown in Table 26.5. 

As we shall see, the value of any vanilla option depends only on the 
value of the variance at the exercise date, ¢ (tex), and on the mean re- 
version function H(t;) at the deal’s pay dates t;. Calibration determines 
the functions ¢(t) and H(t) fairly directly. Obtaining the mean rever- 
sion parameter x(t) requires differentiating H(t) twice, which is an 
inherently noisy procedure. Similarly, obtaining o(t) also requires dif- 
ferentiating ¢(D. This is why calibrating directly on the HW model 
parameters (instead of the LGM formulation of the model) is often an 
inherently unstable procedure. 
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The LGM parameters can be written in terms of the HW paramet- 
ers as 


t 
H =A | aidr 4 B 


0 
t 


1 t 
C(t) = af mes K(T)AT gr. 
0 


where A and B are arbitrary positive constants. Since different A and B 
yield the same HW model, and thus yield the identical prices, the LGM 
model has 2 invariant representations. 

First, all market prices remain unchanged if we change the model 
parameters by 


H(t) >C-H(Q), EA > iO 


for any positive constant C. To prove this, note that if we make the 
previous transformation and then transform the internal variables x 
and X by 


x > x/C,X > XIC, 


we obtain the same transition probabilities and zero-coupon bond 
prices that we started with. Second, all market prices remain un- 
changed if 


H(T) > H(T)+ KOT) > ¢(T) 


for any constant K. To prove this, note that if we make the previous 
transformation, and then transform the internal variables xand X by 


x—>x+Kt,X ~ X+KC, 


we obtain the same transition probabilities and zero-coupon bond 
prices as before. 

It is critical to pin down these invariances by arbitrarily choosing 
some value of H(t) and of ¢(t)) before calibration. Otherwise conver- 
gence would be infinitely slow, with numerical round-off determining 
which of the equivalent sets of model parameters is chosen. 
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26.1.6.1 Scaling 


On average, interest rates in G7 countries change by +80 bps or so 
over the course of a year. Equivalently, the standard deviation of H'(T)x 
should be about 1% or less each year. We choose to use the time scale 
of years (so AT = 1 means an elapsed time of 1 year) and we scale 
H(T) and H'(T) to be O(1). Then x is of order O(1% x yt at date t, 
and ¢(f) is of the order of O(10~“t). More precisely, suppose we have 
chosen H(0) = 0, and have scaled H(T) so it increases by 1 or so every 
year. Then 


H(T) ~O(),H(T) ~O(T),¢c(t) ~ 00.64 x 10t 
x, X ~ (0.8 x 10-7) /t 
H'(T)H(T)c(t) ~ (064 x 10+)T. 


26.1.7 Summary of the LGM Model 


The complete LGM model can be summarized as 


1 
V(t, x) = e2XVPAC AX foranyT <t 


1 
J2m At 
with A¢ = ¢(T) - ¢(t), and with the (reduced) zero-coupon bond 
formula being 


p(t, x, T) = DT) POOO 


and with x = 0 at t = 0. Consequently, ¢(0) = 0. 

These equations are the only facts about the model we need to price 
any security. This model automatically reproduces the discount curve 
D(T). The functions H(T) and ¢(t) are model parameters, which are 
set during the calibration step, where the model prices are matched 
to the market prices of selected vanilla instruments, usually caplets 
and swaptions. Once the model is calibrated, H(T) and ¢(t) are known 
functions, and the price of exotic deals can be determined from the 
aforementioned martingale formula, using the previous zero-coupon 
formula to calculate the payoffs. Later we will present the calibration 
and pricing steps in exquisite detail. 
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26.1.8 Calibration 
26.1.8.1 Calibration and Hedging 


Model calibration is the most critical step in pricing. It determines 
not only the price, obtained for an exotic deal, but also the hedges 
of the exotic. To see this, suppose we have some model M. It in- 
variably contains unknown mathematical parameters which are set by 
calibration. 

To calibrate, one selects a set of vanilla instruments whose volatil- 
ities (prices) are known from market quotes. Let these volatilities be 
O1,02,...,Op. 

The calibration procedure picks the model parameters by matching 
the model’s yield curve to today’s discount factors D(T) and matching 
the model’s price of the selected vanilla instruments to their market 
volatilities, either exactly or in a least squares sense. 

The calibrated model M’ is a function of today’s discount factors 
and these n volatilities. The calibrated model is then used to price the 
exotic deal. The only step in this procedure which uses market inform- 
ation is the calibration step. This means that the price of the deal is a 
function of today’s yield curve D(T) and the n volatilities 01, 02, . . . , On. 
The price of the deal depends on no other market information. 

Consider what happens at the nightly mark-to-market. The model 
is calibrated and deals are priced as mentioned before. Next the Vega 
risks are calculated by bumping the volatilities in the volatility mat- 
rix (cube) one by one. After each bump, the deal is priced using the 
identical software, and the difference between the new price and the 
base price is the bucket Vega risk for the bumped volatility. Unless 
the bumped volatility is one of the n volatilities used in calibration, 
it has no effect on the calibration of the model, so it does not af- 
fect the price of the instrument. An exotic deal only has Vega risks 
to the n vanilla instruments used in calibration. After the Vega risks 
are calculated for all the deals on the books, enough of each vanilla in- 
strument is bought/sold to neutralize the corresponding bucket Vega 
risk. This means that in the normal course of events, an exotic deal 
will be hedged by a linear combination of the vanilla instruments used 
during calibration. If the span of the vanilla instruments provides a 
good representation of the exotic, then the hedges should exhibit rock 
solid stability, with the day-to-day amounts of the hedges changing 
only as much as necessary to account for the actual changes in the 
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market place. If the vanilla instruments do not provide a good repres- 
entation of the exotic, then the hedges may exhibit instabilities, with 
day-to-day amounts of the hedges changing substantially even for rel- 
atively minor market changes. This latter is highly undesirable as the 
increased hedging costs gradually eliminate any initial profit from the 
exotic. (The nice term for this is leaking away your P&L.) Indeed, in 
practice, even small improvements in the algorithm for matching the 
hedges to the exotics pay off disproportionately in the adroitness of 
the hedging. 

One can move most of, and possibly all, the Vega risk from the 
calibration instruments to a different (and presumably better) set of 
hedging instruments by using risk migration. This is also known as ap- 
plying an external adjuster. As a side benefit, this method also improves 
the pricing, often dramatically. This technique will also be discussed 
as it pertains to the different exotics. 


26.1.9 Exact Formulas for Swaption and Caplet Pricing 


Under the LGM model, the reduced value of the swap is 


n 


Vrec(ts2) = J o (RE — Si) (t,x, t) + Px, tn) = p(x, 10) 
i=] 


n 
Vpay(t, xX) = p(t, x, 10) = P(t, X, tn) — > ov; (r - si) p(t, x, ti). 
i=l 


where 
p(t, x, ti) = Die Hao 


is the reduced value of the zero-coupon bonds. Here D; = D(ti), Hi = 
H(t;) are the discount factors and values of H(t) at the swap’s pay 
dates fj. 

Under the LGM model, the prices of vanilla swaptions, caplets, and 
floorlets depend on ¢(f) only through ¢ (tex), its value at the notification 
date. The swaption prices depend on H(T) only through the differ- 
ences H(t;) — H(to) for the pay dates t; of the fixed leg. This will be the 
key to creating lightning fast and stable calibration schemes. 

Under the one-factor LGM model, the exact pricing formulas for 
swaptions are 
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n 


Vee, = Y a; (a E sı) D:N (c + [Hi - Ho] 2 


= Cex 


yt [Hn — Ho] fe) 
+ D N | —————_ DoN 
( k Vv = 


. y + [Hi — Ho] Cex 
Vpay(0, 0) = DoN T= Ye i (R* - Si) Di a ( Cex ) 


( yt [Hn — Ho] ss) 
Cex l 


- D, N 


n 


= Vice, 0)+Do- ys Qj (Rf = Si) Di- Dn 
i=1 


Here y is obtained by solving 


n 
Ya, (r 2 si) Dje~HiHo)y-3(Hi-Ho)” Sex 4. p) e- Hn-Ho)y-3 Hn HoP Sex = Dy, 


i=l 
Newton’s method requires the derivatives of the prices with respect 
to the model parameters. We observe 


ð ð 
gare 0) = zH, ro. 0) 


E jeath (R - S;) Din’ (= [Hi — Ho] s) 


Sex 


oH 
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ð 
3H Jayr 0)= 5H, o. 0) 


. +[H, - H, 
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ð ə 
— Ý 0,0 ay 0,0 
0 Je rec( ) = Je pay( ) 
H; - H, : 
25 ow (R — Si) Di [H; — Ho] N' (= i 0] s) 
i=1 Sex 


+ D, [H - Ho] N' (= [Hn — Ho] =) 
n lAn . 


Sex 
The caplet/floorlet prices are given by 


Veapter(0, 0) = DoN(di)- (1 + & [R° — sı |) DiN) 


1+a[R° si] j 2 
log 1+a[R™ -s1 ] + (Ey Ho) Cex 


d (Hi Ho) ia 
and 
Vioorter(0, 0) = (1 +ã [R _ sı]) DIN(d*) - DoN (dž) 
SP nia; 0) + (1 +ã [R -si) D\-Do 
_ log Faery E Gh -Hotar 
aia (Hi — Ho) Vex 
where 


o _ D(t)-D() 
= ã&D) 


is the forward FRA rate. The caplet/floorlet prices are clearly Black’s 
formulas for call/put prices for an asset with forward value Do, strike 


(1+ aR” qpP1> and implied volatility satisfying 


Simp tex = (H| — Ho) V tex 
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26.1.10 Approximation of Vanilla Pricing Formulas for the 
One-Factor LGM Model 


It is practically useful to develop approximate formulas for the one- 
factor LGM model, even though exact closed form formulas are 
known. Recall that market prices for swaptions are usually quoted in 
terms of Black’s formula 


rec 


v0) = [R Na) - R?N(A9)} Lo 


mkt 


po) = [RON(=d9) - RNN) Lo = Vre (0) -Lo {R -R° } 


pay rec 
where 
Rie 1 12 
d l So + 50 tex 
n o fex 
and 


_ Do -D1 + Vie GSiDi 

Lo l 
Here D; = D(t;) are today’s discount factors at the pay dates. By using 
equivalent volatility techniques (or direct asymptotic), one discov- 


ers that under the LGM model, the implied (Black) volatility of the 
swaption is approximately 


Vta YL, a; (R? — Si) Di (Hi — Ho) + Dn (Hn - Ho) 


/JRFERo ie &Di 


This provides a good way to use market quotes of the implied volat- 
ility to obtain initial guesses for calibration. One can rewrite these 
quotes more simply in terms of the implied normal volatility. Under 
the Gaussian (normal) swap rate model, the value of the swaption is 


R Ri _ RO Ri _ R? 
viet) = {cat —R°)N (5) - ON JTexG (~~) | Lo 


R? 


OB tex © 


ON 4/ Tex ON4/ Tex 

P R? _ R R? — R&® 
O fæ - R®®)N (==) - on J/TexG (=) | Lo 

ý ON 4/ Tex ON4/ Tex 
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The equivalent volatility work shows that the implied normal or 
absolute volatility oy, is approximately 


Vta Xi- ai (R? - Si) Di (Hi — Ho) + Dn (Hn — Ho) 
~ RhxRO pera: l 


OBN tex © 


26.1.10.1 Forward Volatility 


Forward volatility is a key concern of calibrated models. Suppose that 
we calibrate a model and then ask what the swaption volatilities will 
look like at a future date f in the future. If the volatilities are increasing 
with t, we may be buying future volatility at too dear a price, and if 
volatilities are decreasing with t, we may be selling future volatility too 
cheaply. 

If we repeat the previous equivalent volatility analysis at a date ¢ in 
the future, then we discover that the (normal) swaption volatility at 
that date is 


VE Cex) EO Ly a (R° - Si) Di (H) — H(to)) + Da (H (tn) — H(to)) 
Tex —¢ Xei &iDi 


On~ 


where 


_ Do- Dit Vie GSiDi 
Lo 


R? 


If H(t) is decreasing exponentially, then ¢ (Tex) should be increasing 
exponentially to compensate. 


26.1.10.2 Calibration Strategy 


The most critical aspect of pricing is choosing the right set of vanilla 
instruments for calibrating the model. Even small improvements in 
matching the vanilla instruments to the exotic deals often lead to sig- 
nificant improvements in the price and the stability of the hedge. For 
each type of exotic, the best calibration strategy often cannot be de- 
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termined from purely theoretical considerations. Instead, one needs 
to determine which method leads to the best (the most “market fit”) 
prices and risks. 

We will now briefly discuss calibration strategies, illustrating the 
different strategies with a simple Bermuda swap (Bermudan swap- 
tions and callable swaps are considered much more carefully in a later 
section). 

Consider a Bermudan receiver with start date tọ, end date t,, and 
strike R™. Let the fixed leg dates be fo, f,...,fn, and let the exer- 
cise dates be 1, T2, . . . , Tn. If the Bermudan is exercised at 1;, then the 
holder receives the fixed leg payments 


ajR®™ paid at t,i=j,j+1,...,n—-1 
1 +a, R paid at tņ, 


where 


aj = cv8g(ti-1, ti, P) 


is the coverage (day count fraction) for interval i computed according 
to the appropriate day count basis 6. In return, the holder makes the 
floating leg payments, which are worth the same as 


1 paid at t;-1 
aS; paid at t; i =j,j+1,...,m-— 1, 


where 
æ! = cv8g(Ti-1, Ti, P^. 


Here we have adjusted the basis spread to the fixed leg’s frequency 
and day count basis as discussed earlier. Therefore, if the Bermudan is 
exercised at Tj, one receives/makes the payments 


-1 at tı 
a; (R — S;) at tfori=j,j+1,...,n— 1 
1+, (R = Sn) at iy 
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Clearly at any point t, x the jh payoff is worth 


n 


VF aa) = J o (R* - Si) Px t) +P, tn) -Px fe). 


i=j 


The first step in calibration is to characterize the exotic, extracting its 
essential features.»> If the Bermudan is exercised on exercise date Tj, 
one receives a swap worth 


n 


F(t) = J o (R -S;) 1, x, t) +x, tn) -EC x, tj) 


i=j 


The first step in calibration is to characterize the exotic, extracting its 
essential features. If the Bermudan is exercised on exercise date Tj, one 
receives a swap worth 


n 
X a(R — SPCT, x, ti) +P, x, tn) -PC x, tj-1) at T. 
i=j 


Suppose we evaluate this swap using today’s yield curve with a parallel 
shift of size y, 


P(t, x, ti) > Dte" = Die". 


The shift y; at which the j’ h swap is at the money is found by solving 


Fa; (R — S1) Die MH 4 Dye MOHD = Dp, 
i=j 

The Bermudan is characterized by 

(a) the set of exercise dates T1, T2,..., Tn; 

(b) the set of parallel shifts y; for j = 1,2...,n; and 

(c) the length t, — to of the longest swap. 
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The second step is to select a calibration strategy and choose the calib- 
ration instruments. As we shall see, since we have 2 functions of time 
to calibrate, we can calibrate 2 separate series of vanilla instruments. 

Under the LGM model, the prices of vanilla swaptions, caplets, and 
floorlets depend on ¢ (t) only through ¢ (tex), its value at the notification 
date. The swaption prices depend on H(T) only through the differ- 
ences H(t;) — H(to) for the pay dates t; of the fixed leg. This will be the 
key to creating lightning fast, stable calibration schemes. 

The trick is to calibrate on vanilla instruments whose pay dates line 
up exactly. We now go through the various calibration strategies for 
this Bermudan. 


26.1.11  swaptions 


From 
Vax) _ VT, X) oh if VEO 
T E K X |X(t) = J =] NTD (T, X|t, x)dX 
we have today 
_ f VT,X) 
V(0,0) = NCT, xÝ Tx, 0)dX. 


The payoff at tex is given by (Fig. 26.1) 


F(T) + xE (T) 


Fig. 26.1 Cash flows for a swaption 
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n 


Ry », atjP(texs X, tj) + D(tex, X; tn) — P(tex, X, to) 
j=l 
Plte x, T) _ D( ppe POIO Se, 
N (tex, X) 


If we integrate this for a receiver swaption we get 


n 
x-Hj -=H -H 
Viec = Rf Ý ajD;N (es s) + Dn N (s) -DoN (s) 
` ` N Sex sN Sex Sex 


j=l 


where x solves 


n yl 2 
Ry D ajDje” 2 Sex + Dyetn® 2H Sex = Doeltox-2Ho Sex 
FA 
Sex = G (tex), Hi; = H(t), Ho = H(to) 


As we can see, the price only depends on the accumulated volatility at 
the expiry date ¢ (tex) and on the mean reversion on pay days H(t;),j = 
0,1,...,n. From a “put/call” parity, we also have the price of a payer 
swaption 


n 
Voay = Vree = Ry be ajDj — Dy + Do. 
jel 


26.1.12 Bermudan Swaption 


We study a “2 into 4, annual pay, receiver at 8%, 5 days’ notice”. The 
fixed leg is given by 


Raj paidat tj;=1,2,...,n—1 
1+Rf&n paidat tn, 
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where Rẹ is the fixed rate and with notification dates qj, j = 0,1,...,n- 
1. If we exercise at t; we get a swap starting at t; with fixed leg pay- 
ments at 41, t)42,...,t, and a floating leg at t; worth 1. The current 
value is given by 
n 
V=Rp D> appt, X, te) + P(t, tn) -plj X, tj). 
j=1+1 


26.1.13 Calibration, Diagonal + Constant x 


In the calibration process we choose the model parameters ¢(t) and 
H(t) to match the LGM prices to the market. For a “2 into 4 Bermudan 
receiver @ 8%” we use for the underlying risk: “2 into 4”, “3 into 3”, “4 
into 2” and “5 into 1”, i.e. the diagonal. For the smile risk we calibrate 
to swaptions @ 8% (Table 26.2). 


k(t) =-H" (H(t) => H(t)=A+B-e*! 
H(t) > H(A +C 
H(t) > H(p/eC, cA > C. ct) 


Set H(t) = (1 —e*")/k Then calibrate 


2 into 4 ¢ (t2), H(to), H(t3),..., H(t) => C(t2) 
3 into 3 ¢ (t3), H(t3), H(t4), . . ., H(t6) => C(t3) 
á into 2 ¢(t4), H(t4), H(ts), H(t) => ¢(t4) 
5 into 1 ¢(ts), H(ts), H(t) => C(t5). 


Then use ¢(0) = 0 and the above and use linear interpolation. 
Choose « = 2%, which gives H(t) and then use the diagonal to get ¢;. 
The calibrated LGM prices are shown in Table 26.3 and their predicted 
implied volatility in Table 26.4. The error in the At-The-Money volatility 
with diagonal k = 2% are shown in Table 26.5. 
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Table 26.3 LGM prices 
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ly 


1y 0.58% 
2y 0.88% 
3y 0.89% 
Ay 0.89% 
5y 0.83% 
6y 0.81% 
7y 0.78% 
8y 0.75% 
9y 0.72% 
10y 0.66% 


2y 

1.62% 
1.68% 
1.71% 
1.65% 
1.59% 
1.55% 
1.49% 
1.44% 
1.35% 
1.26% 


3y 


2.33% 
2.42% 
2.40% 
2.35% 
2.28% 
2.22% 
2.14% 
2.03% 
1.92% 
1.81% 


4y 5y 

2.97% 3.48% 
3.03% 3.60% 
3.03% 3.61% 
2.99% 3.57% 
2.91% 3.48% 
2.84% 3.37% 
2.70% 3.20% 
2.57% 3.05% 
2.44% 2.90% 
2.30% 2.74% 


6y 


3.96% 
4.11% 
4.13% 
4.10% 
3.97% 
3.84% 
3.65% 
3.49% 
3.33% 
3.19% 


7y 


4.39% 
4.57% 
4.61% 
4.54% 
4.40% 
4.26% 
4.05% 
3.89% 
3.74% 
3.59% 


8y 


4.77% 
4.99% 
5.01% 
4.93% 
4.80% 
4.64% 
4.43% 
4.37% 
4.12% 
3.96% 


9y 

5.13% 
5.33% 
5.35% 
5.28% 
5.14% 
4.99% 
4.78% 
4.63% 
4.47% 
4.32% 


10y 


5.93% 
5.61% 
5.65% 
5.59% 
5.46% 
5.32% 
5.11% 
4.96% 
4.81% 
4.66% 


Table 26.4 Implied ATM Volatilities (from LGM prices) 


1y 


1y 14.65% 
2y 14.67% 
3y 14.53% 
4y 14.40% 
5y 14.12% 
6y 13.83% 
7y 13.47% 
8y 13.27% 
9y 13.11% 
10y 12.97% 


2y 

14.46% 
14.53% 
14.38% 
14.20% 
13.95% 
13.70% 
13.32% 
13.15% 
13.00% 
12.94% 


3y 4y 5y 


14.30% 14.15% 13.99% 
14.38% 14.21% 14.05% 
14.19% 14.02% 13.87% 


14.02% 
13.81% 
13.55% 
13.19% 
13.04% 
12.94% 
12.91% 


13.87% 13.72% 
13.66% 13.53% 
13.42% 13.30% 
13.07% 12.98% 
12.96% 12.91% 
12.90% 12.86% 
12.88% 12.85% 


6y 


13.83% 
13.90% 
13.72% 
13.59% 
13.41% 
13.20% 
12.92% 
12.86% 
12.82% 
12.83% 


7y 


13.69% 
13.76% 
13.59% 
13.46% 
13.30% 
13.13% 
12.86% 
12.81% 
12.80% 
12.82% 


8y 


9y 


13.55% 13.42% 
13.63% 13.50% 
13.46% 13.35% 
13.35% 13.26% 
13.22% 13.15% 
13.06% 13.00% 
12.81% 2.77% 
12.78% 12.76% 
12.78% 12.76% 
12.81% 12.80% 


10y 


13.30% 
13.39% 
13.26% 
13.18% 
13.08% 
12.95% 
12.73% 
12.74% 
12.75% 
12.79% 


Table 26.5 Error in ATM volatility (diagonal, kappa = 2%) 


1y 


1y 1.03% 
2y 0.31% 
3y 0.35% 
4y 0.30% 
5y 0.03% 
6y 0.00% 
7y -0.11% 
8y -0.07% 
9y 0.00% 
10y 0.10% 


2y 
0.81% 
0.36% 
0.24% 
0.10% 

-0.04% 
0.01% 

-0.06% 
0.00% 
0.08% 
0.25% 


3y 


0.65% 
0.20% 
0.16% 
0.06% 
-0.05% 
0.03% 
0.00% 
0.08% 
0.22% 
0.42% 


4y 5y 


0.59% 0.40% 
0.10% 0.27% 
0.02% 0.11% 
-0.06% 0.02% 
-0.11% 0.00% 
0.00% 0.08% 
0.00% 0.08% 
0.14% 0.30% 
0.32% 0.53% 
0.54% 0.81% 


6y 


0.32% 
0.16% 
0.05% 
0.00% 
0.00% 
0.09% 
0.09% 
0.33% 
0.59% 
0.89% 


7y 


0.26% 
0.05% 
0.00% 
-0.02% 
0.03% 
0.11% 
0.11% 
0.36% 
0.65% 
0.97% 


8y 


0.13% 
0.00% 
0.03% 
0.01% 
0.06% 
0.14% 
0.14% 
0.44% 
0.77% 
1.13% 


9y 
0.00% 

-0.05% 

-0.05% 
0.01% 
0.10% 
0.18% 
0.19% 
0.53% 
0.90% 
1.29% 


10y 


-0.12% 
-0.08% 
-0.04% 
0.04% 
0.15% 
0.24% 
0.24% 
0.62% 
1.02% 
1.44% 
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26.1.14 Calibration to the Diagonal with H(T) Specified 


Suppose that H(T) is specified a priori. (A possible source of such 
curves H(T) is indicated next.) Typically H(T) is given at discrete 
points H(T|), H(72),...,H(Ty). In that case, piecewise linear interpol- 
ation is used between nodes. This is equivalent to assuming that all 
shifts of the forward rate curve are piecewise constant curves. 

With H(T) set, we can use the preceding procedure and formulas to 
calibrate on the diagonal swaptions. This determines the value of ¢(f) 
at T],T2,..-,T,. AS before, one adds the point ¢(0) = 0, ensures that 
the ¢; = ¢(7;) are increasing, and uses piecewise linear interpolation to 
obtain ¢(f) at other values of t. 

Origin of the H(7). Suppose one had the set of Bermudan swap- 
tions 30 NC 20, 30 NC 15, 30 NC 10, 30 NC 5 and 30 NC 1. Wouldn't 
it be nice if the same curve H(T) were used for each of these Ber- 
mudans? The 30 NC 10 Bermudan includes the 30 NC 15 and the 30 
NC 20 Bermudans. It would be satisfying if our valuation procedure 
for the 30 NC 15 and 30 NC 20 assigned the same price to these Ber- 
mudans regardless of whether they were individual deals or part of a 
larger Bermudan. 

One could arrange this by first using a constant x, let’s call it k4, to 
calibrate and price the 30 NC 20 Bermudan. Without loss of generality, 
we could select 


H'(T) = ek4(T30-T) 
1 — ek4(230-T) 
A(T) = ————— 
K4 
for To9 < T < T39. We would calibrate on the diagonal to find ¢(t) 
at expiry dates Tm, Tn+1,... beyond 20 years, and then price the 30 
NC 20 Bermudan. Selecting the right value of «4 would match the 
Bermudan price to its market value. Neither the swaption prices nor 
the Bermudan prices depend on A(T) or ¢(t) for dates before the 20- 
year point. 
To price the 30 NC 15, one could use the H(T) obtained from «4 for 
years 20 to 30, and choose a different kappa, say «3, for years 15 to 20: 
H'(T) = e3 (T2077) ok4(T30-T) 
1 — e3 T207) 


H(T) = e“4(T30-T20) 4. 
K3 K4 


1 — e“4(T30-T20) 
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Calibrating would produce the same ¢ (t) values for years 20 to 30 as be- 
fore. In addition, for each «3 it would determine ¢(f) for years 15 to 20. 
By selecting the right «3, one could match the 30 NC 15 Bermudan’s 
market price. 

Continuing in this way, one produces the values of ¢(t) and H(T) for 
years 10 to 15, for years 5 to 10, and finally for years 1 to 5. These 
¢(t) and A(T) would then yield a model which matches all the di- 
agonal swaptions and happens to correctly price all the liquid, 30y 
co-terminal Bermudans. These «(tys turn out to be extremely stable, 
only varying very rarely, and then by small amounts. Typically a desk 
would remember the «(f)’s as a function of the co-terminal points, 
relying on the same x(f)’s for years. 

In general, if T„ is the co-terminal point and Tọ, T),...,7,_1 are the 
“no call” points, then A(T) is 


— e&i(Tj-T) — eka) 2 
na = [] etto y y LEA CAA [] ett 
i=j+1 k=j+1 i=k+1 


for Ty S <T< Tj. 
After H(T) and ¢(t) have been found, one can use the invariants to 
re-scale them if desired. 


26.1.15 Calibration, Diagonal + Linear ¢(t) 


This is an idea pioneered by Solomon brothers. Let us use a constant 
local volatility œ. Then 


t 


C(t) = J adu = at 
0 
is linear. By using the invariance ¢ (t) > ¢(t)/ C?,H (T) > CH(T) we can 
choose « to be any arbitrary constant without affecting any prices. So 
we choose 


C(t) = agt, 


where f is measured in years, and the dimensionless constant ao is 
typically 10°. For this calibration, we use the other invariant to set 
H, = H(t,) = 0. We now determine the values of H; for other values 
of i by calibrating on the diagonal swaptions, starting with the last 
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swaption. Recall that the price of the jth diagonal swaption is 


; n * + [Hi -Hjls 
9,04 (0,0) = R° D o; (RM — Si) DiN ( —— f 
. = JG 


yi + [An — Aj-als; yj 
+D,N | 1—4 ) -pa ( 2 


and its derivative with respect to Hj_; is 


* 4 [H -Hlc 
— J/GDn|Hn- Hj- |G ieies) 
Ta 


Here y; is given implicitly by 


n 
YO a(R, )Dje PA] Na Sppe D a Ar MS = D 
i=j 
Consider the last swaption, j = n. It depends on ¢, = ¢Ẹ(Tn), on Ay, 
and H„—1. Of these, ¢, is known, H, has been set to zero, so only H-1 
is unknown. Since ymod is a decreasing function of H„-1, there is a 
unique value of H,_; which matches the model price to the market 
price. This can be found easily using a global Newton’s method. We 
can then move onto the j = n — 1 swaption. This swaption depends 
on Hy_2, which is unknown, and ¢;-1, Hn-1, and Hn, which are known. 
Working backwards like this, we can calibrate all of the swaptions, and 
for each calibration there will only be a single unknown parameter, 


Hy}. 

J 

This calibration procedure will yield Ho, M,...,H, on the dates 
to,t1,...,t,. One uses linear interpolation/extrapolation to get H(t) at 


other values of t. 
Infeasible values. In deriving the swaption formulas, we assumed 
that H(T) was an increasing function of T. (This assumption was 
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stronger than we needed: inspection of the previous argument shows 
that one only needs to assume that there is a unique break-even point 
y* , with in-the-moneyness on the left.) Since we are calibrating the H; 
’s separately, it may happen that H,_; may exceed H; . An practice, this 
has never happened to my knowledge. Still one must be prepared.) 
After each Hy_; is found, one should check to see that Hj_) < Hj 

If this condition is violated, one should reset H,_| = H; . This means 
the jth swaption would not match its market price exactly. Instead it 
would be the closest feasible price. 

This kind of calibration is used by Solomon & Brothers. The hedging 
instruments for the underlying risk are again “2 into 4”, “3 into 3”, “4 
into 2” and “5 into 1”, i.e. the diagonal. 

We then assume that ¢(f) is linear 

s(t) = at set æo = 107 


H(t) > H(t)+C 
H(t) > HOC, c(t) > C. cA) 


Set H(t¢) = 0 Then calibrate 


5 into 1 ¢(ts), H(ts), H(t) => H(t2) 
4 into 2 ¢(t4), H(t4), H(ts), H(t) => H(t3) 
3 into 3 C(t3), H(t3), H(t4),...,H (te) => H(t4) 
2 into 4 C(t2), H(t2), H(t3),...,H (te) => H(ts). 


Use a global Newton-Raphson method and linear interpolation to find 
H(t ti). 
The error after calibration are shown in Table 26.6. 


Table 26.6 Error in ATM volatility (diagonal, constant alpha) 
ly 2y 3y 4y 5y 6y 7y 8y 9y 10y 


1y -1.29% -0.52% -0.35% -0.13% -0.08% 0.06% 0.10% 0.05% 0.00% -0.06% 
2y -0.39% -0.33% -0.32% -0.24% 0.11% 0.04% -0.01% 0.00% 0.00% 0.00% 
By -0.46% -0.32% -0.18% -0.11% 0.01% 0.00% 0.00% 0.01% 0.02% 0.05% 
Ay -0.16% -0.16% 0.01% -0.13% -0.01% 0.00% 0.01% 0.03% 0.07% 0.11% 
5y -0.13% 0.02% -0.08% -0.14% 0.00% 0.02% 0.06% 0.10% 0.14% 0.20% 
6y 0.24% 0.00% 0.00% 0.00% 0.10% 0.11% 0.14% 0.18% 0.22% 0.28% 
7y -0.29% -0.14% 0.00% 0.04% 0.15% 0.17% 0.19% 0.24% 0.29% 0.34% 
8y -0.14% 0.00% 0.10% 0.17% 0.34% 0.37% 0.41% 0.50% 0.58% 0.68% 
9y 0.00% 0.08% 0.22% 0.32% 0.53% 0.59% 0.65% 0.77% 0.90% 1.02% 
10y 0.04% 0.19% 0.36% 0.48% 0.75% 0.83% 0.91% 1.07% 1.23% 1.35% 
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26.1.16 Calibration, Diagonal + Row 


The hedging instruments for the underlying risk is again “2 into 4”, “3 
into 3”, “4 into 2” and “5 into 1”, i.e. the diagonal. And exercise or 
wait (forward volatility risk), use “2 into 1”, “2 into 2”, “2 into 3” and 
“2 into 4” (i.e. a row). 

Set ç(t2) = aĝt2H(t2) = 0 and set a = 107 


A(t) > H(t)+C 


Then calibrate 
2 into 1 C(t2), H(t2), H(t3) 
2 into 2 ¢ (t2), H(t2), H(t3), H(t4) 
2 into 3 C(t2), H(t2), H(t3),..., HA (ts) => A (ts) 
2 into 4 C(t2), H(t2), H(t3),..., H(te) => H (te) 


H(t) > H(tV/c, 


c(t) > C- s(t) 


=> H(t3) 
=> H(t4) 


3 into 3 ¢(t3),H(t3), H(t4),...,H (te) => (ts) 
4 into 2 €(t4), H(t4), H(ts), H(t) 
5 into 1 (ts), H(ts), H(t6) 

Use a global Newton-Raphson method and linear interpolation. As 
we will see next, over-calibrating gives worse result. Similar result will 
be found if we use a column instead of a row. 

The error after calibration are shown in Table 26.7. 


Table 26.7 Error in ATM volatility (diagonal and row) 


=> C(t4) 


=> 


C(ts). 


ly 
1y 0.00% 
2y -0.47% 
3y -0.27% 
4y -0.56% 
5y -0.13% 
6y -0.58% 
7y -0.81% 
8y -0.20% 
9y 0.00% 
10y -0.01% 


2y 
0.00% 
-0.27% 
-0.43% 
-0.31% 
-0.38% 
-0.53% 
-0.30% 
0.00% 
0.08% 
0.14% 


3y 
0.00% 
-0.41% 
-0.23% 
-0.33% 
-0.42% 
-0.19% 
0.00% 
0.12% 
0.22% 
0.32% 


4y 
0.00% 
-0.30% 
-0.33% 
-0.41% 
-0.28% 
0.00% 
0.12% 
0.20% 
0.32% 
0.43% 


5y 
0.00% 
-0.06% 
-0.20% 
-0.14% 
0.00% 
0.21% 
0.27% 
0.37% 
0.53% 
0.70% 


6y 
0.00% 

-0.13% 

-0.10% 
0.00% 
0.11% 
0.30% 
0.33% 
0.41% 
0.59% 
0.78% 


7y 
0.00% 

-0.09% 
0.00% 
0.10% 
0.22% 
0.39% 
0.37% 
0.45% 
0.65% 
0.86% 


8y 

0.00% 
0.00% 
0.08% 
0.19% 
0.31% 
0.46% 
0.43% 
0.54% 
0.77% 
1.02% 


9y 

0.00% 
0.06% 
0.15% 
0.27% 
0.39% 
0.54% 
0.49% 
0.63% 
0.90% 
1.18% 


10y 


-0.02% 
0.11% 
0.22% 
0.35% 
0.47% 
0.62% 
0.56% 
0.72% 
1.02% 
1.34% 


26.1.17 Calibration, 


Caplets + Constant x 


In the calibration process we choose the model parameters ¢ (t) and 
H(t) to match the LGM prices to the market. For a “2 into 4 Bermudan 
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receiver @ 8%” we use for the underlying risk “2 into 1”, “3 into 1”, “4 
into 1” and “5 into 1”, i.e. a column. 


k(t) =-H"()/H'(t) => H(t)=A+B-e* 
H(t) > H(t)+C 
H(t) > H(IC, cA —> C. cA 


Set H(t) = (1 — e™*)/k Then calibrate 
2 into 1 ¢(t2), H(t2), H(t3) => ¢(t2) 
3 into 1 ¢(t3), H(t3), H(t4) => ¢(t3) 
á into 1 ¢(t4), H(t4), H(t5) => ¢(t4) 
5 into 1 ¿ (ts), H(ts), H(te) => ¢ (ts). 


Then use ¢(0) = 0 and the above and linear interpolation. You can 
also use caplets + linear ¢ (t). 


26.1.18 Calibration to Diagonals with Prescribed ¿œ 


Suppose ¢ (t) is a known function which is increasing and has ¢ (0) = 0. 
We could carry out the preceding calibration procedure to determine 
H(T) from the diagonal swaptions; the procedure does not depend on 
C(t) being linear. 


26.1.19 Calibration to Diagonal Swaptions and Caplets 


A Bermudan swaption can be viewed as the most expensive of its com- 
ponent European swaptions, plus an option to “switch” to a different 
swaption should market conditions change. The component swap- 
tions are just the diagonal swaptions, so calibrating to the diagonals 
accounts for this part of the pricing. On any exercise date, “switch” 
option is the option to exercise immediately, or to delay the exercise 
decision until the next exercise date. Since these delays are short, typ- 
ically 6 months, one may believe that the switch option can best be 
represented by short underlyings. Accordingly, one could argue that 
one should calibrate to either a column of caplets or a column of 1 
year underlyings, as well as the diagonal swaptions. Here we calibrate 
on the caplets and swaptions simultaneously; in the next section we 
calibrate to the diagonal swaptions and the swaptions with 1 year 
underlyings. 
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26.1.20 Calibration to Diagonal Swaptions 
and a Column of Swaptions 


One could argue that caplet and swaption markets have distinct iden- 
tities, and that mixing the 2 markets introduces small, but needless, 
noise. Instead one could calibrate on the diagonal swaptions and a 
column of swaptions with 1 year tenors. dn most currencies, these 
are the swaptions with the shortest underlying available.) 


26.1.21 Other Calibration Strategies 


There are many other simple calibration strategies; although they are 
not overly appropriate for pricing a Bermudan, they may well be 
appropriate for other deal types. 


26.1.21.1 Calibrate on Swaptions with Constant x 
or Specified H(T) 


Suppose we have chosen a constant mean reversion parameter «K, 
or have otherwise specified H(T). Then the calibration procedure 
just needs to find ¢(t). Suppose we have selected an arbitrary set of 
n swaptions to be our calibration instruments. In LGM valuation of 
each swaption the only unknown parameter is ¢(t) at the swaption’s 
exercise date. Using a global Newton’s method to calibrate each swap- 
tion to its market value thus determines ¢(t) and the exercise dates 
T1, T2,...,Tn Of the n swaptions. After obtaining the ¢j = C(t/), we 
need to ensure that ¢(tj) are non-decreasing, altering the offending 
values if necessary. We then include the value ¢0 = ¢(0) = 0, and use 
piecewise linear interpolation to obtain ¢(t) at other dates. 

Note that this method fails if 2 swaptions share the same exercise 
date t; calibration would either yield the same ¢, in which case one 
of the swaptions is redundant, or differing ¢, in which case our data 
is contradictory. If the exercise dates of any 2 swaptions are too close, 
say within 1-2 months, the results may be problematic. For this reason 
one usually ensures that the swaption exercise dates are, say, at least 
21/, months apart, excluding instruments from the calibration set to 
achieve this spacing, if necessary. 


26 Pricing Exotic Instruments 713 
26.1.21.2 Calibrate on Swaptions with Specified ¢(t) 


Suppose we have chosen a linear ¢(f), or otherwise specified para- 
meter ¢ (t). The calibration procedure just needs to find H(T). Suppose 
we have selected an arbitrary set of n swaptions to be our calibration 
instruments. We can then arrange the swaptions in increasing order of 
their final pay dates. Let these final pay dates be T1, 72, ..., Tn. Suppose 
we use our invariance to set H = H(0) = 0, and we use piecewise linear 
interpolation 


H(T)=A\T for <T}, 
k-1 
H(T) = > A(T =Teg)t# he Tes) “for Ten TS Ty, 
i=l 
n-1 
H(T) = 3 AilTi-T1)+ An(T-T,1)  forT,<T 
i=1 


where Tp = 0. 

For the first swaption, the slope A; determines the value of H(T) at 
all the swaption’s pay dates. Since ¢ (t) is known, the LGM value of the 
swaption depends only on a single unknown quantity, A. It is easily 
seen that the value is an increasing function of Al, so one can use 
a global Newton scheme to find the unique Al which matches the 
swaption’s price to its market value. The value of H(T) at the second 
swaption’s pay dates is determined by both A; and Ao, of which only 
A2 is unknown at this stage. Again a global Newton scheme can be 
used to find the A» needed to calibrate the swaption to its market 
value. (In rare cases it may occur that Az < 0; in this case we need to 
set A> = 0, its minimum feasible value.) 

We then continue in this way, calibrating the swaptions and obtain- 
ing the Aj;’s in succession. This method will fail only if 2 deals have 
the same final pay date, and will work poorly if the final pay dates are 
too near together. For this reason one usually ensures that the final pay 
dates are, say, at least 2!/, months apart, excluding instruments from 
the calibration set to achieve this spacing, if necessary. 
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